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8 PARTICLE INTERACTIONS AND 
CENTRAL FORCES 


8.1 MULTI-PARTICLE CONSERVATION LAWS 

Multiparticle conservation theorems now. A generalization of results for momentum and energy 
reached in Ch.2 for a single particle. 


Latin indices (a,b,c, etc.): vector indices 
Greek indices (a,p,y, ...): particle number 


Reminder of conserved quantities: 


Invariance of L 
r a -> r a + 8r 

r« -> f a + (80 x f a ) 
t —» t + 5t 

Consider an arbitrary system of particles: 


conserved quantity 

P = Z m « r a 

a 

L = Z X Pa 

a 

H ( = T + U) 



Center of mass located by 


/ 

M 

V 


R 



i 


rrR 


) 


( 8 . 1 ) 


We have 


X m a r a = Z m a( r a ~ R) = MR - MR = 0 . 


( 8 . 2 ) 
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Particle a: 


; internal force on 
a due to P 

) 

: external force 

Total force on a: 

t = F« e> + X tp • (8.3) 

P 



Call f a = ^ f a p . From Newtons third law 

4p = -V (weak form) (8.4) 

Also assume 

4a = 0- (8.5) 


(Newtons mechanics not equipped to handle self-interactions which, however, really do exist!) 
Newtons second law: 


n = F (e) 4- f 

tr a a c 


(8.6) 


or 


_di 

dt 2 


( m a r a) 


n(e) 


+ 


X f p 


Sum on a: 


_dl 

dt 2 



X *£ e) + X f P 

a a ,p r 


X f ap - f 12 + tl, + ••• - 0 ' 

a,P ' 5 ' 

=> = Xr = H" e) . 

a 


(8.7) 


( 8 . 8 ) 
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Center of mass moves as if the total external force were acting on the entire mass of system concentrated 
at the center of mass. (Internal forces have no effect on CM motion.) Of course, if F (e) = 0, then 

(mr) = P = 0 , 

dt 2 v ' 

and total linear momentum is conserved. 


Likewise, the angular momentum of the a— particle about the origin is 


L = r x p . 

a a -tra 


( 8 . 10 ) 


Summing this: 


L = X L a = X r a X Pa = X fa X Va) 

a a a 

= 2 f« + *) X m « (^a + R) * 

a 

L=RxP + ^r a x p a . (p a = m a r a ) 


( 8 . 11 ) 

( 8 . 12 ) 


Total angular momentum about a coordinate axis is the angular momentum of the system as if it were 
concentrated at the center of mass, plus the angular momentum of motion about the center of mass. 
Now the rate of change of L a is 


L = r x p = r x 

a a a 


c + 




'aP 


L = X L a = X r a X + X r a X f a P * (8.13) 

(assume = 0 ) 


a, P 

© © 


- z 


© 

© 

Therefore, we may write 


x Fa"’ 


N (e) , 


Z r x f, 

a 

a, P 


= - X r p x f < 


a, P 


'aP 


(8.14) 

(8.13) 


® = \ X ( r « ■ r p) x f «p = 0 • (see fi 9 ure 


below) 


(8.16) 


a,p 
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I am showing an attractive force case. I am also assuming the “strong form” of Newtons second law, 
which says that f ap and fp a lie along the line connecting the two particles. (Violated, for example, 
in electromagnetism.) Under these circumstances then 


L = N (e) . 

So L = const, in time if N (e) = 0. 


(8.17) 


For the kinetic energy, we have 


T = Z \ m a r a 2 = \ Z m a ( R + r aJ 

= \ M ^ 2 + \ Z m a f a 2 • (8.18) 

Kinetic energy, like angular momentum, consists of two parts: the kinetic energy of the center of 
mass, plus the kinetic energy of motion about the center of mass. 

Now, what about work being done on the system? Define 

2 2 

W 12 s E J f a • = X 1 m a r a • dt 

a i a i 

= f «(^2a - Vi a ) / 

a 1 a 


W 


12 


T — 
± 2 




(useful in a bit) 


(8.19) 
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Assume conservative external and internal forces: 


Fa"’ = -V a U a (Xj, 


^a Vtp p j 


Then, alternatively 


© 


© 


© 


w, 


12 


= E / ' df « = E / ^ ‘ df « + E / f «P ' df « 


a 1 


a 1 


a,|3 1 




a 1 


a,i 1 ^^a,i 


( 8 . 20 ) 

( 8 . 21 ) 


( 8 . 22 ) 


- -2 f du « -2 ( u .« - u ») • 


a 1 


(8.23) 


(2^ is more complicated. On one hand 

2 - 2 - 2 - 

^ = E j f »P ' df a = “E 1 f Pa ’ a = “E J f a P ‘ d ?p ' 


a / P 1 


a,|i ]_ 


a / P 1 


® = ^ E / ■ ( df « ■ df p) • 


a,|3 1 


On the other hand the chain rule gives 


dU ap = 


J ^ dr„ • + 


0U 


a|3 


d r a, i 


a,i 


dr. 


dr, 


P,i 


P,i 


(v„u ap ) 


•dr 


(V«p) 


■dr Q 


dU ap = (V a U ap ) • dr a + (Vp U ap ) * df p , 


-f 


a|3 


^ a ^ap 


=> dU„„ = 


a|3 


© 

© 


- f «p • (dr a - dr p ) . 

\ E f £ «P ' ( df « “ df p) = “ \ E / d0 «P ' 

Z cx,[3 1 

“ ~ E ( U 2ap " U lap) • 


a,|3 1 


a,P 


(8.24) 


(8.25) 


(8.26) 


Download free eBooks at bookboon.com 


n 













MODERN INTRODUCTORY MECHANICS PART II 


PARTICLE INTERACTIONS AND CENTRAL FORCES 


Therefore 


W = 
VV 12 


-2". -5 2”. 


a|3 

a, (3 

(a*P) / 


Define the total potential energy, 


U - I U a + z u raB = I u a + - z u a p, 


1 a(3 

a a< (3 


a,p 

(a*P) 


then 


W 12 = -u I = u, - u 2 . 

1 


Combining this with W = T — T 


12 2 1 ’ 

=> Tj. + Ui = T 2 + U 2 . 


(8.27) 


(8.28) 


(8.29) 


(8.30) 
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Total energy, kinetic + potential, is conserved. The n-particle Lagrange function in CM coordinates 
can then be written: 


1 -2 1 
L = T-U = —MIT 

2 2 


+ ^2 m a r a 2_ 2 U «" E U ap(| r a _r p|)- 


(8.31) 


a a a<(3 

8.2 TWO-BODY RELATIVE COORDINATES 

Let us now specialize to the 2-body problem with no external forces: 


l - 2 1 -2 1 -. 2 

L = — MR + — m.r, + — m,r, 
2 2 2 


hh - *d) 

t 

drop the bar, subscripts 


(8.32) 


Of course r x and r 2 are not independent, but satisfy 

Z m « = 0 => = 0 . ( 8 .33) 

a 

(There are always 3 + 3 independent coordinates for a 2-body system) Let us introduce 


r = r*i - r 2 = r, - r 2 


(no subscripts^ 
or primes on r) 


Picture: 



fixed origin 


Then 


(8.34) 


r = r. 


( ™ \ 

1 + ^ 

V m 2 ) 


r, = 


m 0 


r , 


m 1 + m 2 


(8.33) 
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or 


r = -r. 


r 

1 + 

V 


—\ -r — 

f 

®i 1 

mj 

r 2 — 


+ m,J 


Therefore 


1 -2 1 -2 1 
— m.r, + — m 9 r, = — 
2 2 2 


m, 


m. 


\ 2 / 

+ m 9 


\ m i + m 2 / 


m, 


so 


-2 1 -2 

m 1 r 1 + — m 2 r 2 


jxr 2 , 


where 


\ m i + m 2 / 


f 2 


(8.36) 


(8.37) 


(8.38) 




m 1 + m 2 


1-2 1-2 

=> L = — MR + - Lir - U(r). 
2 2 


(“reduced mass”) 
(8.39) 


Now find the Hamiltonian for the above L. Adopt R., r as generalized coordinates. 

dL 


OR; 


dL 

dr. 


= MR ± = 


(p = mr) 


= ^ r i = Pi, 


(p = fir) 

small p; no indices 


(8.40) 

(8.41) 


Z dL • x - ' dL 

R, +> — r - L , 

i dR, 1 r dr L 1 

P 2 p 2 

H = — + — + U(r) . 

2M 2p 

H cm 

H is a constant of the motion since it has no t dependence. 


(8.42) 


(8.43) 


Hcm = 

t 

const. 


P 2 

T ! 

H re i = 

*1 + 

U(r) 

2M 

2 p 



t 



of motion 

also 

const. 

, of 
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The total angular momentum is 

L = R x P + f x x Pl + r 2 x p 2 

i J r'iwr /) 


(relative angular momentum) 


Can rewrite l : 


Pl = = m. 


p 2 = m 2 r 2 = m 2 


( m \ 

rci 2 

vm 1 + m 2 ) 


( ™ \ 

, m i 

\ m i + m 2 / 


r = |ur 


r = -|ir 


= r x x Pl + r 2 x p 2 = 


m 0 


r x 


m 1 + m 2 


+ 




r x 


m 1 + m 2 




(8.44) 


(8.45) 

(8.46) 


(8.47) 
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~r itl + m 9 _ / _l\ _ _ 

6 = — -- r x hir = r x p , 

m 1 + m 2 v ' 

^>L = RxP + rxp. 


(8.48) 

(8.49) 


Already know L = 0. But L CM = R x P = 0 also => t = 0 . This has 2 important 
consequences: 


1. Since C is time independent and because C ■ p = t • r = 0 (_L to both P and r), 
it is clear that the motion is in a plane _L to C . 

2. I = | r x p | = p,r 2 | <j) | 


~rd(|) 

1 2 , fl . 4 

dA = — r d<p , base x height of tnanglej 

dA i 2 : t 

=> — = — r (j) = — = const, 
dt 2 2p, 

Radius vector sweeps out equal areas in equal times. Called Kepler’s Second Law. This is also true 
in the CM frame, but the rate is different. 



8.3 RUNGE-LENZ TREATMENT OF COULOMB FORCE 

Let’s examine the vectors which may be formed from C , p and r. The scalar products 

r • C = r-(rxp) = 0 , ( 8- 50 ) 

p • I = p • (r x p) = 0. ( 8 . 5 i) 

vanish as pointed out before. The vector products are 


M = t x r , 

(8.52) 

III 

X 

(8.53) 


They are simply related: 


/ 


[l M = [1 


\ 


l x r + l x r 

Vo / 


xp = i 


(8.54) 
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Next, consider 


A s i x p. 


Hamilton’s equations for r., p. are: 

3H. 


r, = 


Pi 


rel 


3 Pi 

5H rel 


(identity) 


dr L 


dU 

dr L ' 


But 


dU(r) dU dr 


dr. 


dr dr 


V 2 2 2 

r i + r 2 + r 3 f 


SO 


dr 

diT 


4 = 


2r 


r 


(r; + r 2 2 + r 3 f 


- 7 r du 

x p = -e x- 

r dr 


or 




M / du\ 

r IdrJ ’ 


Now 


M-*xr-(rxp)xr-rx(pxr), 


(A x (§ x C) = B (A • C) - C (A • §)) , 

=> M = p (r • r) - r (r • p) = pr 2 - r (r • p) , 
r • p = r (jir) = [irr , 


(8.55) 


(8.56) 

(8.57) 


(8.58) 


(8.59) 


(8.60) 
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In addition, 


e = 


_d_ (r 
dt 


r 

r r 


r . 


= ~y (rr - rr) 


(e r is unit vector along r). This means 
M = |Ltr 3 e r , 
and also that 


A = -^ir 2 e r 


\dr / 


All this has been done without reference to the form of U(r). Now assume 


(8.61) 


(8.62) 


(8.63) 


U(r) = - - , 
r 


k = G rriim 


r LL 2 


t 


Newton's gravitational 
constant 
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=^> F = -VU =-- e 


2 r 


Then 


A = -|ike r or — + |ike r j = 0 

A constant of the motion! This is the Runge-Lenz vector, A. 


(8.64) 


A = A + |ike r = (r x p) x p + |ik e r . 


Notice that 

dt 

A • e 


= 0 only for U =-. Can now show that 

r 

= [(r x p) x p] • l + |ike r • l = 0. 


(8.65) 


1 to rxp 


T 

rxp 


The vector A is in the plane of motion along with r and p . Can get the equation of motion 
as follows. 

A • r = A • r + p,kr , 


A ■ r = (l x p) • r = l ■ (p x r) = -C 2 , 
A • r = -t 2 + p,kr. 

Let (|) = angle between - A and r. Then (A = | A|) 

- Ar cos (|) = -i 2 + [Tkr , 


1 

r 


kp 


A \ 

1 + - COS (j) 

pk 


Define 


( 8 . 66 ) 


(8.67) 


a = — , 
pk 


g = -. eccentricity 

pk 


(8.69) 
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Gives equation of conic sections. (Assume t * 0) 


s > 1 hyperbola) 

runbound motion 
s=l parabola J 


0 < s < 1 ellipse) 

r bound 

s=0 circle J 


Most important case: ellipse 


focus 



Notice when (() = n, 0 (angle between - A & r), r = r max , r min : 
a 



1 + 8 


Thus | A | determines eccentricity, A (direction) determines direction of r max . 
“Major axis”, a: 


◄- 2a -► 
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01 1-8 + 1 + 8 v a 

a = r^7 2 = 


(8.72) 


Let’s write the equation of motion in more recognizeable terms. 



a. 1 A 

— = 1 + 8 cos (p , 

r '—-—' 

x 

r 


—(a + ex) = 1, 

r 


(a + ex ) 2 = r 2 = (x 2 + y 2 ) , 


a 2 + e 2 x 2 + 2aex - x 2 - y 2 = 0, 


0 2 as y 

=> x 2 - --- x + — 1L 


1-e 2 


a z 


1-e 2 1-e 2 


= 0 , 


=> (x-ea) 2 + 




a z 


1-e 2 1-e 2 


a 2 e 2 


:i-e 2 ) 


2 


= o. 


a z 


1 + 


l-e 2 v l-e 2 y (1-e 2 ) 2 


a z 


= a / 


(x - ea) 2 + - - - a 2 = 0, 
1-e 2 


(x-ea) y2 

a 2 a 2 (l-e 2 ) 
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Equation of an ellipse (if 0 < 8 < 1; just define x'= x - sa ) . Can read off b, “minor axis”: 



I 2 a 

b = aVl - s = , • (8.73) 

VI - 8 

Have confirmed that equation of motion is an ellipse when r is plotted (Keplers 1 st law). What 
it represents: 




r 



m 2 


It is an ellipse when viewed from a coordinate system centered on ny or m 2 because 



when viewed from the CM, motions still form ellipses. 


(8.74) 

(8.75) 
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^ r are measured from the CM, so the CM point is the common focus point for both ellipses. 
The eccentricity, 8, may be determined from the energy and the angular momentum as follows: 

A 2 = |^ + [ike r j = A 2 + (i 2 k 2 + 2[ik^ • e r , 

2 2 = [(r x p) x p] 2 = [C x p) • x p) = C • |p x {I x p)), 

?p 2 -p(?-p) 

.-. i 2 = e 2 p 2 . 


From before (A • r 



•d • e r = - -C 2 , 

r 

A 2 = <? 2 p 2 + jx 2 k 2 - ^ 


= 2\iC 2 


=E 



r j 


^ 2 k 2 , 


( 8 . 76 ) 
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(E is the total internal energy of the system) 

=> A 2 = 2p/ 2 E + [i\ 2 , 


8 = 


[ik 


( 2^ 2 E V2 

1 + -- 


\ 


pk 


Notice that E can be negative, zero, or positive. Again, for C * 0: 


(8.77) 

(8.78) 


8 > 1 

hyperbola 

E > 0 

8-1 

parabola 

E = 0 

0 < 8 < 1 

ellipse 

E < 0 

8 = 0 

circle 

E = 


uk 1 2 

i — 5 - =-pv 

It 2 2 

(v = total velocity viewed from ny or ny.) 

8.4 LAG RAN GIAN EQUATIONS OF MOTION 

Let s do it again from a Lagrangian point of view to see the more standard treatment. 


L = — MR 2 + — pr 2 - U(r) 
2 2 __ 

^rel 


Dependent variables: R(t) , r(t): 


dL d ( < 3 L^\ 

R : —-i — r = 0 

5R dt V5 r7 

r: 

dr dt \dr 


=> R = 


0 , 

n dU , 

[ir = - — e r , 
dr 


t = r x p = r x 




dU 


— e rl = °- 

dr , 


(8.79) 


(8.80) 

(8.81) 

(8.82) 
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Motion is in a plane, as before. Pick l along z: 

(cylindrical 
coordinates) 



r = r e r + r e r + xr4> + r(j> + r(j> e^. 

We have (Ch.2): 


e r - 4 >e d) t 


®<j> - • 

=> r = (r - ref) 2 ) e r + (2rtj> + rcf>) 

Substitute into (8.81) above: 

(- 12 \ dU 

e r : P ( r - j = - — , 

e ( k : 2r(f) + r<|) = 0 . 

These are the ecp-s of motion. Then, 


( 8 . 88 ) => 2 


4* (r 

T=> 2 In— 

4 > \px 


= - In 


\ 

v c 2y 


\ 2 

r ' 




\ c i/ 


-4 => r 2 d> = const. 

\4>J 


P, 


This is Kepler’s second law, derived again. Also 

_ 1 dU 
p dr 


(8.87) => r - 


p 2 r 3 
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( 8 . 88 ) 


(8.89) 
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Solving this nonlinear differential eq a gives r(t). Can derive a different for r(4>) as follows. 


dr _ dr dxj) _ 1 dr 
dt d(|) dt jar 2 d(|) 



|a d(j) 


Likewise 


jar^ 


(— r 2 ) 



t 



d(f» 



d 2 r 


dt 2 


: t 1 

<t> = - — 

ja r 

1 



[a dt d (|) 2 


t_ U jO 

^ \M« r 'y 



1 


„ * fi) 

d 2 r _ t 2 1 ^r/ 

dt 1 " " J7 r 2 d()) 2 


(8.91) 


(8.92) 
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Substitute into (8.90) above: 


e i 


d 2 i- 


e i 


2 2 -.1 2 2 5 

[i r d(J) p r 


1 dU 
p dr 


or 


(■ 


\rj + I = ji r 2 dU 

d(|) 2 r l 2 dr 


(8.93) 


Given U(r), can find r(<|)). Conversely, can find U(r) given some r((|)). 

Example 1: 

k . 1 

U(r) =-. Define u = — 

r r 


d 2 u 


d(() 


ku 

2 + u = -7T- 


(8.94) 


A simple linear nonhomogeneous difF. eq Q . Same as for Hooke s law oscillator with a constant force. 
Homogeneous solution: C cos((|) - (|)1) 


Particular solution: 

General solution: 


kp 


(like const, displacement of spring for const, force) 


u = 


1 

kp 


unknown 

(l + C cos ((|) - (j)J) 


(8.95) 


E q~ of conic sections again. Supply the initial conditions. Lets agree to measure: 



(([) measured 
from pericenter) 
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Value of C? 


—— = — (l + C cos (I)]), 

* • or 

nm *n ^ 


C cos 


(7 


- 1 . 


Taking derivative of u = — (l + C cos((|) - (jq)), we get 

a v ’ 


1 • c ■ a A \ d(|) 
- — r = - - sin((|) - <y — 


a 


(8.96) 


Because r = 0 at r = r min , and assuming C ^ 0, we get 


„ c ■ I d(J) 

0 =-sm ©•, — =^> (b-| = 0 

T 1 ,, T 1 f 


a 


C = 


dt 

a 


r min 


1 . 


From previous solution we know that r min = 


a 


1 + 8 


a 


=^> C = — (l + e) -1 = s. 


a 


(8.97) 


=^> — = — (l + 8 COS (|>). 

r a 

p 2 k 

[Note: Can also get r orr . directly from E = - — ~ 2 

o max min J r\ 


(8.98) 


r mi 

and then solving for r from the resulting quadratic eq~. See prob. 8.15.] 
Example 2: r = — e d ^ , u = ce - ^ (c > 0) 


P* = 


min T min r Tmin ? 


Exponential spiral: 


(d>0) 
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Force law? 


d 2 u 

d(|) 2 


cdV d * 


d 2 u 


f 


d u ur dU 

—7 + u = -, 

d(|) t dr 




pr 2 dU 
dr ' 


=> F(r) 


dU 

dr 




(attractive) 


(8.99) 


Inverse cube law. This is not the only type of motion in such a force field. This is a special type of 
solution of the inverse cube law corresponding to E = 0. 
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8.5 CELESTIAL MECHANICS 

Let’s derive some other useful results for the inverse square law, U(r) 
t 2 1 dU 


r = 


|i 2 r 3 dr 


k 

r 


( 8 . 100 ) 


f ^ t 2 i 1 dU 

J dt rr " — — + - — r = 0 # 

J ^ ^ r dr j 


1 l 

— ^rr 2 + -- + U(r) = E. (const.) 

2 2[ir 


This is just the statement of energy conservation since 


( 8 . 101 ) 


r + —— = r + r (|) = v , 


2 2 

[i r 


1 2 

E = - [iv + U(r) , 
2 


2 2 ( 

v = —Eh-1 . 

h v 


( 8 . 102 ) 


E is related to “a” by 


a = 


a 


a 


1 - 8 2 


1 - 


1 + 


ajtk 2 


2rE) 2rE 


^k 


k 

=> a = — — > 0 , 
2E 


t 

since E < 0 for an ellipse 


_ 2 k ( 2 1 } 

^ v = - ( - - “J , 


p v r a 


(8.103) 


for elliptical motion. 
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There is a connection between x, the period of elliptical motion, and a, the semi-major axis, that 

is well known for U(r ) =-. The connection is stated roughly as 

r 

2 3 

x oc a 

t 

“proportional to” 


and is known as Kepler's third law. We actually derived this rough form last semester from an 
invariance argument. Recall that 


■ J rel = \ VC" ~ U(r) . 


with U(r) = c r n under the transformation 

fr —» ar, 
t —» a 1 ^ t. 


(8.104) 




(8.105) 


behaves as 

L 


'rel a -^rel * 


^ E T of motion from 8 j L rel dt = 0 unchanged. For the inverse square force law 

U(r)-- , n-1 , 

r 



as we had before. To get the full form of Kepler s 3rd law consider 

T = J = 2j d*, (8.106) 

full V <|> 0 V <t>( < l > ) 

orbit 

( = r x p = r x ( = rp # = 

t 



T = 2 7 f d< t> r2 W' 
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/ S 2 \ 


r 2 = 


\[ikj 


(i 


+ 8 COS 


( t ) ) 


2 ' 


=^> X 


2 tJLf ^ 

_l M- 2 ^- 2 o (l+£ cos(j)) 2 



r 


SO 


X 


2 




_ 

COS(|>) 2 


(8.107) 


Useful integral: 


f 


d(j) 


s sin (|) 


(1+8 cos(|)) 2 (e 2 - l) (l + e cos <j>) _ £ 2 j 3/2 


tan 


-1 


(1 


\ <j> ^ 

e tan — 

’ 2 
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Plugging in the limits, <|) = 0, 71 , we get 


T 



(tan 1 (oo) - tan 1 (o)j, 
n ( sta y on \ 

2 single brach J 


X = 



\k/ 



f 


x 2 = 




Keplers feelings about discovering his 3— law: 


(8.108) 


“Now, since the dawn eight months ago, since the broad daylight three months ago, and since a 
few days ago, when the full sun illuminated my wonderful speculations, nothing holds me back. 
I yield freely to the sacred frenzy; I dare frankly to confess that I have stolen the golden vessels of 
the Egyptians to build a tabernacle for my God far from the bounds of Egypt. If you pardon me, 
I shall rejoice; if you reproach me, I shall endure. The die is cast, and I am writing the book - to 
be read either now or by posterity, it matters not. It can wait a century for a reader, as God himself 
has waited six thousand years for a witness.” 


From “The Discoverers” by D. Boorstin 
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Lets do the above integration from (f> = 0 to an arbitrary angle (|): 


f 3\y 2 ^ . 

pa J d(|) 

V k J o (l + 8 COS (I)') 2 


27Tt 


t 


- (i - 


dcj)’ 


o C 1 + 


8 COS 


<n 


i\2 


“mean anomaly” 


(8.109) 


Using the previous integral gives: 


2 Tit 
x 


2 tan 


-l 


^(1 - e) tan —^ 


JiT 


e-Jl - e 2 sin <|) 
(l + s cos <))) 


( 8 . 110 ) 


Careful use of this equation gives t((|)). But what we usually want is c|)(t)! The need to invert 
this relationship gives rise to Kepler’s equation. We can derive it as follows. Adopt the following 
coordinate system: 



... (x + ea) y 

ellipse: --—— + ±- = \ 

a b 

V]/ = “eccentric anomaly”. 


First, get relationship between (|) and \|/. 

x + as 


cos \|/ 


x = a(cos \|/ - s) , 


( 8 . 111 ) 


sin \|t 


i 


a 2 - (x + as) 2 _ ^ b 


£=> y 

b 


= b sin \|/ , 
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or 


= ai/l" 


s sm \|/ , 


so 


r 2 = x 2 + y 2 = a 2 (l - e cos \\i) , 


=> r = a(l - s cos \j/) 


Remember: 


11, 

— = — (1 + 8 COS (pj 

r a 

, a 1 a 

=>1 + 8 cos <p = — --- , — = 1 

a (1 - s cos t|/j a 

Taking the differential of both sides of this gives 

f , > 


- 8 


d(l + s cos (j>) = (1 — £ 2 ) d 


\1 - 8 COSVJ/ J 


d(|) = (l - e 2 ) 


2 x sm \|) 


d\)/ 


sin (|) (l - s cost)/) 


Moreover, by definition 


sin (|) = — 
r 


b sin \|/ 


a(l - s cost)/) A 


d(|) = Vl - 


d\|/ 


1-s cost)/ 


We now have 


d(|) 


(l + 8 COS <)>) 


VTT 


d\)/ 


1-s cost )/J 


(1 - 


(l - 8 cost)/) 2 


T 


use (8.114) and (8.116) 


7- ^7377 l 1 " 8 cost)/) dtp , 

t 1 - 8 ) 
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(8.114) 


(8.115) 


(8.116) 


(8.117) 
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so now 


2jit 


2 Tit 


. (l-e 2 p f - **- - T , 

1 ’ i (1 + 8 cosf) 2 

¥ 

= J (l - 8 cos\|/')d\|/' = v|/ - s sin v|/ . 


(8.118) 


Eq (8.118) is Kepler’s equation. One can also show (integrate (8.116)) that 

(b /I + s vi/ 

tan — = ,/- tan — . 

2 VI - £ 2 

So, the method of finding (|)(t) using Kepler’s eq~ is: 


(8.119) 


1. Solve Kepler’s eq- approximately for V|/(t). 

2. Use (8.99) to find (|)(t) from 


4>( t) = 2 tan 


-i 


.£ 


+ 8 V|/(t) 

- tan - 1 -— 1 
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8.6 GENERAL RELATIVITY MODIFICATION 

I now want to talk a little about a small modification to these orbits caused by Einsteins general 
theory of relativity. Start with the effective (approx.) potential for perihelion (Mercury) precession. I 
am leaving out terms in effective potential which do not cause precession. This potential is observer 
dependent. The effective potential with respect to an observer at r = oo: 


U^f(r) = - 


GMm 

r 


3 (GMmA 2 

me 2 ^ r ) • 


( 8 . 120 ) 


This is not the same as the effective potential as, for example in Marion, which is with respect to 
the observer in orbit. 


Let’s put this into our orbit equation: 


d(|) 


r / 1 

+ — 


[ir 2 dU e 


C 2 dr 


Let u = — 


d 2 u 

df 


Solve: 


+ u = 


( 

GMm + 

C 


me 


(GMm) 2 u 


d 2 u 

dip" 


6 = 

Solution: u = 




6G 2 M 2 mp 

<? 2 c 2 

kp 


(1 - 8) 

1 

a (l - 8) 

Most important change: 

<|) = 0 


|l + e' cos |VI - 8(|)jJ . 


analog of eccenticity 
constant 


4 > = 


271 

Vl - 6 


1_ 

r 

1 

r 


r min 


again. 


( 8 . 121 ) 


( 8 . 122 ) 

(8.123) 

(8.124) 
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Notice ^ 11 > 2 x => orbit precesses positively (8 > 0). 

■Jl — 8 



2n 




L , 8} 

2; 


2n^l + —j = 2n + A r 


. „ 6jtG 2 M 2 mu 

A = Jt8 = -——— 

l 2 c 2 


From a def": 

^ 2 = apk = a (l - e 2 ) pk = a (l - e 2 ) GMmp , 
6 jtGM 6jtGM 


A = 


ac 


ac 


(i - 

Amounts to "43 sec. of arc/century for Mercury. 


Newtonian effective potential: 


rel = - ^ 

2 2pr 


+ U(r) . 


= uNewton^ 
ef f 


Define, 


L centrifugal 


d 

dr 


/2 


2pr 


pr 


(8.125) 

(8.126) 


(8.127) 


(8.128) 


(8.129) 


Looks like: 


T ..Newton 
U eff 
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The full general relativity effective potential including centrifugal term: 



New G.R. effect: capture (realized in “black holes”) 
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8.7 ORBITAL STABILITY 

Equation which describes orbital motion (from (8.90)): 


r - 


[i 2 r 3 


1 dU 
\i dr 


Let r Q be the radius of a nearly circular orbit, and let r 

orbit (f = — ) 

L 


x 


= r 0 + x, where— « ] 


r o = " 


f(r 0 ) 


Then 


a -® ( 


r o + x )" = f ( r o + X ) • 


Assume 


f (r 0 + x) = f(r 0 ) + (r - r 0 ) £ (r 0 ) + ... . 


X 


Also 


( r 0 + x )" 3 = r o 


-3 


X 


1 + - = r c 

V r J 


-3 


1 - - 

xj 


Then 


x - 


-1 r r 


-3 


VEy 


r 0/ 


= f ( r o) + xf '( r o) 


t 


substitute the above 

(-) fVo) 


X + 




-T 


r 0 / 


= f(r 0 ) + Xf'(r 0 ) , 


Let 


a = - 


3f(r 0 ) 


- £’(r 0 ). 
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(8.131) 

(8.132) 

(8.133) 

(8.134) 

(8.135) 

(8.136) 

(8.137) 
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Then 


a 2 > 0 : x = c 1 sin at + c 2 cos at - stable 


a 2 = 0 : x = c ± t + c 2 ] 


2 ^ r\ |a|t , —10t|t j 

a < 0 : x = Cje 1 1 + c 2 e 1 1 J 

Condition for stability: 


unstable 


3f ( r o) 


+ f'(r 0 ) < 0 , 


or (notice f<0) 


F 'Co) 

F W 


> -3 . 


Try this out for F(r) =- 


U = - 


(a - 1) r-‘ 


, ak 

F ( r > = 3CT ' 


F ' ( r °) 

F(r„) 


stable if -a > -3 or a < 3. (Must also have a>lforas U^O as r^oo 
interpretable (k > 0 for all cases): 


a < 3 case: 


Newton 
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a = 3 case: 


Newton 



8.8 VIRIAL THEOREM 

Now consider (x = period of motion) 

(^ (£<r) r)) t =ijA(f ( r|r)dt, ( 8 . 14 1 ) 

= - [f(r)r] | = 0 . 
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for a periodic system. (Even if it is not strictly periodic, can let x —> and make or assume 


lim — [f(r) r 1 

T —^ 00 T L J 


= 0 . 


(8.142) 


if the system is bounded .) Now we have 


— (f(r)r) = f'(r )r 2 + f(r)r , 
dt 


(8.143) 


(f'(r)r 2 ) + (f(r)r) = 0 . (averaged over time) 


(8.144) 


Can now use 


r - 


(i 2 r 3 


1 dU 
[i dr 


so 


f ’(r) r 2 + 


2 . f(r) 


[i 2 r 3 


- ( f v) f- 

[i \ dr 


(8.145) 
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First example: f (r) = const. 

=> <0 = 0 , 



or, in words: 


(8.146) 


Average centrifugal “force” = Average radial force 


Must be true for any periodic system. Another example: f(r) = 


1 


— pr + 


2pr 


= - r 


du' 

dr/ 


du' 


(T) = - (r 
x ' 2 \ dr/ 

This is just the “virial theorem” applied to central-force motion. Let’s take U(r) = 


(8.147) 


— —. Then 


<T) 


= \§- 


<T) = - - <U) . 


But, T + U = E => (T) + (U) = E , 
<U) = 2E <0 

[<T) = -E >0 

(For circular orbits, can take off ( ).) 


(8.148) 


for bounded motion 


(8.149) 


Application: galaxy clusters => must be “dark matter” in the universe. Outweighs visible matter 
by a factor of 5! 
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8.9 


PROBLEMS 

1. A system of point particles interact via forces which follow the “strong form” of Newtons 
second law. That is, the force of P on a, f a(3 , points along the instantaneous line connecting 
them, as shown. 



Given the usual connection between fixed (r a ) and center of mass (r a ) coordinates, 
4 = r a + R, 


(r = — \ f 5 M is total mass) and the total force on a, 
M ^ a 


p = F (e) + 

rr u. a ^ 


2 f «p ' 


t t 

external 1 /interna! 
force on a J v force 


Show that the total torque, L = ^ L a ( L a = r a x p a ), for an external force of 
the form, a 


= ni a g , 


is simply given by 

L = R x F (e) , 


where F (e) = ^ F^ e) is the total external force. 
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2. Show that 


3. 


d 2 J 

dt 2 


- - £ , (^ = C x p ) 

[i 


when the potential U(r) 


1 

2 


kr 2 . Comment on the expected solution. 


Investigate pure radial bound motion ( £ = 0 , E<0 ) in an inverse square force field, 

F (r) = - ^2 5 f° r two point masses. Let “h” be the maximum amplitude for pure 
radial motion. Imagine the particles move through one another as they pass. How is the 
period of pure radial motion (x rad ) related to the period of a circular orbit of readius r (x b ) 
when h = r? 


4. 


As a follow-on to prob. 3, show that the equation of bound radial motion is a cycloid in 
time; that is, 


r 


h 

2 


(1 - cos 


4 » 


t 


h 

2 



(<|> - sin <|)) 


r 


where h is the maximum radial distance between the masses. What is the relationship 
between r and t when E = 0? (E = total energy.) [Hint: Try making a change of variable, 
r = 2 (1-cos <|>)> directly in the r integral.] 


> Apply now 
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5. You are given that the distance of separation of two planets is given by 
r = d(l + sin (|>), 

where d > 0 is a constant and (|) is the angle. (Ignore the fact that the planets strike one 
another at 0 = 270°.) 

a) What central force law is responsible for this motion? 

[Ans.: F (r) = -3 

pr 

b) Find the total kinetic plus potential energy of this system in the center of mass frame, 
given that U(oo) = 0. 

[Ans.: E = 0.] 


6. Show that the maximum radial velocity of two planets in an elliptical orbit about one 
another is 


2nea 


r max T (l_ g 2)l/2 


where 8 is the eccentricity, T is the period and “a” is the semimajor axis. 

7. Show that the product of the magnitudes of the maximum and minimum relative velocities 
of two bodies in elliptic motion is (27ia/x)2. 

8. Starting from (8.110) of the notes, derive to first order in 8 only (by a Taylor series expansion) 
the approximate result, 


, /L v 2jtt _ . 2jrt 

(f)(t) « - + 2e sm(-) . 

T T 


9. (a) Do a simple-minded solution to Keplers equation, Eq.(8.118) of the notes. Assume 8 

is small and do an expansion of the form 

V = Vo + 8Vl + s >3 + ••• 

« „ 2 Tit 

Find V|/ 0 , V |/ 1 and \|/ 2 in terms of the “mean anomoly”, M = . 

(b) Show that to order 8 you get the same result for (|)(t) as in prob. 8.8 above. 
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10. Show, for the inverse square law of gravitation (“a” is the major axis, “s” eccentricity), 


11 . 


12 . 


13. 


< 72 > 


(b) <^> 

(c) <A> 


2 y/l^2 ' 


1+ 2 


a 2 (1 —s 2 ) 


1 


2 ^Jl- 


; < \ r n-2^ > <|>/ 


where 

< . . . > 

< . . . > 

(a) Given 


=“/••• dt, 


2tc 


2k 


/• • • d<|). 


In - 


2n cltj) 

J ( 1+8COS (j)) n ' 
0 


1 S 

A)t 

2 <r n_1 > t 3 

l 1 + n+1 

dsj 

v (1-s 2 ) 3/2 J 


can you find a relation between In and In+1 ? [Hint: Differentiate in 8.] 

(b) Use this to relate <r n > t to <r n-1 > t . [See prob. 8.10 for notation. 

/ 

Ans.: <r n > t = a(l-s 2 ) 3/2 

Verify directly that 

<U(r)> t = 2E , 

k p 2 k 

where U ( r) = - “ and E is the total (internal) energy (=- r ). You may use any 

result from prob. 8.10 above to do this. 


Show that 


1 

a 


< r n >t = fy <r n+1 > 




where the < • • • > notation indicates an average over an orbital period with respect to the 
subscripted quantity. V|/ is the angular quantity used in the Kepler equation derivation, 
given by Eq.(8.113) of the notes: 

r = a (1 - s cos \|/), 

where “a” is the semimajor axis and “s” is eccentricity. 
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14. Show that 



( l =| r x p | ) for any periodic central force motion. You may need a result from 
the text. 


15.(a) Using the idea of effective gravitational potential, show that the turning points in an 
elliptical orbit under Newtons force law are given by 


-k ±. k* + ^ 


r = 




2E 


where E < 0 is the total center of mass energy of the system. 

(b) Show that the r min , r max you find agree with Eqs.(8.70), (8.71). 
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16. A paricle in a nearly circular orbit about a planet suddenly has its tangential velocity 
changed by a small amount, Sv. Show that 

^ 6^ _ 5a 6b 8v 8 e 

C a b v E 

where a is the semimajor axis, b is the semiminor axis, v is the magnitude of velocity, E 
is the (total internal) energy and l is the angular momentum. 

Extra credit: The particle s velocity is now changed by an amount 5v when it is at r = 
rmax and 8 is no longer small. Show that the changes in the orbital parameters are now: 

5^ _ 6v 6a 6 e 6b 2 6v 8 e _ 1 - e 6v 

l v ' a E f b 1 + e v f E l + ev 

17. According to the stability analysis, stable circular orbits should exist for a force law 
of the form, 

k 

F(r) = - ^ , (a<3 , k>0). 

a) Find how the period, x, is related to the radius of the circular orbit, a (i.e., the analog 
of Kepler s third law. The text presents two ways of doing this; extra credit if you get 
the proportionality constant.) 

b) Apply the virial theorem to this force law and find how the time averaged kinetic and 
potentials energies, < T > and < u h are related to the total energy, E, for circular orbits. 

18. A physicist is standing on the Earths equator when a very strange thing happens (the sort 
of thing which occurs only on physics homework sets). Suddenly, the Earths radius shrinks 
to zero, although its mass remains unchanged. This leaves the physicist in orbit about the 
point Earth. Calculate 8 (eccentricity) and r min of the physicists orbit. (R £ = 6.37x103 km, 
M e - 5.98xl0 27 gm.) 

Extra credit: Get formulas for 8 and rmin when the physicist is located at a latitude X. 
What happens when he or she is located at one of the poles? 

Other Problems 

k 

19. Let us investigate the attractive inverse cube force law, F ( r ) = - r 3 (k> 0 ), a bit more. 

f 2 

a) Find the general solutions, r(t) and r((|)), when k = —. {£ is the magnitude of the 

[i 

relative angular momentum and p is the reduced mass.) 

b) Describe the qualitative motion of a planet for —> k> 0 . In particular, is there any 
bound motion? 
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20 . 


21 . 


Two point planets (reduced mass p) are in a circular orbit about one another at a 

k 

distance r Q with a force law F(r) = - (l<a<3). The orbit is slightly perturbed 

by a passing comet. Find the period, T, of small oscillations induced by the comet about 
the circular orbit. Show that in general, this period is not equal to the orbital period, 


^orbit ~~ 2ft 


1/2 


( r o) 


(a + l)/2 


Consider two planets which experience a repulsive Coulomb potential, U ( r ) = ^ (k>0 ). 
The total internal energy is E (>0) and the magnitude of the relative angular momentum 
is £. 


a) Show that the distance of closest approach, r min , is given by 


k 

2E 


1 + Jl + 


2E l 2 
pk 2 


b) Show that the radius as a function of angle, (f>, measured from the point of closest 
approach, is 


1 

r 


kp 

a2 


-i + Jl + ^4- cos ® 

\ pk 
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22. Prob. 8.11(b) resulted in the connection, 


<r n > t = a(l-s 2 )3/2fi + A Y <r n ~ 1 >t 

t aU 8 ’ Y n+1 dsAd-£ 2 ) 3/ Y ' 

where 8 was orbital eccentricity. Using this result twice, show that 


<r 2 > t = a 2 


(l + 2 s2 ) 

( 1 - 8 2 ) 2 


23. The central force responsible for the given trajectory in prob. 8.5 above was 
C 2 d 

F (r) = -3 —-. I also showed that the total energy of the system, E, was given by E = 0 
pr 

(given that U ( cc ) = 0 ). 


a) Draw a graph of the effective potential, U e f f (r) , vs . r , for this force law. 

Given E=0 for this trajectory, discuss the types of motion possible. 

b) Based upon the U e f f ( r) graph, and keeping £ fixed, find the minimum amount of 
additional energy, E min , to break free from the force center and move to r = oo. 

24. Investigate pure radial motion (£ =0) for the F(r) =-- central force law (k > 0) again, 

r 2 

but this time for E>0(U(°°) = 0). Show that a simple parametric representation, 
similar to prob. 8.4 above (which was for E < 0) is possible. Set the initial conditions to be 
r=0 at t=0. Solve as completely as possible. Show that your solution approaches the correct 

velocity, , as t —> 00 . 

p 

25. Given an attractive inverse cubic force law, F(r) = - — (k > 0 ), the plot of the effective 

r 

potential looks like the graph in the text: 


k i 

2 2[i 

(spirals in 
or out) 



Find the orbit solutions, r((|)), for cases (T) (total internal energy, E > 0) and (2) (E < 0). 


26. 


Given an attractive force law (k, a > 0) 


F(r) 



r 
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show that circular orbits for r sufficiently small are stable, but for r large they are not. 
Find the equation which determines the value of r where the transition from stable to 
unstable orbits occurs. 

27. In the text I derived 

dA _ C 
dt 2[i 

where l— \C\ {C= r X p ) and |U is the reduced mass. “A” is the area swept out by 

the radius vector r from the point of view of an observer on ny or ny. Now show that 

dA' 2 JA dA 

dt ~ ^ + m2J dt ' 

the area swept out by r in the center of mass frame of reference. Note this 

< 1 . 


where A’ is 
means that 

A' 

A 
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9 SCATTERING AND COLLISIONS 
OF PARTICLES 

9.1 COULOMB SCATTERING 

I will discuss scattering next since it builds on the results of the last chapter. 


Consider E > 0 (s > 1) motion in an attractive Coulomb field. Picture: 



Follow the course of r during scattering from m 2 ’s point of view: 



Equation of orbit (an hyperbola, £ > 1): 


11 / \ 
— = — (1 + £ COS (p) , 

r a 


1 

r min 



+ 8 ), 


(9.1) 

(9.2) 
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8 


a 


( 

1 + 

\ 


2E l 2 ' 
^k 2 y 


1/2 

r 



Notice: 25 + 0 = n. Alsc 


5 +0= 71, 

Jt + 0 

=> © = - 


(9.3) 

(9.4) 


(9.3) 


One gets 0 = 


it 


focus of the hyperbola.) 


for the repulsive force case, (where the scattering center is now the exterior 


Scattering event looks completely different from CM frame: 



direction 

Paths are also hyperbolas here (r ', r ' are just rescaling of r.) No matter which picture you prefer, 
it is clear that 0, 0 are quantities relating to the direction of r and therefore can be thought of 
as being measured in any frame. However, we will actually define scattering angles with respect to 
velocity vectors, which will therefore take on different values in alternate inertial frames. 

0 = positive angle between initial and final velocity vectors for either ny or m 2 in CM frame. 


Download free eBooks at bookboon.com 


55 








MODERN INTRODUCTORY MECHANICS PART II 


SCATTERING AND COLLISIONS OF PARTICLES 


Things are simpler in the CM; however, this is not the usual experimental situation. Eventually, 
we will learn how to translate our results in the CM to other frames of reference. Notice that only 
a certain range of angles for (|) are now permitted for s > 1: 

1 1 , v 

— = — (1 + e cos q)j, 

r a 

11, ^ 

=> — = — (1 + 8 COS©), 

oo a 

=> cos© = - — (allows 

8 

of 0 ; 

=> only angles cos c() >-are allowed. (Values greater than this would say that r is negative.) 

Go back to the Runge-Lenz vector to see it from another viewpoint: 

A = A xp + pke r , 

=> A • p = pke r - p . 


2 symmetric values (9.6) 

take the + value) 
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Remember, A is directed along the symmetry Saxis, in the direction opposite to r . In particular, 
if we consider the initial situation with r — 1 * °° , 

A • = pke roo - p OT . 

A • p m = Ap OT cos © , e roo - p^ = -p„ , 

=> Ap x cos 0 = -pkp m , 

=^> cos © = - — = - - , as before. 

A 8 

9.2 DIFFERENTIAL CROSS SECTIONS 

Need some more concepts for scattering. Let 

I = flux in the incident beam (# particles per unit area per unit time) 

dN = number of particles through a ring of radius b and width db in the incident beam per unit time. 



«-large distance-> 


dN = (2jtb|db|)l . 


(9.7) 


The particles passing through the ring are scattered through the angles between 0 and 0 + d0. dN 
becomes a scattering concept when we assume that b = b(0): 


dN (0) = (2jtb (0)|db (0)|) I . 


(9.8) 


dN (0) 


number of particles scattered into (0, 0 + d0) per unit time per incident flux. 
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This is now a quantity which is independent of I. We now define the differential cross section as 
(take d0 positive) 


da (0) dN 1 

d0 d0 I 


= 2k b(0) 


db (0) 
d0 


Intrinsically positive. We usually use solid angle. 


(9.9) 


d£2 s sin 0 d0d(|). 


(9.10) 


We will assume azimuthal symmetry here, so 
d£2 = 2k sin 0d0 , 


^( 9 ) - -A_ 

d£2 v sin 0 


db (0) 


d0 


(9.11) 


9.3 RUTHERFORD SCATTERING IN THE CENTER OF MASS FRAME 

For the Coulomb problem, remember 


A = p 2 k 2 + 2p/ 2 E, (A = pks) 
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At great distances, 


* 


0 


E = kinetic + potential , 

Poo 

=> E = - (two body problem form), 

2p 

f = |r x p| = b p m . 


sin 0 = 


_ iy2pE _ fp„ _ bp" 


But 


tan © = - 


bp" 

A 

pk 

A 


bp" 

pk 


© = 


0 + jt 


Q 

tan 0 = - cot — , 
2 


cot 


0 bp" 


pk 


( 9 . 12 ) 
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Provides the necessary connection between b and 0 for the Coulomb problem. Now get 


db uk d 0 

— = ^ — cot - , 

d0 p 2 d0 2 


jik 1 

2p 2 ■ 2 0 

sm — 


Therefore 






\P 1 ) sin 0 


sin — 
2 


. 0 

e cot 2' 


. „ , . 0 e 

sin 0 = 2 sin — cos — , 

2 2 



(9.13) 


True also for repulsive Coulomb force. Indeed, unchanged by quantum mechanics (non-relativistic), 
except that k —> ±(zze 2 ) for electrically charged particles. (Interpretation of it in quantum 
mechanics is completely different, however.) 


First measured experimentally by H. Geiger and E. Marsden. By careful observation in a darkened 
room, they found (by counting!) that approximately one in every eight thousand a-particles (Helium 
nuclei) was backscattered from a thin gold target. Lets see if we can understand the —-— factor 
backscattering fraction from the above formula. 8000 
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Model: (each atomic layer) 


atom 


§ 

• 

unit cell 


1 a | 

| (can 


assume other shapes) 


Assume layers are randomly oriented with respect to each other (i.e., one atom in front does not 
“shadow” an atom behind.) 


Each a-particle can be imagined to pass through the unit cell shown as may times as there are 

layers of atoms. We will take “backscattering” to mean — < 0 c n 

2 


do _ dN 1 

dQ dQ I 


1 (tik^ 

4 Ip i, 


sin 


0 ' 


dN = I da = 


pLj 


dQ 


.40' 
sin — 

2 


^back ^^back ^ ^ 


§ 1 / 


dQ 


VPoo/ t sin 4 — 


0 

For this rough estimate, we will set sin 4 - -> 1 over this range of 0, and since the solid angle 

jt ^ 

corresponding to — < 0 < ji is 271 is 2 tt, we then have approximately, 


Jt 

^back ” “ ^ 


f \ik) 

\ ) 


For scattering as off gold foil: 

k Zze 2 , |i - m a , 

Z = 79 , z = 2, 

/ 2 \ 

Zze m. i 


N. 


- *1 


m aVa / 
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Let us assume that (nonrelativistic formula) 

(I looked it up) 


1 2 

- m «v« - 5MeV 


( 1 MeV - 1.6 x 10 6 erg ). We also need 

I e | = 4.803 10" 10 esu 

m a = 6.68 10 24 gm 


P. = m a v a = 1.034 x 10“ 14 gm 


cm 


sec 


Still need I. In our case (imagine a single particle passing through the sample; also imagine 
multiplying both sides of the above equation for N back by the total time of the experiment so that 
N back is a pure number): 

I = number scattering processes per alpha particle / unit cell 


I = 


n , 
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where n is the number of atomic layers in the sample thickness. Sample thickness was ^ 1 \i = 10' 4 cm 
(also looked this up). So 

10“ 4 cm 

n = -, 

a 


I = 


10 cm 


Some other calculations giving a: 


P gold =19.28 


gm 

___ 3 

cm 


N a = 


Avagadro 

1 

6.023 x 10 23 , m 

gm 


gold 


-197 amu 


atoms 

cm 3 / gold 


( # atoms | , „ „„ 

p -“ x Kt" r 


^6.0 2 3 x 10 23X 
197 


= 5.9 x 10 


22 


atoms 


cm 


This means 

a - (5.9 x lO 22 )"^ 3 = 2.57 x 10' 8 cm 
=^> I - 5.90 x 10 18 cm" 2 . 


Putting all the pieces together now gives us 


=> N 


back 


4.76 x 10‘ 5 


1 

21,000 * 


too small " 2. One reason: did not integrate the cross section. This means we have underestimated 
the number of interactions. (But see the homework problem!) 
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9.4 SIMPLE TREATMENT OF LIGHT DEFLECTION 

Since we talked a little about one general relativity effect last chapter (perihelion procession), will indicate 
another result here. For very small angular defections, the above relationship between b and 0 becomes 


b 


„ jik 1 



=> 0 ~ 



1 

b 




m i m 2 

m 1 + m 2 


, k 


\ 

Git^m., 


Use this to model a light ray grazing the radius of the Sun: 



Therefore, take b = R s , m 2 = M s (Suns radius, mass.) What to do about p ~ m 1 , p x ? For light, 
hv 

P = —• 
c 


Use the relativistic connection between mass and energy to replace (E = pc) 


m, 


hv 


=> 0 


2 liU GM S OU 


hv 

c 


R= 


2GM s 

c 2 R c 


Almost correct, but too small by a factor of 2! (This was, in fact, Einstein’s original result, which 
he modified later.) Einstein’s general relativity gives the correct result: 


0 - 


4GM S 

c 2 R s 


1.75" (sec) of arc . 
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Another new thing about light scattering in general relativity: glories. 



A tiny amount of 
of light is actually 
backscattered. 
Appears as a sort of 
"halo" 


9.5 CROSS SECTION COOKBOOK 

We have been discussing a special case of scattering, the inverse square force law. We need some 
cookbook formulas for doing other force laws. Will give differential cross section in the CM frame. 
Need general connection between b and 0. Start: 



l 2 

2pr 2 


+ U(r) . 
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Can “reduce the problem to quadratures.” Solve for r: 


r 



(E - U(r)) - 


2„2 


m r 


(9.14) 


t 

have to pick correct root 
(U e ff(r) picture good for this) 


d(j) dt 

d(b = — 1 — dr = 
dt dr 


<l> 


dr , but (j) = 


pr 


d(j) 


1 / 

/r 2 dr 

2u (E - U(r)) - P — 


(9.15) 


Assuming E > 0 , E = f-h and integrating on dr from r . to r = oo , we get 
\ 2 IV 

r = oo, we get (£ = b p, x ) : 

0 defined positive 


1 

oo 

0 = +b f 

Pnin 



(9.16) 


Warning: r min itself is a function of b, in general. 


Cookbook steps: 


1. Evaluate r .(b). 


2. Do integral. 

3. Use 0 = 71 + ® or 0 = —-— in the attractive or repulsive cases, respectively, to 


find b(0). 


4. Plug in 


da 


dQ sin 0 


db 


d0 


Another warning: does not work for all potentials, U(r). Most require lim U(r) = 0 . 
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Just to get a feeling for using the cookbook method, do the repulsive Coulomb case (different from 
attractive case). Picture: 



rrij - m 2 


1 

20 + 0 = it 


Jump to result: (you will verify this in a problem) 

kp 


cos 0 = 


Can write as 


P»b 


1 + 


kp N 

\P»b 7 


cos 2 0 = 


kp 

P»b 


1 + 


kp ^ 

P»b 


2 ' 


(9.17) 


=> tan 2 © = 

M 

2 

or b 2 = 

fkp) 


kp j 


(P 2 J 


tan 2 © 


Choose b = + 


^kp N 

v P ») 


tan © , © = —-- , 


t 


different from attractive case 


, kp 0 
b = —t- cot - , 


same as before. 
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9.6 CONNECTION BETWEEN LAB AND CM FRAMES 

Only problem: cross sections usually not measured in the CM frame. 

Lab Frame CM Frame 

(usual expt. situation) (P = 0) 


initial: 


u„ 



m 2 


u 


1 




u 


2 
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Only showing initial and final velocities; not discussing dynamics but kniematics. We will assume 
elastic collisions - no heat generated or mass/energy lost. We will conserve both momentum and 
energy. We will also assume that all motion takes place in one plane (azimuthal symmetry as we 
assumed before in the cross section discussion). We will take the angles \|/, etc. as positive; if 
not, we can always re-orient our axis so that they are. 

Only difference between the two viewpoints: observed by two people who have a relative velocity, 
V. All primed and umprimed quantities are related by V: 


initial 


final 


u 1 = + V 

(9.18a) 

v 1 = + V 

(9.19a) 

u 2 = u 2 + V = 0 

(9.18b) 

v 2 = v 2 + V 

(9.19b) 

=> u 2 = -V 





notation 

u - initial velocity 
v - final velocity 
1,2 - which particle 

prime, unprimed - CM, Lab frame, respectively 


By definition, 


lab coordinates 


■£2»A' 


R = -—-(m^ + m 2 r 2 ), 


+ m 2 


V = ---(m.u, + m 2 u 2 ), (u 2 = 0) 


m 1 + m 2 


ntu, 

V = - 


m 1 + m 2 


(- <) 


( 9 . 20 ) 
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In the CM coordinate system (P (after = 0 . Therefore 

J ' 'before 

nqCq + m 2 u 2 = 0 => nquq = m 2 u 2 , 


m 1 v 1 + m 2 v 2 = 0 => m 1 v 1 = m 2 v 2 


( 9 . 21 ) 

( 9 . 22 ) 


(u 15 v 1 etc. are magnitudes only.) Assume any potential that exists between the particles 0 0 as 
distances 0 oo; then, far enough apart, the energy is purely kinetic. Let 

T 0 = total energy in CM frame 

Mbefore = Mafter ' 


1 , 2 1 < 2 1 2 1 2 

— HULL + - m 9 u 9 = - IILV, + - m 9 v 9 , 
2 2 2 2 


irqu/ + 


\ 

2 



’2 ’2 

V m 2 ) 

m 2 u 1 = m 1 v 1 + m 2 

\ m 2 / 


Vi 2 , 


(9.23) 


Also 




2 

’2 '2 




u 2 + m 2 u 2 = 



\ m 2 ) 


\ m 2 ) 


m 1 v 2 2 + m 2 v 2 2 , 


u 2 = v 2 . 


(9.24) 


Thus, to summarize: 


. . m 9 m 9 

u i = Vi = —- u 2 = —- v 2 


m 1 


m, 


=> only 1 unknown velocity magnitude in the CM frame. Other unknown: 0 ((|) = n - 0). Lets 
say we measure these 2 things in a given collision. How are they related to quantities in the Lab 
frame? From before, 


V = 


m 1 + m 2 


= -u 2 , 
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but 


m, ■ 

u 2 --U lf 

m 0 


m. + m 9 m. 

U 1 = —- 1 • — u, = 

m i m 2 


m \ 

1 + ^ 

m 2 / 


u, 


t 


known or measured, say 
Not specified yet: V , V 2 - We have, 

Vi = vj + V, 


v 2 = v 2 + V. 


Have to start invoking angles now: 


v 2 = Vj 2 + V 2 + 2V • vj, 
v 2 = v 2 2 + V 2 + 2V • vj. 


V = -u. 


V = u 2 = 


m, 


m 


u i 

2 / 


refer to 
the figures 


m, 


2V • v x = 2Vv 1 cos 0 = 2 — u 1 u 1 cos 0 


m. 


2V • v 2 = 2 Vv 2 cos(jt - 0) = -2Vv 2 cos 0 


m i . m, 

= -2 — u x — c 
m 2 m 2 


2 ' 2 

Vi = Ui + 


fm V 


m 2 / 


m 


u x + 2 — u x cos 0 , 
m. 




\ 

2 

\ 


2 ' 2 
=> v, = u, 

1 + 


+ 2 


COS 0 



V m 2 ) 


\ m 2 ) 
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V, = 


V 2 = 2 




\ 

2 

/ \ 

2 


\ 

= 


U 1 

+ 

m i ’ 

— Ui 

- 2 

™ 1 

U 1 


v m 2 

) 


\ m 2 ) 


V m 2 

/ 


cos 0, 


/ 'm 1 . ^ 

— 

\ m 2 / 


v 2 = 4 Ui 2 


/ m \ 
'm,' 

m 2 / 


(l - COS 0) , 


. 2 0 
sin — . 

2 


Therefore v l5 v 2 are known if u x , 0 are known. Only things left: ip, 'E, 

® v, = v; + v, 


v 2 = v 2 + V . 


© 

x : v 2 cos xp = v x cos 0 + V , 
y : v : sin xp = v : sin 0 . 


(9.27) 


(9.28) 

(9.29) 
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© 


X 


-COS0 

1 

: v 2 cos ^ = v 2 cos (jt - 0) + V , 


-v 2 sin ^ = —v 2 sin (jt - 0) 


(9.30) 

(9.31) 


T 

sin0 


Divide — from (l): 

x 


tan = 


sin 0 


v x sin 0 
v i cos 0 + V cos 0 + Z. 


V, 


m. 




But — = —■— = — => tan tp = 

v x u x m 2 


sin 0 


n m i 

cos 0 + — 


m 0 


(9.32) 


(9.33) 


Says that the connection between 0 and V|/ is not unique under some circumstances. Look at 

m i 

—L » 1* 
m 0 


tan ip ~ — sin 0 


ITT 


(9.34) 


Given \|t, 2 values for 0. Other extreme, « ]_; 

m 2 

tan ip - tan 0 , 


ip ~ 0 . 


(9.35) 


One value of \\i => one value of 0. Geometry helps: 




V, 


rru 


V 

v' 


^ > i 

rru 
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However, \|/ andy^ determine 0 uniquely. Also, one value of 0 => always one value of \|/. Can now 
see what the two extreme cases above correspond to geometrically. 

-» r m 1 

More geometry: —- > 1 case: 
m 2 



There is a max value of \|/ in this case: 


sin ^ max 


Zi 

V 



(9.36) 


In this case one value of \\l => one value of 9. Rutherford case: —- « 1, essentially scattering 

m 2 

individual a-particles off the entire sample since the atomic centers are fixed=> unique determination. 


Special case: m i = m 2 


tan t|t = 


sin 0 
cos 0 + 1 


0 . 0 0 
2 sin — cos — 
_ 2 2 

-i 2 0 -2 

1 + cos-sin 


0 ' 
2 
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„ . 0 0 
2 sin — cos — Q 

- - —o— - - tan -' 

2 cos 2 — ^ 

2 


^ = 


Now consider equation ( 2 ) above. Divide —: 

x 


v, sin 0 


v 2 sin § 
v 2 cos § —v 2 cos 0 + V 


tan Sj = 


sin 0 


Q V 

- cos 0 + — 

V, 


However, V = v 


2 ’ 


(9.37) 


tan ^ 


sin 0 

- cos 0 + 1 


2 sin cos 
2 


2 sin 2 


0 

2 


0 

2 
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, c 0 , (it 6 \ 

tan E = cot — = tan-, 

2 12 2 


& n - 0 i ’^- 9 ^ _ <|> 

^ ~ 2 ^ 2 ' 

0 

For the m = m case we had ip = — , so 

12 2 


S + i|> = | a (m, = m 2 ) 




Let’s find relationships between the various kinetic energies. Definitions first: 

/ , . \ 


T ' 
■*•0 


= total K.E. in 


u / 


T N 

■4 


lab 

CM 


= final K.E. of 1% in 


frame. 


lab 

CM 


frame. 


Similarly for T , T 2 . Of course 


T 0 T : + T 2 , T 0 T + T 2 


( 9 . 38 ) 


( 9 . 39 ) 


( 9 . 40 ) 
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Now 


T = 

J-n 


1 2 


To = \ (m ^; 2 + m 2 u; 2 ) 


( 9 . 41 ) 

( 9 . 42 ) 


= ^ m i v i 2 


T i = - m i v i 


T, = 


m 2 v 2 


T 0 = 


- m 9 v 9 
2 


Lets try to express them all in terms of T . Remember 


u i = u, 


/ ty» \ 

1 + ^ 


\ 


m 


2 / 


m 9 . 

u x = —- u 2 r 

m, 


=> T 0 = 


m, 


m 0 


'm ' 


2 \ 


\ m 2 / 


/ „ \ 2 / 

1 + ^ 


W 


m 


2 / 


1 + ^ 


m 


2 / 


U l 2 f 


=> T = - 


n^m 2 + m 2 m 2 


( m i + m 2 ) 


u, = 


i |XU 2 . 

T 


reduced mass 


Li 

m _ r m 

■ L o - L o * 

m, 


( 9 . 43 ) 


so T ' < T 0 always. Also 


■ 1 . 2 1 . 2 

T, = — m.v, = — m.u, 
2 2 


1 2 

- m 1 u 1 

2 


1 + Hd 


m 0 


^ = 


m 0 


\ m i + m 2 j 


T . 
■"■0 • 


( 9 . 44 ) 
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- m 9 


T 9 = — m 9 v 9 = — m 9 
2 2 


m, 


u, 


\ m 2 


m, 


\ m 2 , 


u, 


1 + ^ 
m. 


. m.m, 

T,- 1 2 . T„. 


( m i + m 2 ) 


2 0 


( 9 . 45 ) 


T v 

We have — = —. Go back to earlier v expression: 

T n u: 


v n = 



(^ \ 

nii 

2 

f m i 'l 

r 

i 

\ m 2 ) 

+ 2 

\ m 2 / 


V, = 


U, 



/ \\ 

ITL 

1 + 

l 


V m 2 )/ 


1 + 


/_ \ 2 

V m 2 / 


+ 2 


\ 

m ' 


\ m 2 / 


cos 0 , 
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M 

2 

/ , 2 \ 

m 2 + m 1 + 2m 1 m 2 cos 0 

UJ 


v K + m 2 ) 2 , 


What about in terms of \|/? Must express cos0 as a function of \|/. 

sin 0 


tan ip = 


cos 0 + x 


/ ^ \ 

m 1 

x = —- 


m 


2 / 


[“Aside 1”: 


sin 0 

x =-cos 0 , 

tan ip 

sin 0 cos ip - cos 0 sin xjj 


x = 


sin xjj 


x 


sin (0 - ip) -| 

= ---- , useful later. J 

sin ip 


tan ip = 


1 - cos 0 


(9.46) 


(9.47) 


cos 2 0 + x 2 + 2x cos 0 
tan 2 ip (cos 2 0 + x 2 + 2x cos 0) = 1 - cos 2 0 , 

cos 2 0 (tan 2 ip + l) + cos 0 (2x tan 2 xpj + (x 2 tan 2 xp - l) = 0 

V V ' 

t 


cos ip 

A quadratic equation in cos 0. Solution to Ax 2 + Bx + C = 0 is 
-B ± x/b 2 - 4ac 


x = 


2A 
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SO 


-2x tan 2 xjj ± 


, „ 4 (l - x 2 tan 2 ib) 

1 4x 2 tan 4 xp + —--- 


cos 0 = 


cos ip 


cos xp 


=>cos 0 = -x sin 2 xp ± cos 2 xp . x 2 tan 4 xp + 


(l - x 2 tan 2 xp) 


cos xp 


(9.48) 


Inside the square root: 


x tan xp 


tan xp--— 

cos xp, 

-T 


= -x 2 tan 2 xp , 


=> cos 0 = -x sin 2 xp ± cos 2 xp ■ -^-x 2 tan 2 xp , (9.49) 

\ cos xp 

t 

would have written |cos \|/ |, but 
because of the ± signs, this does 
not matter. 


From above 


so 


T 0 


i 1 2 Q 

1 + — + — COS 0 

X X 



T 0 


1 + ? 


- 2 sin 2 xp ± 2 cos xp./— - sin 2 xp 


+ i ' 2 

x/ 


(9.50) 


(9.51) 
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Notice the numerator is a perfect square: 


/ 


cos tp ± J— - sin ip 


x 


= cos tf» + — - sin tf) 


1 

±2 cos if). — - sin tf) , 


x 


1 -2 1 .2 
= 1 + — - 2 sm tf) ± 2 cos tf)J— - sin tf) 


x 


x 


1 + 


X 


1 .2 

cos tf) ± — - sin tf) 


T 

Since — 


(v \’ 


\ u i/ 


V, = 


u, 


1 (1 + *) L 


x cos tf) ± ■yjl - x 2 sin 2 tf) 


(9.52) 


(9.53) 
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m 1 


The ± corresponds to the 2 possibilites in the following figure. Remember, for —- = x > 1 


TC 

Tmax < — => cos ip > 0 So we have 


m Q 


or 



v, = 


u, 


■ (l + x) L 


X cos 




x 2 sin 2 iJj 


t 

short" (x > 1) 


Confirmation: “short” case v 1 should —> 0 as x —> 1 + (through larger values of x): 


v, 


u. 


cos 


y 


sin xjj 


= 0 . 


cosi[>-|cos l[>| 


=> Only the stretch case survives for x < 1. (Can see that for x < 1 the negative root would make 

T 

v ] negative.) The same interpretation applies to the expression for —L . 

T 

x o 

["Aside 2": From ©, 

y equation: 

v 1 sin = v 1 sin 0 

[“Aside 2”: From (T), y equation: 

v 1 sin = v 1 sin 0 

x cos xp ± T “ x2 sin 2 


u. 


sin 0 v x _ (l + x) . 

sin v 1 


u, 


(l + x) 
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"stretch" 

l 


cos (0 - ip) = ±<Jl - x 2 sin 2 iJj . ] 

t 


'' short '' 


[“Aside 3”: From “Aside 1” we have 

sin (0 - i[>) = x sin ip . 


Thus, from 


sin 2 (0 - ip) + cos 2 (0 - x[>) = 1 , 


we have immediately 


"stretch" 

1 _ 

cos (0 - ip) = ±yjl - x 2 sin 2 xjj . ] 

t 


"short" 

We get very simple results when ny = m r Then, for \\f * 0, we have only one 0 value: 



There are actually two solutions where V and v x are co-linear in this case: 



0 f = 0 (two 

particles miss 
each other) 
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The solution where the particles miss each other has ip = 0 whereas the head-on collision is 
characterized by V| / —> 90°. (Remember, ip = — for x = 1. We have for x = 1, 

^ = _1 

^ (i+ 

This goes to zero as V| / —> 90°, which means in this limit ny comes to a complete halt after the 
collision. This fact is useful in nuclear reactor modulators, which slow down (or moderate) neutrons. 
It says that the best way of slowing down free neutrons is a material which contains light nuclei 
(like the deuterium in so-called heavy water.) Obviously, 

^ = sin 2 x\> , ( 9 . 55 ) 

in this limit. 


I)- 


[2 cos ij >] 2 = cos 2 ij> 


( 9 . 54 ) 
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9.7 A KINEMATICAL EXAMPLE IN THE LAB FRAME 

Let s do an example. Let s say a particle of mass ny scatters elastically from one of mass m 2 at rest 
in the Lab frame. The ratio — = f is given. Find the angle \|/ through which m is scattered. 
We have Ul 


v i = 


u, 


(1 + x) L 


X cos 


T ± >/i - 


x 2 sin 2 xp 


Thus, 


=^> f (l + x) = x cos xp ± -^1 - x 2 sin 2 xp 
=> f 2 (l + x) 2 + x 2 cos 2 xp - 2f x (l + x) cos xp 


= 1 - x 2 sin 2 xp 

=> f 2 (l + x) 2 + x 2 - 2f (l + x) x cos xp = 1 . 


Solve for cos xp in terms of x and f: 


cos xp 


x 2 - 1 + f 2 (l + x) 2 
2f x (l + x) 


To find the meaningful range for x, write 


-1 < cos xp < 1 


Plug cos xp = 1 into the above (note 0<f<l and that x is positive): 

f + 1 


x( COS xp = 1) = 


1 - f 


.. (largest x) 


The other limit is for cos xp = - 1, which we get by simply letting f 

=> x(cosxp = 1) = --- (smallest x). 

x(cos xp = -1) 


-f. Thus 
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9.8 RUTHERFORD SCATTERING IN THE LAB FRAME 

Lets find the differential cross sections in the Lab frame. Remember the meaning of da: 
dN i i 

da = — = 2it b db . (9.56) 

I 1 1 

This is unchanged going from one frame to another. (At this stage there is no reference to angles.) 
Then 


b = b (0) , CM frame or b = b (ip) , Lab frame. 


The connection between Lab and CM cross sections is, 


Lab 

i 

da 


dQ 


M = 


CM 

1 

da dQ(6) 

dQ { ’) dQ(ip) 


da 


dQ 

Result of “Aside 1”: 


M = 


da / x sin 0d0 
d£2 sin tpdtp 


sin (0 - ip) d 

---- = x . Do - on both sides: 

sin ip dip 


sm 


(0 - ip) 


COS 


(0 - ip) 


cos ip + 
sin ip sin ip 


d0 

dip 


- 1 


= 0, 


Use 


d0 sin (0 - ip) cos ip 

dip cos (0 - ip) sin ip 


sin (0 - ip) = x sin ip , "Aside 1" again 


cos (0 - ip) = +\Jl - x 2 sin 2 ip , "Aside 3". 

T 

specializing to x < 1 . 


(9.57) 

(9.58) 


(9.59) 
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We find 


d 0 _ n x sin x|x cos x|x 

^ I I 

dty -y/l - x 2 sin 2 xJj sin x|x 


(9.60) 


d 0 _ x cos xp + yjl - x 2 sin 2 xp 
dxp ^/l - x 2 sin 2 xp 

so (x < 1 case; using “Aside 2” now in (9.58): 


(9.61) 


da / x da ,_ x sin0 d0 

-(x|)) = -(0)-, 

dQ v 7 d£2 v 7 sinxjx dxjx 




xcosxjx + yjl - x 2 sin‘ 


xi> 


-y/l - x 2 sin 2 xjx 


(9.62) 
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Also 


from Sln ( 6 ~ M , 


X, 


simp 

0 = ip + sin -1 (x sin ip) 


(9.63) 


Notice for “b” large enough in the x = — = — > 1 case, we are in the “stretch” case. 


m 0 


When b is smaller than the value that causes scattering into angle tp M 


= sin 


f m \ 

m 


V m i/ 


, then 


we are in the “short” case. (Formally, the critical value of the impact parameter is given by 


b 0 = b ip = sin 


m 0 


\ m i/ 


) Both, however, cause scattering into the lab angle y, and thus the 


analog for x > 1 is the rather complicated expression 


dQ^ dQ ( 0 (^)) Stretch 


xcosip + 'Jl-x 2 sin 2 xp 


'Jl-x 2 sin 2 xiJ 


+ ^( 0 (^)) l S hort 


xcosxJj - 'Jl-x 2 sin 2 x|) 


yi-x 2 sin 2 ip 


(9.64) 


If x < 1, can do an expansion in powers of x. To 0— order, 
0 = xp, 


dQ 


w ‘ I ( 0 ) L - 


da 

dQ 


) - - 


’ 4 

\ / 


sm 




(9.65) 


Remember, p is measured in the CM frame. Connection with other quantities: 

Poo = ^ ^initial | = 0 (ui + U 2 ) , 


m i u i = m 2 u 2 


m. 


u 2 = —- u 1 , 
m 0 


= 


m 1 + m 2 
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m.m, (, m, 
——— 1 + — 

V m 2 ) 


nq + m 2 


iq = m 1 u 1 = jaxq 

t 


u i = 7 


u, 


m 2 / 


Now 


2 2 2 2 

P» _ 1* U, mq 


2\x 2\x 2 


= T 0 , 


But 


m 2 

T = - — 

■ L o 

m 1 + m 2 


T = 


1 + x 


Therefore for x « 1 we have 


da 


(^) - — 
dQ K ' 16 


k V 1 


Vh)/ sin 4 — 


( 9 . 66 ) 


written entirely in terms of Lab quantities. Now let’s see what the first order correction to this 
expression (in x) is. 

0 = 0 O + X0 : + ... , 

(9.63) => 0 O + x0 1 = sin -1 (x sin xjx) + x|), 

• -i y 3 3 5 iii 

sin y = y + — + — y + ..., y < 1, 

6 40 11 


0 O + x0 : ~ x sin x|x + x|x , 


0 O = xp 
0 1 = sin x|x. 
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Also 


x cos xp + -yjl - x 2 sin 2 xp 


yjl - x 2 sin 2 xp 


[l + X COS Xp ] 2 


sin — - sm 
2 


xp X 


4 x xp 

~ sin — + — sin xp cos — , 

2 2 2 


^P_2 ^P 


= sin — + x sm — cos 


= sin 


1 + x cos 


2 2 

2 


■ 4 0 . 4 Ip f, „ 2 'Vp 

sin — - sin — 1 + 4x cos — 


da 

dQ 


M " \ 

^pik^ 

2 

1 

(l + 2 x cos xp) 

Ip-J 

. 4 IP / 

sin — 

2 

1 + 4 x cos 2 ^ 


~ 1 + x 


1 

/ „ 2 xp^ 

2 cos xp - 4 cos — 


1 + x | -2 cos 2 — - 2 sin 2 — I = 1 - 2x . 


jxk 

_2 




V 2 M-/ 


k k (l + x) 


2T n 


2T„ 


P» , 


k ^ 


\ 2T o/ 


(l + 2x) . 
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Finally, then to first order in x: 


da 

dQ 


M 


i 

16 


\ T 0/ 


• 4 ^ 

sin — 


(l + 2x) (l - 2x) . 


-1 


First order correction vanishes! This means the Rutherford cross section, which is strictly an x=0 
result, holds very accurately for 0 < x < 1. Will not torture you with the next order, -x2, correction. 


9.9 TOTAL CROSS SECTION 

Total cross section: 


a 


f do 

all 

angles 


da 

dQ 


(0) dQ . 


(dQ = sin 0d0d<|)) 


(9.67) 
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For Rutherford: 


a 


1 ( [ik'l 

MpTJ 


„ . e e 

2 sin — cos 
2 2 

1 


2n TC 

J d + / 


0 0 + 


sin 0 d 0 


sm 


0 • 

2 


(strictly speaking, 
the point 0 = 0 is 
excluded in the integration) 


a 


2jr 


(p-j 


f 

n + 


0 

cos — 
_2_ 

• 3 0 
sin — 


d 0 

2 


1 


a 


3k 


\P«y 


1 

2 sin 2 



+ oo. 


Understandable: particles are always deflected regardless of b value. Can see it in: 


da = 2 jt bdb , 

k max 

^ o = 2k f bdb = Jtb 2 ax . ( 9 . 68 ) 

0 

where b max is the maximum impact parameter which suffers an angular deflection ^ 0. Thus, in 

classical mechanics, the only type of force laws for which a is finite are those of the form 

F(r) = 0 , r > a 


where a is some finite value of separation. Example of this possibility (0 = \|/ here since the sphere 
is infinitely heavy): 



radius = a 
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Clearly, this problem has azimuthal symmetry if we take the +z axis as shown. We have 

b = a sin 0 , 

20 + 0 = tc . 


notice: -sign 
(repulsive scattering) 

i 

e 


b = a sin 


jt - 0 

= a cos — 
2 j 2 


da 

dQ 


sin 0 


db 

d0 


db 

d0 


=-sm 


da 

dQ 


0 

a cos — 
2 


0 


2 sin — cos — 
2 2 


a . 0 

- • - sin - 
0 2 2 


-a 2 

4 


a = f — a 2 d£2 = Jta 2 

J 4 


Just what we expect. 


Download free eBooks at bookboon.com 


93 











MODERN INTRODUCTORY MECHANICS PART II 


SCATTERING AND COLLISIONS OF PARTICLES 


9.10 PROBLEMS 

1 . The Rutherford differential cross section is 


da 

dQ 


1 (jjJc'' 

4 \plj 


. 4 0‘ 
sm — 

2 


Therefore, the backscattering cross section is 


a. 


Vk 

\ P °°) 


/ 2n 

, d */ de 


jt/2 


sin 0 

. 4 0 
sin — 

2 


In the text I estimated the integral by setting -yy |e-*= 1 < which resulted in an 


e 

sin — 
2 


underestimate of a back . Now do this integral exactly. By what factor was the original 
estimate off? Will this improve the agreement with Rutherford s experiment? 
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2 . We estimated in the notes the fraction of 5 MeV alpha particles backscattered (- < 0 < jt) 

from a target made of gold foil and found that about one in every 21,000 particles should 
backscatter. (This was changed to about one in every 10,500 after we corrected the integral.) 
Suppose Rutherford wanted to make a lead shield to protect the other experiments 
in his lab. Note that for lead, Z = 82 and m pb = 207 amu, and that the density is 
p pb =11.4 . Estimate the minimum thickness of the lead that would shield against 

these alpha particles (in centimeters). [Hint: Require all particles to be backscattered from 
the lead. Neglect the possibility of multiple scattering.] 

i 

3. Given that (T 0 is total CM energy) 

k 1 


exactly for some unknown central force law, find 

a) the CM differential cross section, — (0). 

dQ jt 

b) the number of particles backscattered (that is, with angles 0 such that — < 0 < jt ) in 
the CM frame. (Call this N back ). Assume a known incoming particle flux, I. 

4. Do the integration in (9.16) leading to (9.17) of the notes. 



5. Consider scattering off of a weak potential U(r) such that U(b)« E, where b is the impact 
parameter and E is the total energy. Show that 


1 + 


U(b)\ 
2E j ' 


which shows that r min < b for an attractive potential (U(b)<0)and vice versa, as one would 
expect. 


6. Look up and plug values in (M s , R s are the Suns mass, radius) 


A0 = , 


to get the angular deflection for starlight in seconds of arc. 
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7. Consider a head-on collision as seen from an unknown reference frame. Assume the ratio 

m i 

x ~~ is known. It is observed that ny comes to a complete stop after the interaction. 

Assuming energy and momentum conservation, find the after/before ratio of ny’s kinetic 
energies in terms of x. 


8 . Derive: 


(a) tan iji = 

(b) sin ip = 


sin 21 


x - cos 2§ 
sin (() 

-y/l + X 2 - 2X COS (|) 


where x = —- and the angles \|/ and (|) are defined in the text. 


m 0 


9. From prob. 9.8(b), or any other means, show that 


tan \\) 


sin (|) 
x - cos (|) 


10 . In reference to Eq.(9.65) of the text, evaluate the critical impact parameter variable, b , for 

k 

the Coulomb scattering potential, u (r) = ± —, in terms of k, p, and x. 

r 

11 . Consider scattering of a point mass ml off of a hard sphere of radius a and mass m 2 . 
(The ratio x = is arbitrary.) 



radius = a 


The angle of the sphere’s (m2 s) deflection in the laboratory frame (m2 initially stationary) 

i /1>\ 

is given by , § = sin 1 — ,where b is the impact parameter. 

\a/ 

a) Show that the cross section evaluated in the CM frame of reference is a constant. 

b) Find the deflection angles 0 and \|/ for ml if the impact parameter is b = -^=. 

V 2 


In addition to the above, find the the cross section in the Lab frame of reference for c) 
x< 1 and d) x> 1. 
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12 . Show that (9.62) of the notes can be written more simply as 


da 


da 


/v V 


— (i|>) = — (0) 1 

dQ dQ 


\ v i/ 


v i v i 

Vi • v' 


13. Given a center of mass differential cross section, 


= A cos 2 0 , 
dQ 


(0 is the deflection angle of ml in the CM frame) and a particle flux, I, find the number 
of backscattered particles ( — < ip < jt ) per unit time in the Lab frame. Assume m 2 
» m 1 where m l is the mass of the incident particle. 

14. With the same cross section as in problem 9.13 above, find the full Lab cross section for 
x < 1. How does one specify the “short” case in the relation between 0 and \|/? 
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15. Hard sphere scattering: 


incoming sphere, radius = c 



sphere, radius = a 


Assuming a symmetrical scattering event, find and the total cross section. 


Other Problems 

16. Escape velocity is the minimum speed a particle needs to escape a planet or star, starting from its 
surface. Previously, we estimated the angular deflection of a light beam traveling near the Sun, 
treating the light as if it were an ordinary massive particle. In the same spirit, find a formula 
for the maximum radius of a star of mass M from which light, traveling at the speed of light, 
c, may no longer escape. (This is called the “Schwartzchild radius” of the star.) 

17. The maximum scattering angle, V|/ max , in the Lab frame for ny in the x > 1 case 

(x — —) was given by sin(ip max ) = —. Show that this happens when the CM 

m 2 x 

scattering angle, 0, satisfies, 

p 1 
cos 0 =-. 

X 


18. Consider a finite range inverse square force law. That is, we have an attractive force 


F( r ) 


f k 

- — , r < r 0 
1 r 

[ 0 , r > r 0 
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Show that the deflection angle 0 is given by (E=total energy,b=impact parameter) 


2 b 2 - ^ 


cos 0 = 


r„ .11 — 1 + 4b 2 


1 - 


r o E / 


19. Concerning the same problem as 9.18, show that the CM scatterimg angle 9 can be written as 


0 = 2 sin 1 


'b' 


\ r o/ 


- sin 1 (cos 0) 


20 . Given the results in probs. 9.18 and 9.19, go as far as possible in finding the scattering 

cross section, -, in the CM frame for this force law. 

dQ 

21. Eq.(9.62) gives the scattering cross section in the Lab frame, —- (tjj) | Lab in terms of ,-^g 

dQ 

cross section in the CM frame, ( 0 (r|j) ) | CM , for the x = m /m < 1 case. Turn this around 

dQ 

by expressing /qx i in terms of the laboratory quantity, also for x<l. [Ans.: 
dQ CM 


da 
— ( 9 ) 
dQ 


CM 


da 

dQ 


(T(0)) 


(1 + X cos 0) 


Lab [1 + x 2 + 2x cos 0] 3 7 2 ' 
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where 

sin 0 

tan ip = -. ] 

cos 0 + x 

22 . (a) Show directly that the equations (x = m /m ) 

0 = ip + sin _1 (x sin ip), 

and 

cos 0 = -x sin 2 ip ± x cos ~ sin 2 ip , 

are compatible. 

(b) For x = 2, what two roots, 0 ,0 2 , are associated with \\f = tt/8 ? 
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23. Given the result for the Rutherford differential cross section in the CM frame for m , 


— ( 0 ) = - 

dQ v ’ 4 


[ik ^ 

VP- / 


. 4 0 

sin — 
2 


find the differential cross section for m 2 (NOT ny) as a function of the appropriate angle: 

a) in the center of mass (CM) frame 

b) in the Laboratory frame (m initially at rest) 

24. Given the result for the Rutherford differential cross section in the CM frame (all for ny), 

da /^v 1 (iik\ 1 


(xjj), in the Laboratory frame, when ny = m 2 (x = 1). 


dQ 4 




. 4 0 

sin — 

2 


da 

find the differential cross section, - 

dQ 
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10 NON INERTIAL 

REFERENCE FRAMES 

10.1 FINITE DISPLACEMENTS AND ROTATIONS 

Need transformation connecting fixed and moving (body) noninertial frames. Picture: 



infinitesimal amount in for the r axes, but 


time, dt finite displacements, R 


(1,1’), (2,2’), (3,3’) axes coincide in direction at some instant of time, t. Clearly, the reason for this 
is to describe, for example, motions relative to the Earth, which is noninertial. This is all done for 
convenience, not any real physics reason. In fact, strictly speaking, this chapter has zero physics 
content! This does not mean, however, that these considerations are not useful or convenient. 


There are two types of transformations that will be involved: 


1. Displacement: r = r' - R . (see above figure) 

2. Rotation: (a “generic” passive rotation) 


r, = 



(passive) 
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Picture: (specialized to rotation about 3,3’) 
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Put them together. Step 1: displacement. At the end of step 1: 



Now wish to rotate. Rotate the r axes: 



A rotation 
about an 
arbitrary axis 


r, 



But 


r j = 


r j - R : r SO 


r, 



( 10 . 1 ) 


( 10 . 2 ) 


Download free eBooks at bookboon.com 


104 








MODERN INTRODUCTORY MECHANICS PART II 


NON INERTIAL REFERENCE FRAMES 


Represents: 



In matrix notation, this is 
r = X(r' - r) . 
It’s inverse is 


r R = X _1 r = X T r . 


(10.3) 


(10.4) 


We really need only the relationships between r and r ' for an infinitesimal rotation. What is X 
for such a situation? Remember 


= 1 , X. -1 = X T , 

2 ^ik = 8 jk * 

i 

Assume 




t 


change in X necessary 
to represent an infinitesimal 
passive rotation 


(10.5) 


( 10 . 6 ) 
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Substitute above: 

2(8 1) +6^ i )(S ll + 6 lk )-8 Jl . 

0- order: J Mi* * 6 u • ^ 

i 

1- order: J ( 6 ij 8 A k + ^iA*) = 0 ' 

i 

=> 5X ±j + 5X kj = 0 . 


or 8A, jk = -SA. ki , antisymmetric 


v kj 


Also implies there are only 3 independent elements: 


bK, = 


ID 


0 (1) (2) x 
-(1) 0 (3) 

-(2) -(3) 0 


(1), (2), (3) arbitrary elements 


(10.7) 


( 10 . 8 ) 
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10.2 INSTANTANEOUS RELATIONS FOR VELOCITY, ACCELERATION 

Now go back to Ch. 1. Representation of an infinitesimal active rotation on a vector: 


50 



or 8r i = S E i* 80 j r k • (10 - 9) 

Now, any active rotation is given by a passive rotation in the opposite direction. Start: 



Passive, - 0: 



Active, 0: 
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Passive , - 0 


have the same 
components 

\ 



This gives us two ways of specifying the effects of an active rotation: 
Way 1: 6r ± = ^ £ ijk^0j r k • (active rotation) 

Way 2: 6r ± = ^ X T ±j Tj - r ± , (passive inverse rotation) 

j 

(final) (initial) 



Must be the same: 

2 8XT ik r k = 2 e ijk 80 j r k' 

k j, k 



or =* 6X ki = J £ kij 80 j • 
j 

Notice, as expected, 6X ki is antisymmetric in k,i. 

Our relationship between the primed and unprimed coordinates are again, 
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Consider an infinitesimal change on both sides: 


© © 

6r i - - y 8*1, r , + 2 '•!; 5r, - (10.13) 

j j 

Why 2 terms on right side? (Q arises from the rotation of the noninertial axes while (2^ is due 
to the independent motion of the particle relative to the ? axes. Found earlier, 

- 2 6e * • 

k 

Of course also 

(10.14) 
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SO 



2 £ i*D 80 k r D + 2 ( 8 U + bX %) br i • 


(10.15) 


t 


can drop (2 nd order in small quantities) 


Thus 


8r i - SR, = J £ ikj 80 k r j + 8r i * 

K j 


Divide by 8t 


'8Q 

v6t 


dQ \ # 

dt / * 


dr! dR ± _ dr ± y d9 k 

dt dt dt p. ikj dt j ' 

D 


But 


co k 


de, 

dt 


so 


dr ' dr dR 

- = — + — + co x r , 

dt dt dt 


or 


(10.16) 


(10.17) 


(10.18) 


(10.19) 


v f =v r + V + (oxr. (10.20) 

v f j is the velocity of the particle relative to the fixed frame. (It is the velocity as measured 

by an observer at rest in r') Must be a constant in magnitude and direction if the particle has no 
real forces acting on it. v r is the velocity relative to the moving frame, whose axes coincide with 
the fixed aces at the given instant in time. 
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Picture: 



Now do the second variation: 


Sr: - SR, - 2 8 K, r, ^ Sr, , 


6 2 r : - 5 2 r, = 


2 6 2 E, + 2^ 8r, 

j 

+ 6 ! r,. 

j 

t 

can replace by 5j_j 


( 10 . 21 ) 


Before, we compared 

6r i = 


to 


and got 


6r ± = J SX T ik r k f 


6^ T ik = 2 £ ijk^0j 


( 10 . 22 ) 


(10.23) 


(10.24) 
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Now 


But 


so 


S 2 ^i = 2 8 ij k 6(60^) 




= 2 8 i* ( 8 V* + 8e : r 0 




8r k = J 8 ^m 80 f r m r 


£,ia 


82r i = 2 ^ 820 :^ + 2 8 W* 80 j 80 ^ 


jfk 


£,m 


[let m->k) 


(10.25) 


(10.26) 


= 2 8 i3* 820 3 r k + 2 ^m^k 80 ; 80 ^ r k 


j,k 


jfkf 

l, m 


(10.27) 


H 


w 


T. 
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On the other hand, compare this to 

6 2r ± = 2 r k • 

k 

Identify: 

J 6 ' e ) + I 6e = 6e ' • 

j 

(Switch indices: k —> j , j —> k everywhere.) Put it all back together: 

6 2 r ; - b 2 R, = 2 EikjSXrj + 2 ^km^j^k^Vj 

j/k, 

£ r m 

+ 2 2 e ikj 8 0k 8r j + S 2 r ± • 


j/k 


(10.28) 


(10.29) 


j/ k 


(10.30) 


Divide by dt 2 


( b l Q d 2 Q 


5t 2 


dt 2 


dr d 2 R ; 


dt 2 


dt 


m2 


e... + d r ‘ 


k /j 


~ ikj dt 2 


^ikm 


k,m 


dO, 

dt 




d0, 

dt 


dt 

\ 


r. 


/ 


_ V d0 k dr i 

+ 2 2 8 ikj 

j/k 


dt dt 


T , ., d0 ± . d 2 0 ± , . . 

Identity co = —- , co. = -- and write in vector notation: 

1 dt 1 dt 2 

r' = R + f + o5xr + cox(coxf) + 2coxf. 


(10.31) 


(10.32) 


R represents the acceleration of the origin of the r coordinate system relative to the r' origin. 
Will be zero if we consider uniform motion. Picture: 



Also have co = 0 in the case of constant angular velocity (magnitude and direction). 
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10.3 USEFUL EARTH COORDINATE CHOICES 

Write in the inertial system 

F = mr. 

Write in the noninertial system 

F ef f = m r- 

which gives 

F eff = m |rR-(oxf-ajx(toxf)-2wxfJ, 
or 

F eff = F-m|R + coxf-(hx(coxf) + 2coxfj. 

t t 

"centifugal" "Coriolis" 


(10.33) 


(10.34) 


(10.35) 


Remember: 

Deflection is to the right in 
northern hemisphere and 
to the left in the southern 
(relative to the initial direction) 


Leads to deflection of air 
masses in a counter-clockwise 
direction in the northern 
hemisphere 
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For the Earth, we often choose: 
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For this choice co = R = 0 , so 

F eff = F - m |(o x (w x r) + 2co x r j • 

Sum up. Because instantaneously our axes coincide in direction, we have 

r' = r + R , 

f'=r + R+ coxf, 

r' = f + R + ioxf + 2(oxf + (ox((oxf). 


(10.36) 


How come we do not get r' as simply the time derivative of the expression r' for example? Because 
?' and r are referred to different axes which are rotating as well as moving with velocity R with 
respect to one another. There is another way of viewing this process more in line with the books 
derivation. When a change in r is considered, calculated with respect to the moving axes, we have 


(6?) f = (6r) r - (6f) passi 

t t 

fixed rotating 

From before: 


(10.37) 


understood in fixed frame 


(Sr) 


passive 


-(Sr) 


active 


-80 x r, 


=> (8r) f = (8r) r + 80 x r . 


Thus 


fdrx 

\dt) f 


(drx 

l dt/ r 


+ co x r. 


/ 

CO 

\ 



(10.38) 


(10.39) 
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This is true for any vector, not just r, as long as the fixed and rotating axes instantaneously coincide 
in direction. Including the effect of translation of the coordinate origin, this now gives 

^dr\ /'drx dR 

Vdt/ f ldt/ f dt 

t t 

make understood in fixed 

replacement frame 

=>r'=r + R + cbxr as before. 


Apply the same reasoning to get r': 


(di') 


(di\ 

dt ) 

f 

Ut] 


dR 

+ — + co x r + to x 
dt 


see above 

1 

'dr\ 

vdtj f ’ 


However 


((di\ _ ^ 

— I = —I + co x r, 

Vdtf V dt/ 


r' = r + coxr + R + coxr + coxr + coxr 


|r + co x r), 


(10.40) 


(10.41) 


r' = r + R + cox 


r + 2co x r + co x (co x r). 


This is also as before. 

Get back to F e f f on Earth (static case, r = 0) 

tff = F - mco x (co x r) . (10.42) 

This means that the effective acceleration due to gravity and the Earth is given by 

§eff = § - CO X (co X R). ( 10 .43) 

t t 

gravity R points from center 

only of Earth to surface 
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Picture: 



As we can see, -co x (to x r) has components along 3 and 1 axes. Another view (corresponds 
to the 13 plane above): 


|co 



has a small southerly 
component in Northern 
Hemisphere 
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At equator, its particularly simple: 


too 


- CO X 


CO X 


*) 



Center of 
Earth 


You will find the angular deviation of a plumb line from the true vertical caused by this effect 
in a problem. 


Exampl e: central force problem. 
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dU - 

F --e, , R = 0 

dr 1 


o 


mr = F - m 


$ + (oxr + (Dx((Dxr) + 2(Dxr 


co = -toe 3 , co = -toe 3 , 

r = re 3 , r = re 3 , ^e 1 = 0 in rotating framej 

r = re lf 

63 x r = -tbre 2 , 

263 x r = -2tore 2 , 

65 x (63 x r) = -w 2 re 1 . 


‘2” components: 

0 = -m (-6)r - 2cor) . 


2 £ £ 

integrate => r to = — , to = -- , as before (conservation of angular momentum). “1” 

„ m mr 

components: 

dU / 2 x 

mr -- - m -to r , 

dr v ’ 


2 1 dU 

r - to r =- 


m dr 


or r - 


2 3 

mr 


- , also as before. 

dr 


Download free eBooks at bookboon.com 


120 









MODERN INTRODUCTORY MECHANICS PART II 


NON INERTIAL REFERENCE FRAMES 


There are other choices of noninertial coordinate systems which can simply motion problems near 
the Earth’s surface. Consider the choice: 



For this choice we have 


F eff = F - m 


R + Q&xr + cox(coxr) + 2coxr 


We need to compute R . Using the general formula 


(—) 

\dt/ f 



+ to x A , 


for any A, we get 



CO x R , 


(10.44) 


since 


dR 
dt, 


= 0 . Again applying the above general equation, we get 


! d 2 R^ 


dt 


= co x 


\ UL It 


/ dR^| 

\dt J f 


= co x 


(co x r) 


Notice that we now get 


usually small 


F eff = F - mco x 


1 

(to x r) - mco x (co x r) - 2mco x r 


(10.45) 
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10.4 DEFLECTION OF PROJECTILES NEAR EARTH'S SURFACE 

Using this F eff , we can now investigate the motion of projectiles near the Earths surface. 
Qualitatively: 


3 (local vertical) 



South 


CO 


7.29 x 


icr 5 


radians 

second 


co = -co cos X e 1 + co sin X e 3 . 
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When projected initially to the East, the particle trajectory will gain a Southerly component due 
to 1. centrigulal force and 2. Coriolis force. However, if initially projected West, the force will 
now deflect the particle to the North. The interesting question is: which deflection will be bigger 
when projected to the West? 

Take 


ef f 


- F - m 


(o) x (co x r) + 2co x r). 


To 0“ order, we have 


R = R E e 3 , 

t 

radius of the Earth 


r i(t) = 0, 

r 2 (t) = - (v 0 cos a) t (to the West if v 0 > 0) , 

r 3 (t) = - — gt 2 + (v 0 sin a) t (does not include effect of rotation) 


Work out -co x (co x r) term: 


co x 


R = (—co cos Xe 1 + co sin ke 3 ) x (R E e 3 ), 

= coR e cos Xe 2 , 

-co x (co x r) = -(-co cos + co sinA,e 3 ) x (coR E cos A,e 2 ) , 


= co 2 R e cos 2 Xe 3 + co 2 R e sin X cos Xe 3 . 

t 

term we are interested in 


Work out Coriolis term: 


2co x r - -2 (—co cos A,e 1 + co sin Xe 3 ) 

(- 


x i — v 0 cos ae 2 


(-gt 


+ v n sm a 


)® 3 ) 
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- - 2(0 cos X v 0 cos ae 3 -2<o cos X (-gt + v 0 sin a)e 2 


-2co sin X v 0 cos ae x 

t 

This is the term 
we are interested in 


Plugging these results back into our F eff equation, we find 

(tff)i = m t - co 2 R e sin X cos X - 2(0 sin X v 0 cos a . 


If our initial condition is that r (t) = 0 at t = 0, then 

r x (t) - (<o 2 R e sin X cos X - 2<o sin X v 0 cos a)t 2 


Now eliminate the time, t, by using the 0~ order equation, 


3 (t) = 0 = - j gt 2 + (v 0 sin a)t, 


t - 


2v„ sin a 


This is the approximate time it takes for the projectile to hit the ground. Therefore 


r^t) - — (<o 2 R e sin X cos X - 2oo sin X v 0 cos a) 


4v 2 sin 2 a 


Question: which effect is larger? Depends on initial velocity, v Q . “Break even” velocity is 


/ \ 0 )R„ cos X 

h.L — 1 - 


2 cos a 


( ooR p _ meters 
1 - £ - 232 - 


l 2 


sec 


Depends on latitude, cos a. (In this problem I have not been very careful about taking care of the 
Earth’s curvature. The above considerations only hold for short range projectiles.) 
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10.5 DEFLECTIONS FOR DROPPED OBJECTS 


Let me now introduce another coordinate system which is useful for calculating the deflection of 
projectiles relative to the gravitational vertical. Remember: 
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Take new, skewed axes along g eff : 


Then since 


3 



tangent to surface 


(2 axis out 
of page) 


F eff = F - mco x (co x r) - m(o) x (oo x r) + 2co x r), 


mg 6 


= ma, 


eff r 


=> a eff = g eff - to x (to x r) - 2o) x r . 


(10.46) 


where now 9 e ±f = 9 e ff e 3 only. This is useful in discussing the deflection of particles relative 
to the local gravitational vertical, which can be established with a plum-bob, say. For example, if 
we had done the projectile problem above in the “skewed” frame, the term proportional to R E in 
the form for r^t) would have been absent. Then, to first order in co: 


a eff - -9eff e 3 - 2co X r . 


(10.47) 


In this new frame, we need to find co . We had 
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After we “skew” it by an angle, it looks like (rotation is around 2 axis): 



Clearly, we have 


co 1 = - 0 ) cos (X + e) , 
< co 2 = 0, 

co 3 = co sin (X + e) . 


SIMPLY CLEVER 
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However, 8 is a small angle (smaller than .002 radians) and will be neglected below. Now imagine 
dropping an object from a height, h. 


In 0 th order: 


r 3 = ~9eff' *1,2 = 0/ 


r = -g eff te 3 , r = - i g eff t 2 e 3 (if r = 0 at t = 0) 


1 st order (in to) correction: 

co 1 - -co cos X, 

co 2 = 0, 

co 3 - co sin X. 


as explained above 


Then, 


-2co x r - 2 (-co cos Xe 3 + co sin Xe 3 ) x (-g eff te 3 ), 
= 2co cos Xg eff + e 2 . 


We now get 


r 3 = -9eff > *i = 0, r 2 - 2 cog eff t cos X 


Integrating twice on r 2 , we get 


r 2 (t) - ^ tog eff t 3 cos X 


Of course, we have 


9eff 

=> r 7 (t) - — co 
3 

Since r 2 > 0, the deflection is to the East. You will study this problem further (to second order in 
oo 2 ) in a further HW problem. 


3 \ 


8h. 

cr ~ 


cos \ 
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10.6 FOCAULT PENDULUM 

Last problem: the Foucault pendulum. Again, use our “skewed” coordinate system. Only changes: 
g F eff = F + mg ef£ - m [to x (to x r) + 2to x rj g eff 

t t 

new external force small (ignore) 

Say 

F = F 1 + F 2 e 2 + F 3 e 3 , 

9eff = — 9 e ff ' 

to = -to cos X e x + to sin X e 3 . 


Put it all together: 

miq = F x + 2mtor 2 sin X , 

mr 2 = F 2 - 2mto (r 2 sin X + r 3 cos X), 

mr 3 = F 3 - mg eff + 2mtor 2 cos X . 

Consider: 



Fi 




F 3 = 


Lllil 

l 


T. 
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Neglect r 3 and 


~ 0 


0 rd n 

3 eq 


above becomes 


=> T - mg eff - 2mcor 2 cos X . 

small but * 0 
(not a constant) 


Substitute this value of T above: 


mr 


y- [2mcof 2 cos 




+ 2mcor 2 sin A , 


t 

small 


-r. 


mr 0 = 


[2mcor 2 cos A] - 2moo[r 1 sin A + r 3 cos A] 


t 


t 


small 


small 
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Set a 2 = ^ eff . Then 

£ 

r x + a 2 ^ = 2<or 2 sin X , 
r 2 + a 2 r 2 = -2<s>± 1 sin X . 

There is a dynamic coupling between the motions. Can solve both at once using complex numbers: 

(r : + ir 2 ) + a 2 (iq + ir 2 ) = 2(0 (r 2 - iq) sin X , 

= -2coi (ij + ir 2 ) sin X , 


or 


q + a 2 q = -2coiq sin X . 

where q = + ir 2 . This is solved by assuming 

q = A e u , 

where A, X are to be determined. Substitute above: 
=> A. 2 + 2i<x>A. sin X + a 2 = 0 , 


X = 


2 L 


-2io) sin X ± x/-4co 2 sin 2 X - 4a 2 


X = -ico sin X ± ix/co 2 sin 2 X + a 2 


But to 2 sin 2 X « a 2 




= — for the Earth, so 

l) 


X - -ico sin X ± ia . 


General solution: (A,B real) 

g(t) = {Ae iat + Be - iat } e' icotsinX 
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A,B are fixed by initial conditions. Write it out: 

A (cos at + i sin at) + B (cos at - i sin at) = C 1 cos at + iC 2 sin at. 
Let C x = A + B , C 2 = A - B , 

=> q(t) = (c x cos at + iC 2 sin at) e~ ia)tsinX . 


There are 2 parts to the motion with different frequencies since co « a . Let ' s say e i(0tsin ^ = i. 
Then, since 

r ± = Re q(t) , r 2 = Im q(t) 

= C x cos at = C 2 cos at 



Eq 11 of an ellipse, fiowever, we made r 2 and r 2 coordinates in the complex plane. The factor e 1C0t sin x 
is just a rotation in the complex plane: 



The real and imaginary parts of the complex number g(t) are just the r l5 r 2 components of the 
real motion. Therefore, a rotation in the complex plane is also describing a rotation of the vector 
r = e 1 + r 2 e 2 in real space. Because of the minus sign in e “ ia)t sin ^ this [ s a clockwise 
rotation in coordinate space. (It would be counter-clockwise in the Southern hemisphere.) Time 
it takes to complete a rotation: 


2jc = coT sin X 
2jt 

co = - => T = 


T = 


2jc 


co sin X 


1 day 


sin X 


days . 
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Goes around once a day at 1 = 90° (N or S poles) and does not precess at all at the equator (1 = 0°). 
Actual motion looks like: 



Looks like the precession of the orbit of a planet under general relativity, but the forces here certainly 
are not central. The “force” that makes it precess, in fact, is purely fictional. 
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10.7 PROBLEMS 

1. (a) I showed in the notes that instantaneously 

r ' = r + R , 
r ' = r + R + coxr, 

r ' = r + R + coxr + 2(oxf + cox((oxf). 


Find the relation between r ' = a ' and r = a. (R = A); assume all derivatives of 
w to be zero.) Ans: 

a' = a + A + 2o5xa + 3o5x(o5xf) + cox(cox(coxf)). 

(b) Continue this process and find the relation between a ' and a. 

2. Show that the angular deviation 8 of a plumb line from the true vertical at a point on the 
Earths surface at a latitude X is 

r n o) 2 sin X cos X 
8 « —- r 

g - r 0 co 2 cos 2 X 

where r Q is the Earths radius and g is acceleration due to gravity. 

3. By balancing centrifical “force” and gravitational force, find the orbital velocity of an object 
in a circular orbit just above the Moons surface. 

= 1-74 x 10 s cm, M = 7.35 * 10 25 gm. 


G = 6.67 x 1(T cm 


dyne x cm , 

2 

gm 

\ Newtons gravitational constant 


4. At latitude, X = 3E3°, how many hours does it take for the plane of a Focault pendulum 
to complete a revolution? 
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11 RIGID BODY MOTION 

11.1 CONCEPT OF A RIGID BODY 

Attention up to now has been focused on individual point particles or collections of point particles 
for the most part. Not very realistic. 

Consider the example: 



into page 


Now have to deal with the disk, not as a collection of individual particles, but as a whole, i.e., a 
rigid body. For the disk, 

» 1 2 
10 = N , I = - mR 
2 2 

t t 

like: ma z = F z like : mass 

“I” is a way of representing a rigid body in a simple manner. Kinetic energy is, 

1 . 2 1 *2 

T = - my + - 10 , 

2 2 

separated into rotational and translational parts. First, let s understand this result for kinetic energy 
and then derive the form for " I". 
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11.2 INSTANTANEOUS KINETIC ENERGY IN BODY FRAME 

Picture: 


3 



The r axes are considered fixed to the body. 
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Connection from last Chapter: 


r'=r + R + coxr. 

t t 

measured components measured 
in "fixed" in body frame 

frame 


(ll.l) 


We are going to be describing the motion from the body frame. (Why not just describe it from 
the “fixed” axes? Wait a sec.) Clearly we have 

r = 0, (r = 0, r = 0,...) (11.2) 

since we are keeping every point in the body fixed relative to body axis. In order to be able to deal 
with the motion of rigid bodies it is crucial that both the kinetic and potential energies separate 
into center of mass and relative motion pieces. We have already seen this is true for a collection 
of point particles or planets interacting via forces derived from potentials. Relative to the fixed 
frame we have 


(11.3) 

(11.4) 

Using Eq.(ll.l), we have 

T = ^ MR 2 + R • J m «(“ x U + \ 2 m «(“ x r a ) 2 . 

^ a ^ a 


R 


2 mara ' (usually use r a for fixed frame) 

-Si 


_1 - 2 


m « r a ' 


But if we choose the origin of the body system to be at the center of mass, we have 


2 m « r « = °' 


1 — 2 

T = - MR 

2 


+ 9 2 Itla ^ t0 X r d * (separated) 


(11.5) 
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Actually, there is another possibility for describing rigid body motion: if R = 0 for some point 
other than the center of mass. This means the origins of the fixed and body systems are not in 
relative motion. We can then, with no loss in generality, assume the two origins coincide. It is 
important to realize what R = 0 means. Since the r axes are fixed in the body, this is a requirement 
that at least one point associated with the body is not moving in some inertial frame. And, if we 
want this description to be valid over finite time intervals, that point must not move (equivalent 
to R = 0, R = 0, ... ). For example, one can describe in this case a top with one point fixed: 


co 



Can not describe a situation where the tip of the top is moved around arbitrarily, however. 

11.3 ANGULAR MOMENTUM AND THE INERTIA TENSOR 

Axes situation for an arbitrary spinning object: 


3 ' ,3 to 
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Point: instantaneous axis of rotation passes through the center of mass if there is no external torque. 
An instantaneous axis of rotation always exists for a rigid body. Sometimes given, sometimes it 
must be deduced. We will learn later that even if L = 0 , this does not imply o5 = 0 for a rigid 
body in general. 

We know that 


r a = <*> X r a . 

The angular momentum associated with this mass element is 
L« = f a x p a = m a f a x f a = m a [r a x (to x 
Sum on a and use a vector identity, 

L s 2 = 2 m « [ r «“ ■ • 


a 


a 



( 11 . 6 ) 


(11.7) 


( 11 . 8 ) 
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Write this out explicily in terms of components. Useful notation (using “x” instead of “r” from 
now on in this chapter): 

x ^i component label 
particle label 

Then 



= 2 m « 2 ^ 8 ^ 2 


X ak " X ai X aj 


= J 2 ^ 2 " X ai X aj 


Define (“Inertia tensor”; more on this later) 


X ij = 2 m « 

a 

Then, in component notation (notice I = I..), 

L i = 2 (0 : ' 

j 

or in matrix notation, 


8 . 


ID 


x„ 


- x„. X 


aj 


L = I w <=— column matrix 

t t 

column 3X3 square matrix 
matrix 


(11.9) 


( 11 . 10 ) 


( 11 . 11 ) 


( 11 . 12 ) 


Terminology: 

1.. : “moments of inertia” 

a 

1.. (i*j): “products of inertia” 
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Because I.. for i * j is nonzero in general, can now understand the above comment that L = 0 does 
not in general imply. Of course, in the usual case that one considers the object to be a continuous 
distribution of mass rather than a collection of discrete elements, one has 

^ m a -» J dv p(r) , so that 

a . V 

f 

(volume element) 


I ij 


f dV P ( f ) 


Si, 


2 x * ~ x i x : • 


(11.13) 


Do a cylinder as an example: (density = p = constant. I usually use p , z, (|) as the cylindrical 
coordinates, but because p is being used as density, I will use the set r, z, § instead.) 



I 


33 


f rdrd(|)dz p 

dv in 

cylindrical 

coord. 


(x 2 + y 2 + z 2 j 

V V 

x 2 +y 2 =r 2 


r 


'~33 


— R 4 2rttp. But p 
4 


m 

KR 2 t ' 


1 2 
- mR 2 . 

2 


Download free eBooks at bookboon.com 


141 
















MODERN INTRODUCTORY MECHANICS PART II 


RIGID BODY MOTION 


What about around other axes? 


In = pj r dr dz [( x2 + y 2 + z 2 ) - x 2 ] t 
( x = r cos (j), y = r sin (|) ) 

© © 

In = P f dr d(j) dz [r 3 sin 2 (|) + rz 2 j , 


= pt J dr d(j) r 3 sin 2 (|) = pt ^ R 4 ^ “ s i n 2c()j | , 


rt 4 

= — ptR 
4 
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1 /2 

2jtp J* dr dz rz 2 = 2jtp - R 2 f dz z 2 , 



4 12 


Similarly for I 22 . In addition, one can see that the I.., i * j, vanish by symmetry. Example: 
I i 2 = -/ dv p (r) xy = 0 . 

n 

sin (J) cos (J) 


So, in this case we have 


L = I33CO303 + + ^22^2^2 • 


Since I n *I 33 , can see that L and 00 are pointed in different directions. Axes for which the I.., i * j, 
vanish are special and are called principal axes. We will find out how to identify them in a little bit. 

11.4 TRANSFORMATION PROPERTIES OF THE INERTIA TENSOR 

It is now clear why we are describing motion from the body axes. It is clear that the I take on 
different values for different orientations of our axes. This is why in general we do not try to describe 
the motion from the fixed frame. As the body moves, the I.. would become functions of time. 
(Sometimes it is convenient to describe motion in the fixed frame, however, if the I., are constant 
there. This is the case for a disk rolling down an inclined plane, say.) 


Can relate the kinetic energy of rotation to I . We had 


T rot 

1 

2 

2 m a 

a 

^ ^ot 

1 

2 

2 m a 

a 


(w X rj\ 
w 2 f a 2 - (to • rj 2 
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or 




in 


2 w i 2 x «* - 2 


= \ 2 2 m « 


a 


X ak - “i^jXaiXaj 


ai 


We have 


42 

COj^COj 

2 m « 

8 ij2 x «k - x . 

i,: 


a 

k 

=> T 

-L rot 

1 

2 

2 = \ 2 L : 

i/j i 

= -2 

I (0 

CD . 



ij i 

j 



In a rotated coordinate system: 

T = \ 2 I i: a) i co : • 


2 t: 


T is a scalar, CD. are vectors under rotations: 

i 

T = T. 

w k = ^ (actually, is oo k a pseudovector. 

i 

= (det X w a>’, ) 


co„ 


(11.14) 


so then 


T = \ 2 I ij ^ik^ j ^k“€ = T ' 


2 £ 

k ,1 
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In order for this equation to reproduce T ' = — ^ I k£ a) k co^, we must have that 

2 k ,i 

I k^ = , a tensor. 

i/j 

Can write this as a matrix equation also: 

4 = 2 ' 


or 


I' = X T I X {X, I are 3x3 matrices) 


Remember that V = X , so this relation may be used to get I in terms off: 

X (x T ix) X T = XI' X T , 
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or (matrix notation) 

I = XT X T (“similarity transformation”) 

or (index notation) 

= 2 = 2 * 

i/j ij 

11.5 PRINCIPAL AXES 

How do we find the axes for which I.., i * j, vanish, and do such axes always exist? 

Let us assume that I ±j is diagonal: 

. (no sum on i) 


As a matrix: 


I' = 


I, 0 0 


0 0 


0 

i' 


(11.15) 


(11.16) 


We have (switching the meaning of primed, unprimed quantities relative to (11.16)) 


M 


The left hand side, which is diagonal, may be written as where I ± are the principal values 

of the inertia tensor. Multiply both sides by X. and sum on i: 


r 


i k ,1 \ i 

=> tjV = 2 V* ‘ (no j sum) 


Can write fhs as 


- 2 Mi 6 - • 
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Can now group these terms together as 


(11.17) 


Free indices: j,m. This represents a set of 9 equations in 9 unknowns. That is, given a set of 
around some given axes, we have derived a set of equations that set conditions on the necessary 
to find principal axes. For these to have a nontrivial solution, we must have the determinant of 
the coefficients of the ^ vanish, for each j value (which is a free index). In matrix notation (the 
“secular equation”): 


det 


-1') 

*12 

*13 

I 2 1 

(I 2 2 - I' ) 

I 2 3 

*31 

*32 

(I 3 3 " I' ) 


= 0, 


or 


det (I - II') = 0, 


t 


unit matrix: 


1 1 0 0 X 
0 10 
0 0 1 


This gives a cubic equation in I'. The 3 roots of this equation then give the values of the 3 
principal moments of inertia: 


f = 


0 0 


0 0 


0 

i’ 


This is a completely general procedure and can be done for any given origin. This procedure gives 
the moments, but what axes are being referred to? 

There is another way of viewing this procedure. General experession for L : 

L i = 2 i i j « j / 

j 
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L = loo. (matrix language) 

What we are requiring above is equivalent to solving the equation 

Ico = Too. (oo: column matrix) (11.18) 

t t 

3x3 matrix number (“eigenvalue”) 

(either I ir I 2 , I 3 from above) 

Called eigenvalue problem, cos which satisfly the above are called eigenvectors. General form: 

number (“eigenvalue”) 

1 

A x = A f x. 

t t t 

square matrix column matrix ("eigenvector") 
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Values of moments of inertia about principal axes correspond to eigenvalues. Directions of principal 
axes correspond to eigenvectors. (One sees the same mathematics, called linear algebra, in many 
other places in physics, including quantum mechanics.) 

Written in matrix language, the above is 


(I - II’)© = 0 , 


or in explicit component language (also follows from Eq.(l 1.17) by multiplying both sides by crq , 
summing on j, and using co^ = ^ C0j ) , 

i 

2 ( z « - I,8 «) - 0 • (11.19) 

j 

Explicitly, 3 equations: 


{fll ^ ) ^1 ^ 12^2 ^ 13^3 ® f 

I21OO1 + (l 22 - I') co 2 + I 23 co 3 = 0 , > 

I 31 co 1 + I 32 oo 2 + (l 33 - I') co 3 = 0 . 


(1E20) 


Again, the condition that there be a solution says det(I - El) = 0, as before. Get something new 
this way: the priciple axis eigenvectors. These are the linear combinations of ce^, ce> 2 , ce> 3 which are 
associated with each eigenvalue, I lf 2 , 3 . How do we find the eigenvectors? After we have determined 
the values of the eigenvalues, l lf 2f 3 , we substitute one of the values back into the 3 above equations 
and solve for the ratio of values of cop cop co 3 . Why do we only get the ratio of values of co i 23 ? If 
wi solves (i = 1,2,3; these are eigenvector labels, not vector components) 


Ico 1 = Ipi) 1 , 


then co 1 = coo 1 , c arbitrary, also solves it. Thus, the normalization of the co 1 column vector 
is undetermined. This makes physical sense since we would not expect the magnitude of co i 
(corresponding to the rate of rotation) to be determined, but simply its direction (given by the 
ratio of components). 

Notice that if co i is given by (vector notation) 

co 1 = co 1 e 1 + co 2 e 2 + co 3 e 3 , (11.21) 
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we can build a unit vector out of it by dividing by co 1 


e , = 


co 1 


co 1 


CO, ^ CD. ^ CD. 

-r e i + e 2 + -r e 3 r 

(0 (0 CO 


( 11 . 22 ) 


J 1 = -y/t*? 


GT = a/00 ^ 2 + 002 2 + OO 3 2 


(11.23) 


We then recognize these components as the direction cosines studied in Chapter 1: 



cos y = 

0)3 

co 1 

cos (3 = 

C 02 

co 1 

cos a = 

0 i >3 


CO 1 


cos 2 a + cos 2 |3 + cos 2 y 


1 . 


So, if we know the 00 ^ : 00 ^ s 0 O 3 ratio for some eigenvalue I , we can always find the direction 

of the corresponding principal axis. This represents a unique axis of rotation since the choice of 
origin fixes a point through which this vector points. 
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There is another important property of the eigenvectors, go 1 . Consider a situation where one has 
distinct (unequal) eigenvalues, I 1 ^ I 2 . Then (matrix notation) 

loo 1 = I-lCjo 1 , loo 2 = I 2 oo 2 , 

=> oo 2 loo 1 = I^u^co 1 r cdTcd 2 = I 2 cdV . 

t T 

row matrix column matrix 


But 


iij = iji 
i 


cotco 1 = ^ C0 2 I i:j C0j = ^ COjl^CO 2 = (JO'IO ) 2 , 


cow = ^ co^co 2 = coW , 
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SO 


H^OTCD 2 = I^CD 1 , 

=> - I 2 ) oLco 2 = 0 , 

I x ^ I 2 => oTco 2 = (b 1 • o5 2 = 0 . 

Thus, eigenvectors of distinct eigenvalues are orthogonal, (point in T directions). If a double root of 
the secular equation occurs so that I x = I 2 , then nothing can be said about the relative directions 
of co 1 and CO 2 , other than that both are T to d} 3 . However, there is no loss in generality if we 
choose oo 1 * o5 2 = 0. One can prove that the principal axes can always be chosen to constitute 
an orthogonal set. Thus, for example we may always choose the eigenvalue orthogonalization/ 
normalization condition, co 1 * o5 j = 5 ij . 

! 

It is conceivable that the 3 could be complex since they correspond to the solution of a cubic 

! 

equation. A slight variant of the above argument shows that the I x 2 3 are in fact always real. 

11.6 PARALLEL AXIS THEOREM 

Lets locate the CM of a body and calculate the inertia tensor about it and another set of axes 
parallel to the CM ones. 

1 , 2 , 3 axes : J ± j 
1,2,3 axes : I ± j 
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On one hand, we have 

( L )i = 2 (Oi = 2 m « x (® x f «))i' 

a a 

for the 1,2,3 axes. We already showed this may be written as 

( E )i -1 - 

j 

=* 2 m « ( f « x (“ x f «))i = 2 ' 

a j 

where J\ . is the inertia tensor for the 1,2,3 axes. Similarly, for the 1,2,3 axes we have 
2 m « ( f « X ( W X ^)) ± = 2 Z i3°V 


which defines the CM inertia tensor I... Question: How are I.. and J.. related to one another? The 

1 j ^ ij j ij 

connection is given once one recognizes that 


r a = a + r a . 


Then 


© 

2 m « ( r « X ( W X r «)) = 2 ^ X (« X ? „ ) + 

a a 

© © 

2 m „ a x (ioxa) + ^m a (a x (w x !j + !„ x (to x a)) 


Lots of simplifications happen: 


(i) = ^ IijCOj , (I±j: CM inertia tensor) 

j 

(D = M a x (co x a) = M^ooa 2 - a (co • a)] , 
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or (in index notation) 

( ©)i = M 2 [6 ±j a 2 - a^jcOj . 

j 


Also 



0 . 


Therefore, we have 


2 J u<°i = M 2 [ 6 iT - + 2 • 


(11.24) 


> Apply now 
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The coefficient of coj on both sides of this equation must be the same: 

inertia tensor 
for a point 
mass at CM 

,- A -\ 

=> J ±j = Iij + M |^5 ij a 2 - a i a j j “Parallel Axis Theorem” (11.25) 


As an example, consider again a cylinder of radius R, but rotated from a point on its edge rather 
than in the center as before. 



rotation 


1 

We know from before that I 33 = “ MR 2 . Since the vector a., pointing from the 3 axis to the old 

( 3) axis has no components along this direction, we have J = — MR 2 . 

2 

11.7 EULER ANGLES 

We have written K.E. and angular momentum as 


T 


\ 2 

^ i. i 


L, 



under the conditions that the body is rigid and that the origin of the body system is either at the 
center of mass or at a point for which R = 0 . The analogous results for linear motion are of course 



P, 


Mf ± . 


There is a major difference in handling the two cases: When the r ± are integrated, we get linear 
coordinates that tell us, at each instant in time, the location of the particle or body. The coi, 
however, when integrated are not quantities which give us the orientation of a rigid body in space. 
In order to describe the orientation, 3 angular quantities will be needed just as 3 linear quantities 
are needed to locate a body’s origin. Many different choices or schemes are possible. A conventional 
and convenient choice are the “Euler angles”. 
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First, I will simply describe the situation. (Better pictures are available in other books.) Final situation: 



0, \|/, (f> are the Euler angles. (As shown, all are positive.) They represent a choice of generalized 
coordinates useful in defining the Lagrangian for a moving rigid body. 

General relation relating rotated coordinates: 

x = A,x\ (matrix notation) (11.26) 


Represents a passive rotation. First, rotate positively (counterclockwise as viewed from above the 
l’,2’ plane) through an angle 0 about x 3 : 


x" = Lx', 


where 


K = 


/ , x „ \ 

cq) sq) 0 

—s(|) c(|) 0 

0 0 1 


Special notation: 

c(J) s cos 4> , s(|) 3 sin (|) 


(11.27) 


(11.28) 


(11.29) 
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! 

Now rotate (counterclockwise) by 0 about X] 


x"' = X 9 x" , 


A.0 = 


1 0 0 x 
0 c0 s0 
0 —s0 c0 


i 

Next, rotate (counterclockwise) by op about 

x = X^x'" , 


\ = 


cap sap 0 
-sap cap 0 
0 0 1 


Put it all together: 

* = \ (X e (X*x')) = (\,X e X $ ) x'. 


(11.30) 


(11.31) 


(11.32) 


(11.33) 


(11.34) 
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Call A, s . Explicitly, 




' (dpC(|) — C0S([)STp) 
(-SXpC(|) - cdstpcxp) 
v (s0S(|)) 


(cXpS([) + C0C([)SXp) 
(-SXpS(|) + C0C(|)CXp) 
(-S0C(|)) 


(sxps0) N 

(cxps0) 


We may now think of the angular velocities tp, 0, cj) as describing the instantaneous state of rotation 
of the body. [Directions come from infinitsmal rotations of V , 0 , o r (|) and the right hand 
rule.] We would like to find the relationship between these quantities and the COj 2 3 . (to represents 
of course the instantaneous angular velocity projected on the body axes.) 


First of all, it’s obvious that 


xp = xpe 


3 * 


(11.35) 


What about 0 and 4* ? It’s also clear from the picture that (0 is in the 1,2 plane) 


0 = 0 cos xp e x - 0 sin xp e 2 


(11.36) 


Let’s now find the components of 0 using the Euler matrix. All we have to do is: 


x = A.x' where x' = 


0 X 


(b , 
y / 


It is easy to verify that 


x = 


/ • \ 

(|) sin ip sin 0 

§ cos xp sin 0 
(j) cos 0 


(j) = (j) sin ip sin 0 e x + (j> cos ip sin 0 e 2 + (j) cos 0 e 3 . 


(11.37) 


Download free eBooks at bookboon.com 


158 








MODERN INTRODUCTORY MECHANICS PART II 


RIGID BODY MOTION 


Adding it up gives 

d) = + (j) + 0 , 

=> o5 = ((j) sin 4> sin 0 + 0 cos ip, (j) cos ip sin 0 
-0 sin xj3 f ijj + (j) cos 0 ) 

or 

r 

(o 1 = ^) sin 0 sin tp + 0 cos tp, 

■ to 2 = (j) sin 0 cos tp — 0 sin 
(0 3 = Ip + (j) cos 0. 


(11.38) 


11.9 EULER'S EQUATIONS OF MOTION 

We are going to study the motions of rigid bodies using 2 methods. One, we already know, is the 
Lagrangian method. This is the reason Euler angles have been introduced. The Lagrangian approach 
is clear. The kinetic energy is 



We now express 00 123 in terms of the Euler angles 0, \|/, (|) and their derivatives and use the 
Euler-Lagrange equations of motion, 

<3L _ d ( dL \ _ 

dqi dt \3q ± J 

where q± G (0? TT ( t ) ). 

There is another approach possible which works directly with the cd 123 . Remember: 

dQ\ 

\dt /fixed 



+ to x Q . 


rotating 
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Take our fixed and rotating axes origins at the CM of the particle and instantaneously coincident 
in direction. Let Q —* L , 


dh\ 

dtjf 



+ co x L . 


But 


\dt/ f 


N and (£,). = ^ ' 

j 


and also 





)■ 


But since we recognize that the body axes are attached to the body 


dt 


I ij co j = I ±j cb. f 



Download free eBooks at bookboon.com 



Click on the ad to read more 









MODERN INTRODUCTORY MECHANICS PART II 


RIGID BODY MOTION 


then 


N i = 2 I i: tb : + 2 * (11.39) 

j j*M 

If we now pick our body axes along the principal axes, then I ± j (i ^ j) = 0 , and 
2 8 ijk “j®/ = 2 E il* I k • 

=ia, j ' k 


Written out explicitly, we have (“Euler’s equations”) 

N 1 = I 1 (b 1 + (0 2 (0 3 (l 3 - I 2 ) / 

N 2 = I 2 cb 2 + (0,(0, (l 3 - I 3 ) , • 

N 3 = 1 3 (b 3 + (0,0), (l 2 - I 3 ) . 


(11.40) 


They are completely equivalent to the Lagrangian equations but in general have a different form 
when co i 2 3 are expressed in terms of 0 f \\f, (f). They are coupled first order differential equations. 
One strategy: solve for co i 2 3 (t) and then invert the above relationships between 03 12 3 and 0 T V|/, (|) 
and their derivatives to solve for 0(t), \]/(t), (|)(t). 

11.9 SYMMETRICAL TOP - EULER SOLUTION 

Lets try carrying out this procedure for torque-free motion (n = 0 ) for a body with = -k * ^3 

\ v J 

(a “symmetrical top’ *). 1 12 
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With no loss in generality, take cd 3 > 0. We can realize such a situation in an orbiting space craft 
or a falling elevator, say. 


Of course, if CD is initially pointed along a principal axis (along 3 or anywhere in the 1,2 plane), 
then the motion is trivial. Then we simply have, 


0 = I 3 cd 3 or 
=h> (D 3 = constant 


0 


I 12 co 


1/2 


or CDi 2 


constant. 


In the general case, we have to solve: 

Xi 2 cbi + cd 2 cd 3 (I 3 — I 12 ) = 0 t 
1 12 cb 2 + o^co^I^ - I 3 ) = 0 , > 

I3CO3 = 0 . 


Notice immediately that <jG > 3 = const, in the general case also. We may now write the other two 
equations as 

(h 1 + Qcd 2 = 0 , 

co 2 - Qu> 1 = 0 , 

where (the “precession rate”) 


Q ss 



t 

(either positive 
or negative) 


t 

positive constant 


(11.41) 


Then we have 


CD, 


+ cd 2 £2 = 


CD, 


+ cd -^£2 = 0 • 


(11.42) 


Solution: 


CD 1 (t) = A cos (Qt + a). 
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Since 


1 . 

co 9 = - — ( 0 , r 
2 Q 1 

=> oo 2 = A sin (Qt + a). 

(A, a = constants; can take a =0 with no loss in generality.) We recognize these as the parametric 
representation of a circle since A 2 = 0) 2 + CD 2 . Describes (note that |co| =constant.): 



prolate 



oblate 
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L, co, e 3 are all in the same plane: 

L = 1 , 2 ( 0 )^ + co 2 e 2 ) + I 3 e 3 co 3 , 

co = (o 1 e 1 + co 2 e 2 + co 3 e 3 , 

so 

L • (w x e 3 ) = L • (-oc^e., + 0 ) 262 ) = 0 

Somewhat pitiful pictures: 



or 



to 


Case 1 : prolate 


Case 2: oblate 


L is a constant in magnitude and direction in the above diagrams. We have 


tan (3 = 


L* a 


l 3 l 3 


/ _ 

co - co 3 e 3 
Jco - co 3 e 3 | ; 


I12CO1 + ^-12^2 

I3CO3 




2 2 

+ co 2 


tan p = hi 'fit + ^ 

I, 


■.3 CO3 


Clearly P is a constant in time. Also 


tan a = — = 

OCL 


+ C 0 o 


ocl 
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is also a constant. They are related by 


tan |3 _ I 12 
tan a I 3 

Therefore: 


If 

If 


13 < I 12 (Q < 0) => (3 > a , 
1 3 > I 12 (Q > 0) => a > (3 , 


(Case 1: prolate) 


(Case 2: oblate) 


(11.45) 


Motion is summed up in the following pictures: 



L w , x 3 in 
same plane 



Case 2 : Q > 0, a > (3 
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Notice: d) rotates clockwise w.r.t. 1,2 (body) axes in Case 1, whereas it rotates counterclockwise 
w. r . t . 1,2 axes in Case 2. Also notice if I 12 = I 3 , L and (So point in the same direction and 
there is no precession. (The case for origin at CM of a sphere or cube.) 


Now complete the discussion by getting 0(t), V|/(t), (|)(t) from co i23 (t). Remember: 


oo 1 
< oo 2 
o ) 3 


(j) sin 0 sin ip + 
(j) sin 0 cos ip - 
1 R H- (j) cos 0. 


0 cos ip, 
0 sin ip, 


or, inverting, 

^ co 1 sin ip + oo 2 cos ip 

— r 

sin xp , 

sin xp + co 2 cos xp). 


sin 0 

0 = co 1 cos xp - oo 2 
xp = oo 3 - cot 0 (co 1 


(11.46) 

(11.47) 

(11.48) 
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If we take x 3 along L, then we identify 
0 = p.(b = const.) 

Therefore 0 = 0 and (from (11.47)) 

0 = co 1 cos if - co 2 sin if , 

=^> 0 = A cos (Qt) cos if - A sin (Qt) sin if , 


or 


0 = A cos (Qt + if). 


If A * 0, then we may choose (\|/(t) gives rotation angle w. r . t. line of nodes) 
if(t) = -Qt + y, 

without loss in generality. And, from (11.46), 

(]) = -(cos(Qt) sin if + sin(Qt) cos if) , 

sin 0 


d) = ——— sin (Qt + if (t)) = ——— , 
sin 6 i_. V 11 sin 6 


(choice) 

4>(t) = , A t + c (e = p) 

sin p 


Eq.( 11.48) is now satisfied also since 


-Q = ( 0 , - cot 0 


sin 0 


sin 0 L 


cot 0 = cot P 


I3CO3 


Il2^ + ^2 


(11.41) 


'A-i 

^12 / 


I3CO3 


ii2^r + ^2 


A = l + CO2 . V 


Download free eBooks at bookboon.com 


(11.49) 


(11.50) 


(11.51) 


167 














MODERN INTRODUCTORY MECHANICS PART II 


RIGID BODY MOTION 


Thus, the time dependence of the Euler angles is given by 


/ v JT / \ + 

0 = (3, Xp(t) = -Qt + — , (j)(t) = -;--- 

2 sm p 


t. 


(|) (t) may also be written as: 




sin a 
sin(a - (3) 


£2t, 


where a, |3 are defined above. Thus, Q is the precession rate of the line of nodes. 

sin (a - (3) 

11.10 SYMMETRICAL TOP - LAG RAN GIAN SOLUTION 

We have used Eulers equations to first get co 123 (t), then use these to get 0(t), \|/(t), (f>(t). Other 
way of proceeding: Lagragian using 0, \| j, c|) from the start. Now do the same problem (force 
free symmetrical top) from the Lagrangian point of view. We have 

T = ^ • 

Z if j 


Again, pick principal axes with symmetry axis of body along x : 


T = 


(<x > 2 + COj) 


I 3 CO 3 . 


But now use earlier expressions for C0j 23 in terms of 0, (j), ij’: 

I 12 (f) 2 + (j) 2 sin 2 0) + I 3 (ip + <j) cos 0) 


1 

T = - 
2 1 


The Lagrangian equations are (L = T here): 


dT 

- — 1 

dT 

d0 

dt 

UeJ 

dT 

d 

'd T \ 

deb 

dt 

v defi / 

dT 

d 

' dT 

dtp 

dt 

A 


( 11 . 52 ) 


(11.53) 

(11.54) 

(11.55) 
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We have 


dT 

ae 

dT 

d<j> 

dT 

dip 

r 

dT 

d0 

dT 

d(J) 


li 2 0, 

I 12 (j) sin 2 0 + I 3 + <i> COS 0 ) cos 0 , 

I 3 + (j) COS 0 ) , 

I 12 (j) 2 sin 0 cos 0 - I 3 (j) (ty + 4> cos 0 ) sin 0 

® - 0 . 

dtjt 


Putting the pieces together gives us the Lagrangian equations of motion: 

(11.53) =* i 12 0 - I 12 <j) 2 sin 0 cos 0 + I 3 <f> (ip + 4> cos 0 ) sin 0=0, 


(11.56) 
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(11.54) => [l 12 4> sin 2 0 + I 3 (ijj + (f) cos 0 ) cos 0j = 0, 

(11.55) => + <j> cos 0 )] = 0. 

Euler’s equations can be recovered from these. For example, substitute in (11.58): 
= co 3 - cot 0 (o) 1 sin \\) + co 2 cos t|t). 


4> = —-— (to 1 sin i|t + oj 2 cos t|)), 
sin 0 


\\) + (() cos 0 = co 3 - cot 0 (003 sin \\) + to 2 cos 

-0)3 sin tjj - co 2 cos t|t) 


(11.57) 

(11.58) 


Thus, Eq.(11.58) just says 

I3CO3 = 0 , 

which is one of the Euler equations above. Likewise, (11.56) and (11.57) above are equivalent to 
the other two Euler equations (after some manipulations). 

Can also apply the Lagrangian approach in describing the symmetrical top with one point fixed: 
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Same T as above, except now L = T — U, where U = mgh cos 0. From the previous results, one 
can show that the total energy, E = T + U, can be written as 


E = ^I 3 «3 + ^Il2 0 2 + U eff (0), 


where the effective potential, U eff (0), is given by, 

(P* - cos 0) 2 _ _ 

E e ff(0) _ ■ 2 n + COS 0' 

2I 12 sin 0 

where (11.57), (11.58) give rise to the conserved conjugate momentums, 

px s = I 12 4) sin 2 0 + I 3 (tjj + <j) cos 0 ) cos 0, 
v 3(j) 

dXj 

s = l3 ^ + ^ cos e ) ' 


(11.59) 


(11.60) 


(11.61) 

(11.62) 


The minimum of this effective potential gives the angle of inclination at which the top precesses 
steadily. 
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11.11 PROBLEMS 

1 . 



A thin solid metal bar (p = const.) of mass m and length 1 is rotated uniformly about one 
end at an angular velocity, w, directed along the 3 axis. 

a) Compute the kinetic energy. 

b) The origin in (a) is given a velocity, v, along the 2 axis. Now compute the kinetic 
energy. (Be careful!) 

2. (a) Given that the potential energy of a small mass m a is given by v a = gm a x a , where 
x a is the height of m a , show that the total potential energy of an arbitrary object in the 
gravitational field is given by V = gMX, where X is x-coordinate center of mass coordinate, 
X = (R C m)x> and M is the total mass. 

(b) Show that the total torque on the object, N, is given by 

N = R cm x F , 

where R CM locates the center of mass and F is the total force. 

3. (a) A thin uniform rod of length b and mass m stands vertically upright on a rough floor 
and then tips over. What is the rods angular velocity just before it hits the floor? 

(b) Repeat for a rod initially standing upright on a smooth floor (no friction). Find the 
final angular velocity about the CM just before striking the floor. 

(c) You should have gotten the same angular velocity for both (a) and (b) parts. Why? 
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4. A solid disk (p = const., radius = R, mass = M) rolls without slipping on the ground at a 
constant rate, 00 = 0 = const. Compute its kinetic energy two ways. 



a) By using Eq.(11.5) of the text. 

b) By using the fact that the point of contact between the disk and ground acts as the 
instantaneous axis of rotation. (You will have to use the parallel axis theorem to do this part.) 

Do you get the same answers for (a) & (b)? 

5. A pendulum consists of two point masses, ny and m 2 , attached to a completely rigid, 
massless bar of length C . m is attached halfway down the bar, m 2 is at the end. 



a) Find the bar s moment of inertia about the pivot point. 

b) Assuming small oscillations, find the bars period of motion. 
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6. A pendulum consists of a thin bar of mass m and length C attached to a pivot point as 
shown. It swings in a plane. 


/ / / / / / 



a) Find the bar’s moment of inertia about the pivot point. 

b) Assuming small oscillations, find the bars period of motion. Does it swing faster or 
slower than a pendulum of the same length with mass m concentrated at its end? 


7. (a) A “physical pendulum” consists of an arbitrary object of mass M whose moment of 

inertia about the pivot point is I. (Motion still confined to a plane.) Given that R is the 
distance from the center of mass of the object to the pivot point, show that this pendulum 

has the same small oscillation period as a simple pendulum of length L = -. 

(b) Find the length of simple pendulums with the same period as problems 11.5 and 11.6. 


8. The moment of inertia tensor, I, for some object has the form, 


I = 


1 X M 0 ^ 

M X 0 

0 0 Y 


where X, Y and M are some given real numbers (X, Y > 0 and |M| < X), find: 

a) the three moments of inertia, I ± , i= 1,2,3. 

b) the angular rotation eigenvectors, co 1 , associated with the principle axes. 
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9. (“Slamming door problem”) A door is constructed of a thin homogeneous slab of material; 
it has a width w and a height t . You are given that the line of hinges bends inward 2° 
with respect to the vertical, in the 3-2’ plane; see the figure. Assuming frictionless hinges, 
what is the angular velocity of the door just before it closes if it starts from rest after being 
opened 90°? [Hint: Find the Euler angles which transform from the body to the fixed axes 
and conserve energy.] 
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10. (a) Given principal axes for which the angular momentum vector may be written 

L = I 1 to 1 e 1 + I 2 co 2 e 2 + I 3 co 3 e 3 , 


where e 1>2j3 are body frame unit vectors, show that the components of L along the 
0 , $ and ip directions are 

( 1 1 = 1 2 s 1 12 ^ I 3 case) 

L 0 = I 12 0 , 

= I 12 (j) sin 2 0 + i 3 (tp + (j>cos)cos0, 

\ = I 3 ( + (j)COS) 0 

(b) Alternatively, given components of co along the 0 , (j) and directions as <a e , co t , and co v 

respectively, show that one may write 

L 0 = I 12 co 0 , 

L * = 1 12 (co^-co^cosO ) + I 3 co^cosO, 

\ = I 3 
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11. Eqs.(11.56) and (11.57) are supposedly equivalent to the Euler equations: 


1 12 

+ (I 3 - I 12 ) © 2 © 3 = 

0 , 

(1 

112 ^2 

+ (I 12 - I 3 ) ©3©! = 

0 . 

(2 


By using the expression for <jc > 123 in terms of Euler angles, show that Eqs.(11.56) and 
(11.57) may be recovered from the above. [Hint: Consider (1) cos \|/ - (2) sin \|/ and (1) 
sin \|/ + (2) cos \|/.] 

12. Show the expression for the effective potential for the symmetrical top, Eq.(11.60), 

(P* - P„, cos 0) 2 ___ 

u e ff (0) - ~ : 2 + cos 

2I 12 sin t) 

follows from the expression for the total energy (T in Eq.(l 1.52) and U = mgh cosG) and 
the definitions of the cyclic momentums, P§ and P\p / given in (11.61), (11.62). 

13. A solid sphere (p = const.) is made into a fixed point top by supporting it on a point as 
shown. Gravity is acting downward. 



Find the minimum angular frequency, 00 min , such that the sphere is stable in a vertical 
position. (Use the result of #12 above.) 


14. Consider a thin rod for which I 3 = 0, making the Euler angle \|/ irrelevant. Show that the 
effective potential for this system has a minimum at an angle 0 Q given by 


cos 0 O 


mgh 

I > 2 
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Other Problems 

15. Given that the l’,2’,3’ axes are principal axes with moments of inertia I 115 I 22 , I 33 
respectively, find the moment of inertia along the new 3 direction, I , in the Euler angle 
diagram, p.11.20 of the text (angles are 0, \|/ and (f>), using the transformation law of the 
I.. under rotations. 

16. Show that the ratio of the magnitude of the instantaneous angular velocity, |co|, to the 
precession rate, Q, for a symmetrical top undergoing torque-free motion, is given by 

|co| _ sin (3 
Q sin (a - (3) ' 

where a is the angle of the vector co with respect to the 3 (body) axis and P is the angle 
of the instantaneous angular momentum, L ? with respect to the same axis. 

17. A “spherical pendulum” is a simple pendulum that is unconstrained in its angular motion; 
it can move freely about it’s pivot point. Consider a spherical pendulum of length L acted 
on by gravity with a mass M attached to it’s end. The mass of the attachment of length L 
is negligible. 
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a) Derive the Euler-Lagrange equations for the pendulum in the spherical coordinates 0 
(polar angle) and (|) (azimuthal angle). Take 0=0 to represent the equilibrium position. 
Show that the z-component of the angular momentum, £ z , is conserved. 

b) For uniform circular motion in f when 0 =const. ^ 0, show that the period of the 
motion may be written as 

\ i/2 

L cos 0 | 

~r~) • 

c) Consider general motions of the pendulum subject only to the condition that 0 « 7t/2. 
Defining the total energy of the system, E, relative to 0=0, show that the turning points 
of the 0 motion are 



0 t 



18. A symmetric top (J l = I 2 = I 12 ) of mass m and constants of the motion and precesses 

steadily, (j) = const. Let = p^ = p. h is the distance from the CM to the fixed tip. 

a) Find the angle, 0 Q , at which the precession occurs. 

b) Show the precession rate is 

• I mgh 

19. (a) A thin uniform rod of length b and mass m stands almost vertically upright on a rough 
floor and then tips over without its end slipping on the floor. The initial tipping angle is 
0 o« 1 with respect to the vertical. Find the amount of time, T, that the rod takes to fall 
over and hit the floor. 

(b) The same problem as in (a) but with a smooth, frictionless floor. 

20. (a) Given a set of principal axes and all axes parallel to them for some object, show that 
the sum of the moments of inertia is minimized when the origin is at the center of mass. 

(b) An object has principal axes for which I n > I 22 > I 33 . Characterize all the points which 
are also principal axes for axes parallel to the original set. 
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Extra: Using the concept of eigenvalue degeneracy (for moments of inertia) and the parallel 
axis theorem, find a larger set of axes which are also principal axes for the object in (b), 
but which are not parallel to the original set. 

21. Consider the symmetric top (In = 1 22 * I33) with one point fixed. We briefly 
discussed how this system can be solved with a Lagrangian approach. However, now attempt 
to write Eulers equations for this system. Can it be done? Remember the components of 
the equations refer to the body system. If these equations can be formulated, discuss the 
manner in which they may be solved (but do not attempt it!). If not, tell me why not. 

22. Consider a solid cylinder of mass M, length t, and radius R of uniform density, p. Let the 
1, 2 plane be located at the top of the cylinder. (This is a little different from the situation 
in the text on p. 11.5.) 



a) Find the elements of the inertia tensor with respect to the unprimed axes. 

b) Find the elements of the inertia tensor with respect to the primed axes. 
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12 COUPLED OSCILLATIONS 


12.1 COUPLED DYNAMICAL EQUATIONS 

In a number of cases we have derived coupled linear differential equations. We have treated such 
systems as special cases, but it is now time to develop some general techniques for handling them. 
Mathematically, we will encounter here both the use of complex numbers (which were used 
in the Focault pendulum discussion and various differential equation solutions) as well as the 
eigenvalue-eigenvector matrix algebra of the last chapter. Examples of systems these techniques will 
cover: conservative, linear mechanical or electrical oscillations, molecular vibrations, approximate 
planetary motions, and many more. We will also use stability analysis and the concept of generalized 
coordinates. I will first introduce the theory and then will work out a number of examples to 
illustrate the general techniques. 


Let us consider a system with generalized coordinates <3j (j = 1,..., n; n = no. of “degrees of 
freedom”), related formally to the x ai by 



- o. 

at 


( 12 . 1 ) 
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Usually, the q. we will consider will be measured from equilibrium positions or lengths. For example, 

/////// 



Let us also only consider U = U yq^ j . This means we are limiting our discussion to conservative 
systems for which H = T + U = const . in time. Another limitation of the following discussion 
will be the assumption that the system is linear. That is, we will only consider potentials that are 
quadratic in the generalized coordinates, or are approximately quadratic. For small oscillations, 
we have 



-<jJ - + 2 


du 


T* a <3k 


0 


I 2 

Z j,k 


d 2 u 


d qj dq k 

:k k: 


o + 


( 12 . 2 ) 


(constants) 

where higher order terms will be neglected. In addition, we will assume or require that (q , q ,... 
etc. are equilibrium q.s): 


1. U 0 = 


2 . 


du 

5q k 


3. 


a 2 u | 

dQ 2 1 ° 


U(qi = qiof T2 = q20f • • • ) = 0 . 

= 0 for each k. 

> o, where Q is any linear combination of the generalized coordinates, q.. 
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Condition 1 is no restriction in the applicability of the analysis since the absolute value of the 
potential U is always arbitrary up to an overall constant. Conditions 2 & 3 above are just conditions 
for stable equilibrium in a system with n degrees of freedom. Notice that condition 3 above (which 
is a sufficient, rather than necessary, condition for stability) is equivalent to the requirement that 
U(q.) > 0 for quadratic potentials, given U Q = 0. We can show this as follows. Assume (i = 1,..., n) 

Q = 2 x i q i ' (12.3) 

i 


where the x. are arbitrary constants. Then 

d U d ^ dq ± 

<3Q 2 dQ 4* dq i dQ 


i dq i 1 
but 1 = — so 

dQ x ± 

d U d ^ 3u(c[i) 1 ^ d U 1 

dQ 2 dQ 4* dq i x L dq^dq. x ± x. * 

Condition 3 now gives, 

V d 2 U I 1 

Tfi dq.dq^ 1 x iXj 


(12.4) 


(12.5) 


(The sum in (12.3) is only over the nonzero x values.) Since the x are arbitrary (but nonzero) 
constants, this is the same as 


U(q±) > 0, 


( 12 . 6 ) 


where = — (see Eq.(12.2)). 
x i 

Because of our condition —— = 0 , we know from Chapter 7 for a system of point particles, 
for example, that we may write (the m here were called a. k there) 


T = o 2 m jk f 

2 j,k 


(12.7) 


where 


m 


jk 


= 2 


m r 


3x ai 3x ai 


a,i 


aqj dq k 


= m 


kj 


( 12 . 8 ) 
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Note that in general m. k = m. k (q.) (functions of the q.), unlike the A. k above, which are just constants. 
However, consistent with the restriction to quadratic terms in U, we will restrict ourselves to 
quadratic terms in the 3 j. This means that we will expand 

ra jk (qj = ra jk | 0 + 2 ?r lo + • • •' (12.9) 

£ 0C T_ 

neglect 

and neglect everything except the first term, which is assumed non-zero. The reason is that terms 
like - q^q.q k i n the Euler-Lagrange equations produce nonlinear terms in the differential 
equations. (Neglecting such nonzero terms give additional restrictions on what is meant by “small 
oscillations”.) Notice also that, in general 

dT 1 v* 3m ik • • 

T- - rX -rf g i q *' 

dq f 2 " dq e 

so that with (12.9) we have 

3T 

- = 0 . 
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After these restrictions, the general form of our Lagrangian in generalized coordinates, q., is 


real constants 

1 1 

1 „ • • 1 - 

L=T- U = —2/ m jk qjqk - - z A jk qjqk • 

^ j,k ^ j,k 

Equations of motion are given by 

<3L d dL 

-— = 0. 

dt dq k 


— 2) ^jk Tj * (Ajk A k j) 

j 

= J m jk • ( m jk = m kj ) 
j 

^ ^ (Ajk Tj + m jk Tj) = ® • 
j 


<3q k 

We have 

dU 

5q k 

dT 

<3q k 


( 12 . 10 ) 

( 12 . 11 ) 

( 12 . 12 ) 

(12.13) 

(12.14) 


Although we have derived these equations in the context of particle oscillations, the small oscillations 
of many realistic rigid body systems can be so characterized as long as the Lagrangian may be put 
into the form of (12.10) above. 


12.2 EIGENVALUE/EIGENVECTOR SOLUTION 

We will use complex number analysis to simplify the solution of this sytem of equations. Under 
conditions 1, 2, 3 above, we know that the solutions are oscillations. Therefore assume 


qj (t) = e ia)t (12.15) 

T T 

complex (involves real frequency 

two undetermined 
constants) 


where the real part of the right hand side is understood. Plugging this back above, we find 

2 ( A :k - w 2 m jk ) aj = 0. (12.16) 

j 

[This is almost the same form as the equations which determine principal axes: 

- I’Sjm, - 0.] 

j 
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Again, the condition that there is a nontrivial solution for the a. is that the determinant of the 

matrix (A.. - CD 2 m.) be zero: 

1 j i r 


det 


A n 

- co 2 m n 

^12 

- co 2 m 12 

■^■21 

- (o 2 m 21 

^22 

- (o 2 m 22 


• • 

• • 


0 . 


(12.17) 




/ 


this gives an n— order equation in to 2 , the roots ofwhich will be labeled at co r = 1,2,... n. 
These are called the eigenfrequencies or characteristic frequencies. In general, the above is more difficult 
than solving the analogous equation for the principal moments of inertial since the 00 2 terms (the 
analog of the priciple moments, I’) enter off-diagonal elements of the matrix also. However, one 
is solving the same mathematical problem in either case. 


Analogy: 


Inertia tensor 


Coupled Equations 


eigenvalues: 


Ii 

(i = 1/2,3) 
(principle axes) 


(r 1,2,...,n) 
(oscillatory modes) 


eigenvectors: 


CO 1 ( COj ) 

j = {1/2,3} 

(vector components) 


a r (a]) 

j = {1/2,...,n} 
(vector components) 


Just as there are directions in physical space which render the inertia tensor diagonal, there are 
directions in mode space, associated with the eigen frequencies to 2 for which the motions uncouple. 
We get the a r by the same procedure as in the last chapter: substitute a known in the algebraic 
equations, and then solve for the ratios 

r r r r 

a i : a 2 : a 3 : . . . : a n 
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for a given eigenvalue, 00 ^ . As before, the overall normalization of the a r (chosen real) are arbitrary. 
This makes physical sense here also since this corresponds to the amplitude of the motion, which is 
not determined by the equations of motion but by the initial condition. Since there are n a r s, we 
can require n condtions to arbitrarily normalize them. In the inertia tensor case, we required that 

cb 1 * cb 1 = 1 , i = 1, 2, 3 . 

Here, we require, 

^ m jk a j a k = 1 f ^ = If 2 f . . . n. (12.18) 

(That this combination is always positive can be shown from the positivity of T, the kinetic energy.) 
Also, we can prove that (just a generalization to an n-dimensional space of the similar proof in 
Ch. 11 for the inertia tensor) 

2 m jk a j a k = 0 ^ (r * s) (12.19) 

j fk 
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for co^ CDg , and that for degenerate roots (co^ = ) we may always choose this equation 

to hold. The above two conditions can be written together as 

2) m jk a j a k = ^rs * (12.20) 

j/k 

Also, we can show that the oo^ are all real. In fact one can show, under conditions 1,2,3 above, 
that co* > 0. 1 Without loss of generality, one may further choose co r > 0. 

The general solution is now given by (real part still understood) 

q,(t) - 2p c a'e‘"- 1 , 

r = l 


where the are now all determined (up to an overall sign) and |3 r are so called “scale factors” 
which, in general, are complex. (We still have not built in the initial conditions, so we still need two 
arbitwrary constants for a second order differential equation.) The [3 r are the amplitudes associated 
with the r th eigenmode as determined by the initial conditions. If only a single |3 r is nonzero 1 , 


|3 r = 0, r * k, [3 k * 0, 

then only a single eigenmode of the system has been excited and the solution of the motion is 
particularly simple. 

/j_\ o k icokt 

qj(t) = |3 k aj e K (no k sum) 


Or, introducing the “normal coordinates”, 

n k s P k e 1<0kt , (12.22) 

we have 

qj(t) = n k aj. (12.23) 

In general, 

qj(t) = 2 n k af . (12.24) 

k 
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Expressed in normal coordinates, one has (using the orthogonality condition, ( 12 . 20 ), as well as 
(12.16)) 


= \ 2 % ' u = \ 


= 0 , 

dn dt l an 1 


n r + to r n r = 0 . 


(12.25) 


This in turn shows that the energy associated with each normal mode is constant since for the r 
mode, 


th 


+ (o’n r 

V, 

II 

1 * 2 

1 2 

r n r + 

— co r n 


= T r + U r 


(12.26) 


Final cookbook recipe for solving a system for small oscillations: 

1. Write the Lagrangian, L, of the system in terms of generalized coordinates, qi, and find the 


2 . 

3. 


Ajk and mjk either by using the explicit formulas (Aj k = - 


a 2 L 


dqjdq k 




d 2 L 


dqjdq k 


or by comparing to the form of the Lagrangian, Eq.(12.10). 

Form the matrix ^Aj k - 0) 2 nij k j and find the n eigenvalues, CD 2 . 


Determine the ratios a* 
to Eq.(12.20). 


a 3 : . . . : a. r n and normalize the according 


4. Write the general solution for the q. as in Eq.( 12 . 21 ). Physical interpretation follows from 
an examination of the motion of the normal modes (setting each P =1 and other P’s= 0 ). 


5. 


Apply the initial conditions and find the P r . 
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If T and U for a system are not given, but the equations of motion are, then obviously we may 
skip step 1 and instead form the characteristic equation using the ansatz Eq.(12.15) for the q.. As a 
faster procedure, note that we may also skip the normalization/orthogonality condition, (12.20), and 
determine q^t) = ^ fcq e lc ° rt (P r absorbed in the new ) directly from the initial conditions. 

r 

Note that in step 3, since the eigenvector equations determine only the ratio of the , we need 
consider only n-1 of the equations in determining all such ratios. Also note that, similar to the 
case of degenerate eigenvectors for principal axes (see discussion on p.11.16 and prob.#11.8 in 
the simple case that M=0), there is an indeterminacy in the eigenvectors corresponding to the 
repeated roots of the characteristic equation. For example, in the case a single repeated root the 
indeterminacy may be removed by any arbitrary specification of one of the roots - the other will 
then be determined by Eqs.( 12.20). 
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12.3 EXAMPLE 

Enough of the theory, let s do some examples to make this more understandable. Coupled masses 
example (this is mathematically the same for small oscillations as the coupled pendulum problem 
at the start of this chapter with M=m =m , x—and k = r £ being the pendulum length; 
see also prob.3.9): 



m i = m 2 = M 

Let x , x 2 (playing the role of the qs) be measured from equilibrium positions. Step ( 1 ): 


T = ~ M (x 2 + x\) . 

=> m n = m 22 = M , m 12 = m 2 i = 0, 

U = - kX ! + - kX 2 + - k 12\ X 2 - X l) r 


a 2 u 

A n = -7 

11 dx? 


o = k + k 12 , A 2 


a 2 u 

dx? 


= k + k i 2 . 


A 12 


a 2 u 


5 x^X 3 


o - k 12 / A 21 - ” k 12 • 


Step (g): 


det 


k + k 12 - to M 


-k, 


-k, 


k + k 12 - 0) M/ 


= 0 , 


(k + k 12 - co 2 m) 2 - k 2 2 = 0 , 


to 1 = 


k + 2k, 


m y 


/ to 2 — 


M 


1/2 
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Step ( 3 ) : 

Equations for eigenvectors, r = 1 case (for n=2, only one of the following eigenvector equations 
is necessary): 

(An — (D 1 m 11 ja 1 + (a 21 — co 1 in 21 ja 2 
(a 12 - w?m 12 ) a i + ( A 22 - w i m 2 2 ) a ! 

. — ^ 12 a i 

— ki2 a i 

r = 2 case: 

k 2 1 2 z z 

12 a i — ki2 a 2 ^ a i = a 2 


k 12 a 2 = 0 

ki 2 a 2 = 0 


a i a 2 



Normalization: 


21 m kj a j a k = &rs * 

jf k 


r = s = 1 : 


r = s = 2 : 


M (a^) 2 + M (a^) 2 
M (a^) 2 + M (a^) 2 


overall sign 
undetermined 

l 



t 


same comment 


Step @: 

r 

Xl(t) = 7m Re ^ ei<ot + p2 ei< ° 2 ^' 

< 

X 2 (t) - -JL Re (-13, e iWlt + |3 2 e 1 "-*). 
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[Note we have ( n r (t) = |3 r e lc ° rt ) 

xi = 2 a^n r = ajn 1 + a*n 2 , 



Physical interpretation of the modes follows from these results. Set = 0 => x 1 = x 2 
for mode 2 (“symmetrical mode”)- Set n 2 = 0 => x 1 = -x 2 (“antisymmetrical mode”).] 
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Step ( 5 ) : As an example, consider the initial condition, 

X], (0) = -x 2 (0) = A, x^O) = x 2 (0) = 0. 

Then 

real part 

of Pi,P2 

1 1 

A = V2M ^ lr + ^ 

“ A = V2M ^ Plr + ^ 

=> P 2r = 0 , |3 lr = V2M A 

Also 


imaginary part 


Xi(0) -J~2M ^ 1 ^ li ^ 2 ^ 2i ) ® r 

x 2 (0) = ^ 1 ^ li _ ^ 2 ^ 2i ) = ® r 


Pli — P 2 i — 0 . 


If we had chosen x 1 (0) = -x 2 (0) 0, then we would have had P 2i = 0, P H =-x(0) , 

co 1 

and mode 1 would still be the only one excited. Any combination of the initial conditions 
x 1 ( 0 )= - x x ( 0 ) =- x 2 ( 0 ) 0 and/or x^O) = -x 2 (0) ^ 0 would have also excited the anti- 

symmetrical mode alone. 

Complete solution for the above boundary conditions: 

Xj = A cos (oc^t), x 2 = -A cos (oc^t). 

The initial conditions x x ( 0 ) = x 2 ( 0 ) 5 * 0 and/or x^O) = x 2 (0) 5 * 0 would have excited 
the symmetrical mode ( 2 ) alone. The general motion is a linear combination of the two modes. 
It is complicated to understand in general, but simplifies in the cases k 12 << k (“weak coupling”) 
and k 12 >> k (“strong coupling”). Lets examine these cases in more detail. 
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12.4 WEAK/STRONG COUPLING 

Let s rewrite the above general solution to prepare to discuss weak coupling. Introduce 


CD, + cd 0 co, - co 9 

CD 0 - 1 , CD b - - L -- 1 


t 

"beats'' 

Can show (mucho algebra!) that the above general solutions for x (t), x (t) may be written as 
(solve for co 1? co 2 in terms of cd q , co b , substitute above and use trigonometric identities), 

A A 

x,(t) = , 1 cos(co 0 t) cos(oyt)- iJ= cos(co 0 t) sin(co b t) 

V2M V2M 

A A 

+ , 3 sin(co 0 t) cos(co b t)- -jJ= sin(co 0 t) sin(oo b t) , 

V2M V2M 


A A 

x 2 (t) =-cos (co 0 t) cos (co b t)- -j== cos(co 0 t) sin(co b t) 

V2M V2M 

* A 

sin(co 0 t) cos(co b t) + , 1 sin(co 0 t) sin(co b t) . 


V2M 


V2M 


where 

Ai — Pir ^2r f A 2 = Pli — P2if 
A3 — — Pli — P2i f A4 = Pir — P2r 

(Please confirm this.) Notice that there are still only 4 undetermined (real) constants, A 1? A 2 , A 3 , 
A 4 . Let us now plug in the initial conditions, 

Xi(0) = D, x 2 (0) = 0, x^O) = x 2 (0) = 0. 


These initial conditions excite both modes 1 and 2, as opposed to the previous set. Plugging in 
above, we find 

A 1 — ■%/2M D r A 2 = A 3 = A 4 = 0 • 

So the general solution is 

x x (t) = D cos (co 0 t) cos (co b t) , 
x 2 (t) = D sin(co 0 t) sin(co b t) . 
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This form is convenient for discussing the case of weak coupling, k 12 << 
^i 2 we have 

C — r 

2k 

small compared 
to OOq 


1 



Of course 



is the decoupled frequency. Motion looks like: 


T b 


x,(t) 
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XJt) 


T 

^ = 10 

o 



H 


2jc 2jc 

T o = — ' T b = — 

<x>„ C0 U 


? T b » T 0 i n W eak coupling. 
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Phenomenon of beats occurs anytime a linear combination of very close frequencies is present. 
(Listen for them the next time you take a ride on a two-propellar airplane.) This corrresponds to 
an approximate multiplicative modulation of the uncoupled motion. 

To discuss strong coupling, let us recast the above solution (still exact for the given BC s) as 


also confirm 


x x (t) = — [cos oc^t + cos oo 2 t] , 

D r i 

x 2 (t) = — [- cos oc^t + cos co 2 tJ 


Then, for strong coupling, k u » k, we have 


/ 

k ] 

1/ \ 
72 

8 s 



\ 

l 2k i 2 J 

/ 


co 1 - 


2k 12 


H M 

2 ) 


, C0 2 - 8C0 1 


We find the motion looks like: 


x x (t) 


^ = 10 

C0 2 
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X (t) ^1=10 

ocl 


D 

D 

2 


— cosovt 
2 2 



Now get an additive modulation. The two frequencies here are understandable from: 
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12.5 EXAMPLE USING MECHANICAL/ELECTRICAL ANALOGY 

Lets do one more example. It is a coupled system we encountered before, but we did not solve it 
because the equations were coupled. We had in Ch.3: 


/ 

/ 

/ 

/ 



///////// 


x x : length of spring 1 
x 2 : length of spring 2 


(These lengths are compared to the equilibrium lengths.) 


m 2 (x, 


= -k 1 x 1 + k 2 x 2 

© 

+ x 2 ) = -k 2 x 2 

© 


(no wall 
on this 
side) 
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Reminder of the mechanical/electrical analogy: 


1 


F V , 

x I , k — 

C 

x «• q , 

m L , 

© => 

* * 1 1 

L iqi + — q x - — q 2 = 0 # 

© =• 

l 2 (% + q 2 ) + -r q 2 = 0 • 
^2 

Circuit looks like: 



L, 



Assume (do not need to know T or U here since we already have the equations of motion) 

, , v icot 

q lr2 (t) = a lf2 e 


-L.oo 

C, 


-L 2 o) a 1 + 


a i - — a 2 = 0. 


- L 2 cd 


a 0 = 0. 


det 


— - L x C 0 2 

Cx 


1 

C, 


-L 9 co 2 — - L 9 o ) 2 


= 0, 


1 


- L 1 0 ) 


- L 2 oy 


L 2 co 2 = 0 , 


or L 1 L 2 co- 1 ^ 0 )-L 2 ol)-L 2 o) + 


1 - 2 ■ = o 


c c c 
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A quadratic equation in go 2 . Roots are: 


1 1 1 

H-i 1 


1 1 \ 

.^ 2^2 ^ 2^1 ^ 1 ^ 1 . 

- 2^1 

\ C 2 L 2 

^ 2^1 -^ 1^1 ) 


4 


Because the above is so complicated, let us just look at a special case: C 1 = C 2 , L l = L 2 . Then 


GO 


1 

LC 



Call 


GO, 


3 + V5 
2LC 


G0 o 


3 - V5 
2LC 


These are the normal mode frequencies of this circuit. We will skip the step of normalizing the 
eigenvectors. The general solution is (real part understood) 


qi (t) = ay %t + 

q 2 (t) = a^e iWlt + a*e ic ° 2t 


8 real constants 


v i ; 1 2 2 

where the a. are in general complex. To find the eigenvector relations, substitute go^, go 2 back 
into the eigenvector equations: 


(1 

\c 




= 0 , 


a, = - 


1 + V5 


Likewise 


a 2 = - 


1 - V 5 
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We now have, 

q x (t) = aje 103lt + a^e 

q 2 (t) = - 


2 ^icovt 


+ S) 

1 icoit 

Un 

1 

\ —1 

2 J 

CL-^ C 

2 ) 


a* e 1 " 2 ' 


4 real constants. 


Using n 12 = e lc ° 1,2t as the normal mode variables, we find 


n i = 


V5 - 1 


<3i - q 2 


V 5 a\ ' 


n 2 = 


1 + V5 


\ 


qi + q 2 


V5 a* ' 


Mode 2 occurs when n, = 0, so 


q 2 = 


s/5-1 


ll * 
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Mode 1 occurs when 


q 2 = - 


i + 


q, 


These are like the modes for the coupled masses we saw in our first example, except the oscillation 
amplitudes are unsymmetrical. We can now build in the initial conditions. Let’s say 

<3i(0) = q 2 (0) = 0 , q,(0) = q 10 , q 2 (0) = q 20 . 

Plugging in above, we get (remember the a\o( 1 , 1 ), a\o( ] , 2 ) are complex) 


0 = Re [ito 1 a!‘ + ico 2 a?]. 


0 = Re 


1 + VH 


ico 1 a 1 - 


1 - S] 


ico 2 a 1 


Likewise 


Li), - ( a i), - °- 

t t 

imaginary part 


real part 

l l 

q io = (<l + HL' 


q 20 = - 


1 + JE) / jX (1 - sfe 




h‘). - i 


Vs - 1 

-q 20 + —z— <3i 


q 2 o + 


1 + Vs 


Full solution: 


qdt) = 


V5 


Vs - 1 


cos + 


V5 


1 + V5 

+ —-— qi 


cos oo 2 t, 


q 2 (t) = 1 = 


1 

V? 


1 + Vs 
+ i q 20 


COS + —j= 

V 5 


+ 


Vs - 1 


cos co 2 t 


Whew! 
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12.6 PROBLEMS 

1. Prove that the squared characteristic angular frequencies, co^, defined from (see Eq.(12.16); 
no sum on r) 


2 . 




cn 


— ^3 


2 


i,: 


r 


are positive, given the equilibrium conditions in the text. 

A double pendulum system is arranged such that at equilibrium the pendulums from 
which the masses, m, are hung are displaced at an angle, 0 O , from vertical as shown. The 
pendulum lengths are l and the spring constant is k. The unstretched length of the spring 
is L and X is the distance between the attachment points, as shown. 



Find the eigenfrequencies of the system for small oscillations about 0 Q . [Hints: First, show that 

the equilibrium angle, 0 () , is determined by Y = X-L = 2^sin0 o + tan 0 O . Then, 

k 

argue that the potential, U, is given for angles 0 1? 0 2 by U(0 1? 0 2 ) = mg^(l-cos(0O+01)) + mg^ 
(1 - cos(0 Q -0 2 ))+ — k(Y - ^ ? sin(0 o +0 1 ) - ^sin(0 o -0 2 )) 2 . Expand for small angles and solve for the 
eigenfrequencies. ] 
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3. A mass M moves horizontally along a smooth rail. A spring with spring constant, K, 
attaches the mass M to the wall. Let x be the distance that the mass M is located from its 
equilibrium position. 


eq. position 



Kw L ' / 
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The Lagrangian is (you do not have to derive this) 

L « — (M + m)x 2 + — (b 2 0 2 + 2bx0) - 0 2 - — x 2 , 

2 2 2 2 

in terms of x and 0 (0 << 1). 

a) Find the squared eigenfrequencies ( O) 2 ) of the system. 

b) Find the conditions on x and 0 which excite each of these modes. 

4. In the first semester we considered a double pendulum, consisting of two equal masses 
connected to each other and a horizontal support 



by weightless rods of length £. For small oscillations, the equations of motion we found 
were, 


6 + i 6 + £ e. 
2 C 


= 0 , 


0 , + 0 , + — 0 , = 0 . 

2 1 t 2 

a) Find the characteristic frequencies of the system. 

b) Solve for the normal coordinates, n, in terms of 0 , and 0 2 . Describe the conditions 
which excite these modes. 
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5. Three masses, arrayed as shown, are coupled together in a straight line with two springs, 
both with spring constant, k. This is a one-dimensional problem, so motion can only occur 
along the x-direction. 


k k 

• 'VW » A/W ♦ 

m 2 m i m 2 


x-direction 


a) Find the squared eigenfrequencies ( 0)^ ) of the system. 

b) Find the corresponding eigenvectors (they need not be normalized). Descibe the 
motion associated with each of the normal modes. 


Other Problems 

6. In two dimensions, a particle of mass m near the origin experiences the potential, 

1 2 1 2 

U(x, y) = — k x x + - k y y + k xy xy . 

Given that k xy < k x k y , and that kx and ky are positive, 


a) Find the eigenfrequencies of the system, co^ 2 - 

b) Show that 


2k 


( (K ~ k y) + >x - k y ) 2 + 4k^ y ) , 


excites one mode (which one?), and 


Y 


x 


2k 


xy 


UK 


excites the other. 


M - k y ) ! + 4k^), 
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7. Prob. 3.8 of this text considers a coupled system of two masses, as shown. 


/ 

/ 

/ 

/ 


kl 

m 2 

k 2 

—w— 

—mrv| 




/ 

/ 

/ 

/ 


The friction between masses 1 and 2 provides a coupling of the motion proportional to the 
relative velocity between the masses, x x - x 2 . Find the equations of motion and discuss 
coupled oscillations of the system. (There are actually no normal modes here since the 
friction dissipates the energy of the system. Look for complex characteristic frequencies.) 


8. Prove that Eq.(12.19), 

2 m Dk a j a k = 0 r ( r * S) 


holds for non-degenerate eigenfrequencies, co^ ^ C0g. [Hint: Look at the similar proof 
in Ch.l 1.] 
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COUPLED OSCILLATIONS 


9. 


Consider the hoop-mass system with angles 0 and ^ as shown. Both the particle and 
hoop have mass m; the radius of the hoop is R. 



a) Show that for small oscillations, the Lagrangian is 

1 . 3 ^ 

L « mR 2 (30 2 + - ii) 2 + 20110 - mgR(- 0 2 + - ip 2 + 0i|0 , 

2 2 2 

and the equations of motion are: 

— (cj) + 0) + (2cj) + 0) = Of 
a 

— (3<|> + 0) + (6(|) + 20) = 0 . 
a 

b) Find the normal mode angular frequencies of the system, oo i and co 2 . 

c) Find the time-independent ratio, (|)(t)/0(t), for each of the normal modes of this system. 
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13 SPECIAL RELATIVITY 

13.1 INVARIANCE AND COVARIANCE 

Newtons equations are: (a = 1 , , n) 


in 


dt 2 


= f„ 


(3n equations) 


(13.1) 


This is for an n-particle system. Showed under certain circumstances they could be written as 


dL 

dx . 


_d_ f dL 

dt 1 dx . 


= 0, 


where L = T - U. We then have 


(13.2) 


t = \ 2 m « x «' 


(13.3) 


u = I u «p 


a < |3 


X a " X p 


)• 


- 

f 


= -v a u aP , 

= E 

|3*a 


(13.4) 


Let us try making the transformation, 


a a 


x„ + 6x . 


(13.5) 


in L. Describes a change of origin, a translation. Under this change 


T -» Tl 
U -*• u| 


Lagrangian is said to be invariant. We also learned in Ch.6 that this means something is conserved, 
in this case linear momentum. On the other hand, consider the transformation 


x - t6v. 


(13.6) 
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in L. Describes a change in velocity or inertial coordinate system. Called a Galilean transformation 
or a Galilean boost. Under this transformation, 


T 




U -> U. 


Lagrangian is not invariant => nothing conserved. However, the equations of motion derived from 
the new Lagrangian are unchanged. We still get 


in 


dt 2 


/ 


in the boosted frame. The equations themselves are said to be covariant. Def-s: 


Invariant: 


unchanged in value as the result of some transformations. 


Covariant: 


Equation or quantity that is unchanged in form as the result of a transformation. 
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13.2 TWO POSTULATES OF SPECIAL RELATIVITY 

Physics at the end of the 19— century based upon 3 things: 

1. Newtons equations (mechanical phenomena, including gravity and waves.) 

2. Maxwells equations (electromagnetism) 

3. Galilean transformations 

Although we still accept Maxwells equations today in the same form as in the 19— century, their 
interpretation is completely different. Through the influence of the mechanical view of the universe 
from 1. above and the knowledge that light had wave characteristics (diffraction, interference, 
etc.), it was assumed that there was a medium for its propagation, the ether, and that the form 
of the electromagnetic equations given by Maxwell were valid only for a frame of reference at rest 
with respect to the ether since the equations were not covariant under Galilean transformations. It 
was thought that reference frames in uniform motion with respect to each other were completely 
equivalent as far as mechanical properties were concerned, but not with respect to electromagnetic 
phenomena. There was only one frame of reference for which, for example, it was thought that the 
velocity of light had the value predicted by the “rest form” of Maxwells equations. 

The Michelson-Morley experiment (M887) was designed to test for motion of the Earth relative 
to the ether. Of course it failed to detect any. So experimentally: 

i. The speed of light is an invariant under velocity boosts. (That is, it is observed to 
propagate with the same speed by all inertial observers.) 

Einstein accepted this fact even though he was, perhaps, not up to date on the M/M experiment. 
He also realized that the concept of the ether was superfluous. With no ether rest frame, the only 
frame of reference that can have any significance to an observer is the frame fixed to himself or 
herself. Therefore, Einsteins second postulate: 

ii) The laws of physical phenomena are covariant under velocity boosts. (That is, the form 
of the laws are the same for all inertial observers.) 

Einstein thought of “physical phenomena” as either electromagnetism or gravity. This has now been 
enlarged to include other forces; “weak” and “strong” nucleus forces as well as electromagnetism and 
gravity, (ii) above implies modifying Maxwells equations or abandoning Galilean transformations. 
It was the boldness of Einstein which led him to modify the transformation laws between 
inertial observers. 
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13.3 LORENTZ TRANFORMATIONS DEDUCED 

Consider: 



A light pulse is emitted from the common origin of K, K’ when they coincide at a singli 
of time. According to (i), the wavefront is described by: 

K: r = ct => r 2 - c 2 1 2 = 0 , 

K' : r' = ct' => r' 2 - c 2 t' 2 = 0 . 


By definition: 


f 2 


— 2 —2 —2 

x 2 + x 2 + x 3 , 


3 

-2 


X; 


f' 2 


X 1 + X 2 + X 3 


3 

= 2 


X; 


(The bar means coordinates relating to the position of the light pulse.) Or: 

K: ^ x ± 2 - c 2 1 2 = 0 , 

i 

K' : J ^l 2 “ ° 2t ' 2 = 0 • 
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Introduce 


= ict 

, = ict' . 



4 


K: 

O 

II 

" A 

IX 

(13.11) 


[i = l 



4 


K' : 

t K 2 - 0 • 

(13.12) 


|A = 1 


Usually use i,j,k = 1,2,3 (Latin) but |u, v, A, = 1,2,3,4 (Greek), is called a 4-vector. Because 
the form of the equation for the light pulse front has the same form in K or K’, we realize this 
relation is covariant under a change of inertial observers. (Although the values of the individual 
terms change in going from K to K’, the equation itself is unchanged in form.) We want the most 
general (linear) transformation (replacing the Galiliean one) which is consistent with this relation, 
but which reduces to the known Galilean transformation for low velocities (compared to light). (It 

is easy to show that ^ x^ 2 =0 is not covariant under the above Galilean transformation.) 

n-1 
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Thus, we assume the transformation law has the general form of: 


X 1 = > 4 l X l + ^ 12 X 2 + ^ 13 X 3 + f 

X 2 = ^' 21 X 1 ^- 22 X 2 ^ 23 X 3 ^ 24 X 4 r 

X 3 = ^31 X 1 + ^32 X 2 + ^33 X 3 + ^34 X 4 A 

X4 A.44X2 + A.42^2 ^ 43 X 3 ^* 44 X 4 f 


(13.13) 


where the X s are not yet known. (Notice we have dropped the bars on the coordinates - we are 

4 

not necessarily talking about light pulse coordinates at this point, so we do not have ^ x^ 2 = 0, in 

[A = l i 

general.) The final form of the above will be called a Lorentz transformation. Can think of xm, 
as representing the space and time coordinates of an arbitrary space-time “event” from the point 
of view of two different inertial coordinate systems, the common origin of which is at the place/ 
time where/when the K, K’ spacial origins coincide. 


The general form of this transformation looks similar to the orthogonal transformations of Ch. 1. 
Actually, orthogonal transformations are consistent with the covariance of ^ x^ 2 = 0. We had 

= 2 ^ ' 

< i 

t' = t . 

It preserves the “length” of the space, time terms in the sum individually: 



t' 2 = t 2 . 


So, 3-D orthogonal transformations are a specific type of Lorentz transformation. Realizing this, 
we now think of Lorentz transformations as describing velocity boosts and/or rotations. It also 
describes more than this, as we will see. (A Poincaire transformation also describes displacements.) 


for a L.T.: 


x 


[i 


2 V x v 


matrix : 
x' = Xx t 


Therefore 


2< 2 = 2 • 

M- ^V,Y 


(13.14) 


(13.15) 
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We want to transform the statement ^ x^ 2 = 0 into ^ x^ 2 = 0 . The only way this can be done, 

M- H _ 

given the linear transformation in (13.13) above, is if the \ an< ^ A quantities are, in 

ii |i 

general, directly proportional to one another (see also prob. 13.18): 


unknown proportionality factor 
1 



t 

no directionality in 
the universe ( homogeneous 
and isotropic ) 


(13.16) 


Use (ii) above now. Must also have 


E x » - F 


x„ 2 - F 2 


2 


5 X 2 ' 


=^> F = ± 1. 


The minus sign (F = -1) is eliminated for two reasons: 

1. The transformation must be continuously connected to the identity transformation, 

x , = v 

2. Causality of events (seen later). 


Therefore ^ x^ 2 = ^ x ^ 2 and so 

ii ^ 

2 v\ y x v x y = 2<- 

[i 

Coefficients of individual x v x y terms must be equal: 

=> y x x = 6 . 

/ J |iy vy 


(13.17) 


(13.18) 
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Reminds us of (Ch.l) 

I - V 

i 

t 

notice 


We then wrote this as 


X T X=l. 


But 


X' 1 X = 1 => ?i T = k' 1 (defines an orthogonal transf.) 
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Then 


XX -1 = 1 ^ XX T = 1, 

- 2 V ^ - S,K • 

i 

t t 

notice 


Situation is essentially identical here, except /is are now 4-D quantities: 


2 =» - 1 . 

Same chain of reasoning leads to 


(13.19) 


% = X" 1 and 2 ^V|x\|X - 8 


vy 


(13.20) 


Go back to 3-D again. An explicit rotation about the 3 axis is given by 


x^ = x 3 cos 0 + x 2 sin 0 , 

x 2 = -x x sin 0 + x 2 cos 0 , 

x 3 = x 3 , (t' = t also, of course) 


(13.21) 


From which the X.. were identified as the direction cosines of the rotation. The same form must 

>) 

hold in 4-D, except that when we make the physical identification, x 4 = ict, some of the “rotation” 
parameters will turn out to be imaginary. So, “rotate” in the 1,4 plane: 


< 


x. 


x 4 cos ip + x 4 sin ip , 


-x 4 sin tp + x 4 cos tp. 


Tie this in with physics now. For the origin of K ' , x( = 0 , we must have Xj = vt: 

0 = vt cos vp + x 4 sin \p 

t 

ict 


=> v = -ic tan ip or tan xjj 


. v 
1 — . 
c 
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Since v is real => \|/ purely imaginary. Let \|/ = ia (a real) 


tan ip = tan (ia) = i tanh a 

sin xp = sin(ia) = i sinh a 
cos ip = cos(ia) = cosh (a) 



tanh a 


v 

c 


So, in terms of real parameters, this “rotation” is 


x^ = x 1 cosh a - ct sinh a , 

*2 = X 2 ' 

x ; = x 3 , 

ict' = -ix 1 sinh a + ict cosh a 


Then 


v 


sin xp = 


1 + cot xp . c 

' v 


c 2 


= -sinh a 


COS = 


1 + tan ip v 


2 ' 


" = cosh a 


v 


From the above, sin y must be imaginary =^> — < 1. 

c 


sinh a = 


v 

c 


i 1 -? 


cosh a = 


I 1 -? 
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We want this to reduce to 


X 1 

= X 1 

X 2 

= X 2/ 

X 3 

= X 3 r 

t' 

= t 


when — 
c 


« 1. Therefore 


sinh a 



cosh a = 


+ 1 
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So we have 


x; = y(x 1 - vt), 
x' = X 2 , 

X 3 = X 3 f 


t' 



A Lorentz boost 



(13.22) 


The form of this transformation was known before Einstein by Lorentz, although it was Einstein 
who gave it the correct physical interpretation. There is a lot more to physics than just getting the 
equations correct! 


This transformation was based upon requiring the invariance of 

2 K = 0 • 

This in turn implied 

2 K = 2 < 2 b °) • 


This means ^ is an invariant under a Lorentz transformation. Let’s call it “invariant squared 

\x 

distance” and give it a special symbol: 

s 2 - 2< (13.23) 

t 

(can be + or - ) 


There is no reference to the coordinates of a light pulse here; it is a completely general relation 
relating lengths and times in arbitrary inertial reference systems. It actually relates the length and 
time intervals between two arbitrary space-time events if we write 

A s 2 ■ 2 Ax u • (13.24) 

Then the first event is no longer required to happen at the common origin of the K, K’ systems. 
So, we have assumed one invariant, the speed of light, 

c = c '' (13.25) 
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and derived another, the squared invariant distance, 


As 2 = As' 


(13.26) 


Thus, everything is not relative in relavity; there are still absolute quantities. 

13.4 ALTERNATE NOTATION FOR LORENTZ TRANSFORMATIONS 

One often sees an alternate treatment of Lorentz transformations with no reference to imaginary 
numbers. It is often used in particle physics. Define 2 

x° = ct. (no i!) (13.27) 


Then a Lorentz transformation is given as: 


notice: upper index. 

1 

x"* = 2 A\ x v . (|i, v = 0, 1, 2, 3) (13.28) 

V 

t 

different from previous A/ s 
(they are real ) 


We identify from above, 


~~ ^ij' ^ i “ - ^14' ^0 - ^44* 


Now 


all important minus sign missing! 


* 


EM 2 - 2( x ‘) 2 + cV • 

\i l 


Not an invariant! Define (“metric tensor”) 


g 


0 12 3 




(13.29) 
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9jtv: gn - g 2 2 _ 533 “ If 5oo — “1, all other components = 0. 
Then 


notice 

i 

2 x V x- - 2 ( x T - cV • 

fl,V 1 

This is an invariant. Sometimes define 


x „ = 2 




Notice: 


= (ct, x), "contravariant 4-vector' 
x^ = (-ct, x), "covariant 4-vector". 


(13.30) 


(13.31) 


Have to be careful in this notation to recognize whether an index is up or down. Then we may write 
s 2 = 2x'^ v x' v = 2X9^ • (13.32) 


|1,V 


\x,v 
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or, using the transformation law for x ^: 

=> A ^ AVkxXxK = 2 g ^ xXxK ' 

ji,v,y,K y,K 

^ A ^ A k = g^ K • ( 13 . 33 ) 

jl,V 

This is the analog, in this new notation, of the equations satisfied by the s (see Eqs.(13.20) 
above) and can be considered an alternate definition of a Lorentz transformation. 


Einstein summation convention: sum on any repeated upper, lower indices. (As we saw earlier a 
repeated index like x v x v does not produce an invariant quantity; this property generalizes.) 


The above may now be written: 

9hv a ^ aV k = 


(13.34) 


There is also a matrix interpretation of this form. Write (the “T” means transpose; notice the matrix 
version can violate the Einstein summation convention) 

(a T )\9[xv AV k = g X K ' (13.35) 

=> A T gA = g. (Takes the place (13.36) 

of A7 A,= 1.) 


No matter which notation and or conventions we use, it is important to realize the point of either 
is to build in the minus sign in front of the time components. We will primarily use the A^ -type 
notation in which all quantities are real. Famous quote: The start of any calculation is to check 
the authors conventions. 

In order to be able to use this notation, must know the rules. Just as for 3-D rotations, 4-vectors 
are defined by their transformation laws. We have 

a'" = 2 or = A\A\ (13.37) 

X 

where A /\ satisfies A T g A = g. This makes Ap a contravariant 4-vector. To make a covariant 
4-vector out of say, just “contract” with g^v : 

A , s g,v AV - (13.38) 
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Define: 


new quantity 

1 

- % , (13.39) 

t 

Kronecker delta 


=> g^ v A v . (13.40) 

|J,V 

Note from (13.39) that the 9" are just the inverses of the gpv, considered as matrix elements: 
g^ v g VK = => g _1 g = 1- The matrix representation of the g^ V in this case is also given by 

(13.29), which is the so-called flat space metric. 


Define A * = g^ v g^ K A v K ; then the transformation law for covariant 4-vectors is just: 


Al = g, K A' K = gji\^ = q^Al A v , 


*|xk X' 


[IK - 


A' = A V A 


(13.41) 


The matrix interpretation of these transformation laws is based upon the index to matrix 
correspondence, 

(The second statement follows from the definition A^ = q [Xv q XK A v K . ) Notice the matrix 
interpretation ignores the “up” or “down” position of the indices; I have placed these arbitrarily 
down. This implies the equivalence of the statements, 

A'^ = A^ V A V o A ' = AA, (A' a column matix) 

A[x = A pi v A v <^> A ,T = A T (A ) _1 . (A ,t a row matrix) 


Likewise (see prob. 13.2(a)) 

A\A/ = 5^ a|(a -1 ) j = 1 or AA -1 = 1. 
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A way to form a Lorentz invariant quantity, or in short, a scalar , is to form (just like we did for 
x*\ x^) , 


= 


A'^ v A' v = 


g^A; a v 


J/.K 


A'^A' = AfA,, . 


(13.42) 


Of course s 2 — x^x^ is also a scalar. We have the first of our “index jockey” rules: 


1. To make a Lorentz scalar out of two vectors, sum on one upper, one lower index. 
We know that the components of the gradient operator, 


W - 


d 


dx 1 


(13.43) 
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transform as a vector. Can then build vectors (V (|), for § a scalar) or scalars 
( V xA for A a 3-vector) for A a 3-vector) out of it. For Lorentz transformations, 
define 


V 


\x 



t 

lower 


T 

upper 



(13.44) 


Can show it transforms as a covariant 4-vector (thus, the lower index). Also 


, V 1 " = g M ' v V v => = 




5x 


Id- 


C dt 


,V 


t 

upper 


t 

lower 


(13.43) 


transforms as a contravariant 4-vector. So, the second index jockey rule is: 

2. transforms as a contravariant 4-vector 
transforms as a covariant 4-vector. 

Two other rules are: 

v 

3. Lower an index with P|dv : P|dvA = 

4. Raise an index with P^ • g^ A v = A^. 


As before: P^ Pvk = 


( $0 cl c2 c3 -i \ 

Oq — Oi — o 2 — o 3 — 1 , 


all others = 0 


Matrix statement: (g' ! )g = L 


? [.IV 


1 

2 

3 


0 12 3 


-1 


(=> V 11 = g^ v V v ) 
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Note: 

* Using 9|av f raises or lowers an index but never changes index order. 

* Number of free up, down indices on one side of an equation must equal the number on 
other side. 


Exampl es: 

1. Show At- 1 Bt- 1 : At 1 B^ = A CT g°P B T g x ^ = A a B T = Bt 1 . (In 

general A---t J ----B--- (X --- = A • ■ ■ ^ • B • • • t 1 • • • for arbitrary tensors A • • • t - 1 • ■ • 
and B- • -p- ■ ■.) 


2 . Given A' !I = A(‘A V ,%' = %, prove that transforms as a contravariant 4-vector: 

dx 1 

dA ,[l Ali dA v Ali dA^ 1 dx AW dA 11 

- Av - A -r = A v -• 

dx ' dx ’ dx dV dx 

T 


3. Prove that V^cJ)' , cp' = cp , transforms as a covariant 4-vector: 

need 

1 

V' ((>' = V (j) . 

M ax'" ax'" ax v ax'" v 

But 

(*"■ - *y) ■ g to A“ p , 

==• - g, c xl )' g l,; - 


a^a A a x ,A = x Y 

y yxa yv 6 A A r 


i'K 


A/ 


ax Y 

dx' x 


= AJ . (Also 


ax 


?A 


ax K 


a\> 


Therefore 


V>' = A v V v <() 


covariant 4-vector. 
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13.5 THE "LIGHT CONE" AND TACHYONS 

Let s get back to the physics being described. For this, we need a picture of the light-cone (imagine 
rotating about the time axis; the “surfaces” shown are hyperboloids): 


future 
time like 
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timelike” 

s 2 < 0 means 

|cAt| > |Ax |, 

spacelike” 

s 2 > 0 means 

|cAt| < |Ax |, 

lightlike: 

s 2 = 0 means 

|cAt| = |Ax |. 


This diagram describes the relationship of the origin to all other space/time points. Events with 
a timelike behavior can always be given a definite causal order since a light beam can reach the 
second event before it occurs. Also, |cAt| > |AX | is preserved by Lorentz transformations. This 
means that spacelike and timelike type events are not mixed. On the other hand, spacelike events 
are fundamentally different: there is no invariant order of these events seen by all observers. Thus: 

Timelike events have an invariant causal order 
Spacelike events do not have an invariant causal order 

Better not mix up such distinct types of events! This is one reason for discarding the Lorentz 
transformation with F = -1, considered earlier, because such a transformation on the light cone 
would take us from a timelike ordering of events to a spacelike one, or vice versa. Mathematically, 
this means we are prevented from jumping across the lightcone. 

To see this transformation more explicitly, consider (the “rotation” interpretation of this 1,4 
transformation no longer holds because of the factors of “i”): 

x 4 = ix-L sin ip -ix 4 cos ip , (13.46) 

x 4 = ix-L cos T + ^ x 4 s i n T • (13.47) 


If we insist as before that the origin of k ' , x' 

1 


sin ip = 


1 “ ^2 
c 


0 , has x 1? has then = v t, then tan op = i —5 and 


(13.48) 


sin \|/ must be chosen imaginary: 


> 1. 
c 


(13.49) 


Hypothetical particles with — > 1 are called tachyons. These particles can be shown to move 
either backward or forward in time, depending on the direction of motion! So by eliminating 
F = -1, we are eliminating the possibility of transformations to frames of reference with — > l. 
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13.6 MATHEMATICAL PROPERTIES OF LORENTZ TRANSFORMATIONS 

Let s do a little investigation of mathematical properties of Lorentz transformations. Simple facts 
about determinants: 

det AB = det A detB , det A T = det A. 

Since 

A T gA = g , (13.50) 

taking determinant of both sides gives, 

det^A T gAj = (det A) 2 det g = det g , 

=> det A = ± 1. (13.51) 

+1: "Proper" Lorentz transformation 

-1: "Improper" Lorentz transformation 

Will see the physical interpretation later. Also, we have 

x' = A,x, (13-52) 

x " ” A 2 x', (13.53) 

x" = A 2 A x x = A 21 x . (13.54) 

Is A 21 also a Lorentz transformation? Use the definition: 

A 2 g A 2 = g , (13.55) 

=> A 2 A 2 gA 2 A x = Aj'gA 1 = g, (13.56) 
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However, A 21 = (A 2 A x ) T = A*A 2 , so 

=> A 2 i g A 2i = g • (13.57) 

=> ^21 is also a Lorentz transformation 
How many independent elements of A^ v are there? Expect physically: 

3 (rotations) + 3 (boosts) = 6. 

Mathematically: , 4X4 = 16, real quantities. However, the equation 

A T gA = g , 

is also (like 3-D rotation) unchanged by the transpose. Therefore, there are as many components 
as in a symmetric 4x4 matrix: 

/ 

x x 

X 

\ 

This means there are 16 - 10 = 6 independent real quantities, the right number. 

Some terminology: 

Restricted Lorentz Transformation (RLT): Transformations physically possible for a frame of 
reference associated with a material observer. 

What this means is rotation + boosts (I will be more precise mathematically in a bit.) We have 
learned that there are Lorentz transformations with det A = +1 or -1. There is a further discontinuous 
possibility as follows. Take the 00 component of the defining equation: 

g„v A \ AV K = g *,k f x = k = o. 

= ^ > g 00 A o A o gn A o A o g 22 A o A o + g 33 A o A o — goo < (13.58) 


X X 
X X 

X X 
X 


10 
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9oo = 9ii = 1 / 


- - ( A '.) 2 + I ( A 1 .)' - -1 - 

i 

=» (A"„) 2 -1 + 2 (a 1 .) 2 * 1 • 

i 

So, there are two possibilities for A° 0 : 


Can also show : ' 

( A ”.) 2 = 1 + 2 ( A ° i ) 2 • 


(13.59) 


(13.60) 


A o * 1 


'Orthochronous' 


A o * "I 


'Nonorthochronous' 


Now our more precise mathematical definition of a restricted Lorentz transformation: 
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im¬ 
proper, orthochronous Lorentz transformation: 

det A = 1, A° 0 > 1 , A T gA = g . 

Prove: If A and A 2 are both RLT s then so is A = A A . 

Show det A 21 = 1 (easy): 

det A 2 i = det A 2 det A 1 = 1. 


Show that (A 21 )° = (A 2 )° (A 1 ) y > 1 (not so easy):From the above 

or (A 21 )\ - (A,)" (A,)’, 


^21 ^ 2^1 


( a 2 .)"„ - N\(\)\ - (A 2 )° 0 (A,) 1 . + 2(A 1 )\ ( A .f 


From the above 


((Ai )° 0 ) 2 - E((A,)‘.) 2 . K ) 1 * SN' 


“Cauchy inequality”: 


n n 


2^2 b ? 2 2 


k=l 1 =1 


l2 


a k b k 


k = l 


J2N , I )2((A.y.) • 


Cauchy inequality 

t 


SfA.rjA,) 1 , - 2 (A,)”, (A,) 1 


The sign of (A 2 >■ the same as ((a 2 )° 0 (A,) 0 ,). 
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Result: 

No matter how many RLT s are performed, the result is still a RLT. 

Another point: 

Consider two timelike space-time events: 



These events must have an invariant time order since they can have a cause and effect relationship. 
There can not be some observers who measure t > 0 and some who measure t < 0. That is, as a 
result of a RLT (boosts done on a material observer), we must show that the time component of a 
timelike 4-vector has an invariant sign (either + or -). (This will be done in a HW problem.) It is 
clear that spacelike events can have the sign of their time components changed by a RLT since there 
can be no causal relationship between them. As pointed out previously, no LT changes a timelike 
event into a spacelike one or vice versa. 

There are 4 possible combinations of the discontinuous parameters det A and A 0 0 . 
We classify these as 


1. Proper orthochronous (RLT s): det A = 1, A° 0 > 1. 


2. Improper orthochronous: det. A = -1, A° 0 > 1. 


3. Proper nonorthochronous: det A = -1, A° 0 < 

4. Improper nonorthochronous: det. A = -1, A 0 


- 1 . 

< - 1 . 
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We have already talked about 1. What do the other transformations imply? Possibility 2 describes 
space inversions in addition to the usual RLT s. Example: 


0 12 3 


a 2 = 


0 

1 

2 

3 


-1 


-1 


-1 


Gives: x -* -x, t 


t . Any other A with det A 


- 1 , 


A 0 


> 1 can then be written 


A — A rlt A 2 . 


( 13 . 61 ) 


Not new. Saw space inversions previously as part of 3-D orthogonal transformations. Takes us from 
one side to the other on the spacelike or timelike s 2 = const, surfaces in the lightcone diagram. 
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Possibility 4 describes something new: time reversal. Example: 

0 12 3 

' -i 

i 

i 

i 

i / 

Gives x —> x, t —> -t. Any other A in this category can be written as 

A = A rlt A 4 . (13.62) 

Does to the time axis what possibility 2 did to space. Like space inversion, this is not a physical 
operation we can imagine carrying out on a material observer (not an RLT). Never-the-less, time 
and space inversions are valuable concepts, especially in quantum mechanics. To begin to understand 
why, it is helpful to consider the collision of particles (representing atoms or smaller objects). For 
example, consider the collision of two particles A and B: 




If Newton s equations held at all levels of description, microscopically, as well as, macroscopically, 
it would be just as probable to see the time reserved reaction: 


time 


4 



velocities 
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The reason for this is because Newtons equations of motion, 

\ 

= 2 U “P (l*» " *p|) * 

V |3<a / 

are invariant under x a —> x a , t —> -t .Modern microscopic theories of matter are also time- 
reversal invariant to a great extent; however, small violations of this symmetry have been seen 
experimentally. The fact that time-reversed reactions occur essentially as often as the original given 
one is an important fact of nature which is reflected in statistical mechanics and thermodynamics. 
However, on a macroscopic scale we know that many processes are not time-reversible; a glass 
shattering on the ground is one example. This is not a contradiction but an aspect of the behavior 
of large numbers of particles whose motions must be treated in a statistical sense. See the page 
from T.T. Lees book “Particle Physics and Introduction of Field Theory” for a further illustration. 

Possibility 3 represents combined space and time reversals. Example: 

0 12 3 

' -i 

-i 

-i 

-i 

v / 



m„ 


a X o 

dt 2 


“ f(7 ' 


Notice 


a 3 = a 2 a 4 . 

This A gives x —> -x, t -» -t, and any other A with these characteristics can be written as usual, 
A = ^RLT -V. (13.63) 

13.7 CONSEQUENCES OF RELATIVITY 

Next, some consequences of special relativity. Consider again two space-time “events.” Invariant 
associated with them is 

As 2 = ^ (Ax ± ) 2 - (cAt) 2 . 

i 
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Divide by -c 2 : 


At 2 




(13.64) 


Ax (or just x): “proper time”. It is invariant and > 0 for timelike events. What is its physical 
meaning? When two events happen at the same spatial position for some observer, Ax 1 = 0 and 


Ax = At. 


Since the value Ax is invariant, this means we can interpret proper time as follows: 

Proper time: time interval measured by an observer who is at rest with respect to two events. 
(The two events take place at the same position.) 

This means 

Ax = At, "rest observer" 

t 
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same numerical value 


1 



any other inertial 
observer 


Conclusion: At’ > At. 

The rest observer always measures the shortest time interval between these events. The name for 
this effect on time intervals is time dilation , because from the point of view of the moving observer, 
the time interval has increased or dilated. Can get a quantitative measure of this effect from our 
Lorentz transformation boost along x : 


x i : 

= Y 

( X 1 

- Vt), 

X 2,3 


X 2,3 i 

r 



/ 

V 

t' = 

= Y 

t - 

x i — 



V 

c 


Reminder: 



Situation : 

First event at 0,0' 
when the axes 
coincide; then a 
new event at 0 

every At seconds 
(as measured by K) 
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Event Is 

t' 

= i 

t = 0 


Event 2s 

t; 

= Y 

0 ) 
t, - t i — 

O 

II 

\—1 
X 




\ C t 


Event 3: 

t; 

= yt 2 



Relation between time intervals: At’ = yAt. We know that At = Ax, so this gives 




At ' = yAxl- 


\ 


Y = 


•J^i 2 


(13.65) 


Time dilation is well established experimentally. 


Another consequence. Consider an object at rest in the K frame directed along x : 



The observer in K measures the length of the object (which is at rest in his frame) and finds it to 
be £. K\ by moving along its length, will measure it to be 


V = vAt\ (13.66) 

where At’ is the time interval K’ measures for passing the object. Now, the passing of O’ past the 
two ends of £ defines 2 space-time events, both of which take place at the same place, O’, according 
to K\ Thus, K’ is measuring the proper time interval between these events, 

At' = Ax. (13.67) 


Download free eBooks at bookboon.com 


242 











MODERN INTRODUCTORY MECHANICS PART II 


SPECIAL RELATIVITY 


Therefore from before, we know the time interval measured by K will be dilated: 


< 

II 

-p 

<1 

(13.68) 

=> £' = vAt = — At. 

(13.69) 

Y 

£ 

Of course, K has At = — , so 


V 


e -Li. tjrrp. 

(13.70) 

Y V 

=>£'<£ , “length contraction” 



Proper length: length of an object which is at rest in an intertial observers frame. 

Interesting point: how an object moving with respect to an observer looks like is quite a different 
matter. Although the above considerations suggest the objects will appear 
contracted, in reality because of the finite speed of light, objects can appear to 
be rotated or have other distortions in appearance. For example, consider a cube 
moving directly away or toward an observer. It would actually look something 
like (consider light pulses which arrive at the same time): 


Toward: 


shape 

seen 




earlier 

position 

later 

position 


Away: 


shape 

seen 



later 

position 

earlier 

position 
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13.8 VELOCITY ADDITION LAW 

Another consequence: velocity addition law is changed. After boost 1 by 


i = Yi ( x i " Vit), 


X 


^ 2,3 ^ 2,3 f 


V, 


t' = Vi t - Xj -j 


Boost it again: (v 2 ) 


K = y 2 K - v 2 t'), 


X 2 2 ^2 2 ^2 2 f 


t" = Y 2 

\ c / 


(13.71) 


(13.72) 
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In detail, this gives 


x i = Y2Y1 


x i = Y2Y1 


Make it look like a single boost: 


* I V 1 V 2 

x i = Y1Y2 1 1 + 

c / 


I 

v,v 9 ^ 

1 / \ , 

X 1 

1+ 2 

- ( V 1 + v 2 ) t 


( C J 


1 

{ v.v,) 


X 1 

1+ ! 

1 - K + v 2 ) t 


1 c ) 



X 1 - 



\ 


V 1 

+ v 2 

"h 

1 + 

ViV 2 

c 2 j 

L* 


t" = Y2Y1 


111 + Vlv2 


2 I - X 1 


( V 1 + v 2 ) 


Exactly the same form as a single boost if 


Yl2 = Y1Y2 1 + 


ViV 2 


V 12 = 


Vj + v 2 
1 + ^ 


new velocity addition law. 


Are they consistent? 


(13.73) 

(13.74) 


(13.75) 


(13.76) 


(13.77) 

(13.78) 


Y12 


1 _ Xu 

° 2 1 - i 


; just as we wanted) 


Vi + v 2 
1 + ^ 


V 


c / 


Exampl e: Vi = .95c, v 2 = .95c (both in the same direction): 
.95c + .95c 


V 12 = 


(.95c) 2 

1 + --^ 


= .99868c . 


Can not boost through v = c in a continuous fashion. This gives another perspective on the need 


to eliminate the earlier tachyonic F = -1 transformation, which required 


> r 
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13.9 MOMENTUM AND ENERGY UNITED 

Up to now, we have concentrated on the space-time aspects of relativity. However, relativity also 
has important implications for momentum and energy. 

Use one of the above results: 

chU 

- is a contravariant 4-vector 

dx 
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Describe the motion of a particle in space and time: 



dx 11 = ( dx°, dx ) 


ds 2 

By definition (remember dx 2 =-) , 


dT = - V " ( dx ) 2 + ( dx °) 2 


(13.79) 


Must be a timelike vector if it describes the motion of a material observer. It is: 

dx^ dx , ds 2 2 

dx dx dx 

Write the components in a more ordinary notation: 


(13.80) 


dx dt dx dt _ 

dx dx dt dx 



dx 

=> — = yv -* v . 

dx v<<c 


Time component: 


dx° dt 

- = c — = yc -* c . 

dx dx v<<c 


(13.81) 


(13.82) 

(13.83) 


(13.84) 
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Re-verification of timelike character: 


' dx j 

V dx / 


1 

f— T 

v 2 

1 

l dx J 

c 2 


f dX °l 

2 


dx J 

" !_ T 

c 2 

dx° 

V _ 

V 2 

v dx 

) 

1-^ 1 

2 


-c 2 . 


V 


dx^ — 

Since the space component of- reduces, in the v « c limit, to just v, a natural definition of 

i . . . . dx 

relativistic momentum is 


P 


M- 


dx" 

m-. 

dx 


(13.85) 


=> p = ymv , p° = ymc . 


Meaning of p°? Notice that p° > 0 no matter which Lorentz frame we are in. (The sign of the time 
component of a timelike vector is unchanged under Lorentz transformations, which was a HW 

problem.) A particles K.E. is also always > 0, so does _ = p°? Expand p° in a power series in T: 

c c 


P = 


me 2 + — mv 2 + 


(13.86) 


Forces on us the realization that there is an energy associated not with motion , but with 
mass. Identify 


E n o « 55 

u = me , rest energy 
E = p°c . “total energy” 

Einstein (“The Meaning of Relativity,” p.47): “Mass and energy are therefore essentially alike; they 
are only different expressions of the same thing.” Kinetic energy, T, is just the difference between 
total energy, E, and rest energy E Q : 

=>E = E 0 + T , T = E - E 0 = me 2 (y - 1) . (13.87) 
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Now 


u, 2 2 

P Pn = -me. 


so 


22 -2 E" 

-me = p- - 


E 2 = p 2 c 2 + m 2 c 4 , 


(13.88) 


or, choosing the positive root, 


= +4 


— 2 2 2 4 

pc + m c , 


/ , . \ 
old connection : 


E = 


p! 

2m 


(13.89) 
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to be consitent with p° expressions. Now we realize we have done for p and E what has already 
been done for x and ct: unify them: 


xX = x 2 - x u = x z - c z t z , 


-2 2,2 


li - 2 0 -2 

P Pu = P“P = P“ - 

11 c 2 


The particle of light is called the photon; the connection between it s energy and momentum is 
given by (13.89) with m=0: E = |p| c . It follows a “lightlike” trajectory (dx = 0 in (13.79)) and 
therefore its inertial reference frame speed is always c. 

13.10FOUR SHORT POINTS 

1. Transformation law for pfy same as for x^ 5 x^: 

P* = A\ P V r 

p; = a; P v . 


2. Conservation of both momentum and energy is now contained in the single statement: 


2 M 

incoming 

particles 


2 M 

outgoing 

particles 


3. Because of equivalence of mass and energy can express mass in units of 


energy 

c 2 


nass = 


i tic; v 

V 

t 

common unit in particle 
or nuclear physics 


(1 MeV = 1.602 x 10' 6 erg.) 


[Often, one drops the _L and masses are often quoted simply in MeV.] 
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4. Equivalence between mass and energy is verified countless ways. For example, the “binding 
energy” of multiparticle systems appears as a deficit from the constituent particle masses. 
For the deuteron, 


deuteron mass: 1875.7 


MeV 


c 


2 


proton : 938.3 Mev/c 2 
neutron : 939.6 Mev/c 2 


} 1877.9 


MeV 


(1875.7 - 1877.9) = -2.2 MeV/c 2 (the binding energy) 
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13.11 PROBLEMS 

1. Show that the 3-D wave equation (c=speed of light), 


/ 1 \ 

v 2 - JLJL 

c 2 at 1 , 


f = o, 


(f is a scalar) is covariant under the Lorentz transformation ( y = _-_ p = iXi ) 

s/i - P 2 c 


x 1 = y(x - vt), 


y' = 


V 


t' = y(t - x —) . 

c 

[Hint: Use the chain rule, 

d_ _ dx^_d_ dt' d 
dx dx dx' dx at' 


ax' at' 

evaluating -, -, etc. from the transformation.] 

dx dx 

2. Consider Lorentz transformations using index notation. Given 

g„v A>\ A\ - g lK , a ” 5 - A”, A', AJ - 


dx 


dx 


A 

k 


3x 




dx K 


show: 


(a) 

A\ 

A k 

iV u 



(b) 

g^K 

A/ 

> 

II 

g^v. 

(c) 

g XK 


II 

> 

< 

g" v 
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3. Consider Lorentz transformations using matrix notation. Given that 
A T g A = g, g T = g, gg = 1 (“1” is the unit matrix), show that 

(a) (A -1 ) T g A -1 = g, 

(b) A g A T = g . 

(c) Write (a) and (b) in component language. 

4. For the tachyonic transformation, 

= iXj sin ip -ix 4 cos ip , 

x 4 = ix 1 cos ip + ix 4 sin ip . 


sin 2 ip 



v 

c 


> 1 , 


(a) Show that x' x 2 + x 4 2 = - Xi 2 - X4 2 . 

(b) Demonstrate that under ip —* id , where a is a real parameter, these become 

= -x x sinh a -ix 4 cosh a , 

x 4 = ix x cosh a -x 4 sinh a . 


(c) Show ^0 ~ V - 1 • (Similar to A° 0 a 1 _ f or Lorentz boosts.) 

(d) For the above, find det(A) = ?. 


5. 


Given the relativistic acceleration 


du 

dx 


where 


u 


dx _ dt 
dx dx 


v 
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( v = is the ordinary velocity; |3 


v 


show that 


a = 


v + — v x (v x v) 
_c_ 

2 v 2 


(1 - p 2 ) 


6. (a) Show that the relativistic form of Newtons second law is ( F = -=- ) 

dt 


F = 


mv 


(1 - ( 3 2 ) 3/2 ' 

if v and v are co-linear, and 


b) F = 


mv 


(1 - |3 2 ) 1/2 

if v and v are perpendicular. [Hint: Prob. 13.5 above.] 


7. (a) A particle known as a muon is generated high in the Earths atmosphere with a speed of 

0.96c relative to the Earth. The muons average lifetime, measured at rest, is 2.2 x 10' 6 
sec. How far does such a muon travel through the Earths atmosphere before decaying? 
Briefly explain your calculation, (c = 3 x 10 8 meters/sec.) 


(b) A 30 yr. old physicist travels to a star, at rest with respect to the Earth, which is 10 
light years distant. (1 light year = distance light travels in 1 year = 9.47 x 10 15 meters.) 
His/Her rocket ship travels at 0.85 c relative to the Earth and star. How old will the 
physicist be on arriving at the star? Briefly explain again. 

8. An astronaut brings an atomic clock on board the International Space Station (ISS), which 
is in orbit around the Earth at about 18,000 mi/hr. This clock is synchronized with an 
atomic clock on the Earth before the flight. If the clock is aloft in the ISS a month (take 
this as 30 days, measured lets say from the Earths point of view), by how much time, in 
seconds, is the ISS clock advanced or delayed (which?) when compared to the Earth clock? 
(The speed of light is about 186,000 mi/sec.) 

9. Show, by using conservation of energy and momentum, that a photon (a massless particle 
of light) with momentum P can not decay into two massive particles with momenta P a 
and Pb . (For simplicity, you may assume the masses of particles a and b are the same, 
although the result holds in general.) Extra: Can a photon decay into one massive and one 
massless particle? 
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10. Show that the time component of a timelike 4-vector t^ 1 , retains its sign under a restricted 
Lorentz transformation. [Hint: Combine the timelike condition, (ct 0 ) 2 > ^ (x 1 ) 2 , 
with a condition for a proper Lorentz transformation using the Cauchy inequality.] 

du^ 

11. Define u 1 * - — , a 1 * - —. Show: 

dx dx 

a) = 0. 

(i.e., the relativistic acceleration 4-vector is perpendicular to the relativistic velocity.) 

b) is a space like 4-vector (assuming oC ^ 0 ) . 

[Way 1: Use (a) and construct a proof by contradiction by assuming ah is timelike or 
lightlike. Way 2: (Brute force) Construct explicitly and show that > 0 . ] 
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12. Given an elastic collision a + b —»■ c + d. 



(Each line represents 
a particle with a mass 
m a , m b , m c , or m d ) 


show that the variables 

S ■ (Pa + Pb) (Pa + Pb) M 

t « (Pc “ Pa) (Pc " Pa) (j 

U = (Pd - Pa) (Pd - Pa) M 


satisfy the relation, 

« . ■ 2/2, 2, 2, 2v 

s + t + u = -c (m a + m b + m c + m d ) . 


[These relatistic invariants are useful in describing the results of scattering events in 
different inertial frames.] 

Other Problems 

13. Show V|uA^ transforms as a scalar. 

14. A particle of mass m is accelerated from rest with a constant force, F, in a straight line. 


v 


a) Show that if the final velocity is v ((3 = — ) , then the distance traveled, D, is given by 

c 


D = 


me 




- 1 


b) The acceleration takes place over a time period, T. Show that this time period is given by 

m me 2 [~ FD ~ 2 7 

cT = - J 1 + — j) ~ 1 • 

F V me 2 


[Hint: integrate the result of part (a).] 
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15. Consider a relativistic particle of mass m and charge q in a circular orbit of radius R. The 
force on the charge is central and is given by 



where the unit vector r points from the force center to the charge. Find the relativistic 

formula giving (3 = ixi of the particle as a function of m, q, and R. 

c 

16. A relativistic particle of mass m and charge q is moving with velocity v in a plane 
perpendicular to a uniform magnetic field, B . 



Using the Lorentz force law, find the radius of the orbit in terms of v, q, B and m. In 

E 

particular, show that the radius for ultrarelativistic motion (E » me 2 ) is R = r —, where 

qB 

E is the total energy (rest mass plus kinetic) of the particle. 

17. A coordinate system K’ moves with velocity v relative to another system, K. In K’ a particle 
has a velocity u' and an acceleration a' , whereas in K these quantities are u and a . Show 
that in K the observed component of acceleration of the particle along v, a ||, is given by 
Ivl 

( P = — ) r 

C 


an = = 


(i - m 


2 v 3/2 


(1 + 


V 


u\ 


a'i 


[Hint: First get the relation between u y and u' y . ] 
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18. There are actually three classes of transformations consistent with the invariance of the 
statement (|Lt = 1,2,3,4 notation; x^ = ict' ) , 

2 - 0 ' (D 

for the coordinates, , of a lightfront. These are: 

1 • < = 2 V x v ' where J = 6 vy , 

p. p 

II. , ("dilations") 

x„ + x 2 c„ . 

III . x„ = ----—— , c„ arbitrary. ((c • x) s > ex,, , 

! (1 + 2(c • x) + x 2 c 2 ) ! L, * * 

X 2 ^ 2 x^x^ , and where (1 + 2(c • x) + x 2 c 2 ) ^ 0 ) ( "conformal 

transformations") 
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Show that cases II and III are consistent with the invariance of equation (1). (Case I 
already considered in the text.) Why are Cases II and III not considered further in the 
text for relativity? 


19. In prob. 13.9, we saw that a photon can not decay into 2 massive particles (or into 1 
massless and one massive particle). Now examine two other cases: 

(i) m A —> y + y (particle of mass mA decays into two photons) 

(ii) ni A —> m B + y (particle of mass mA decays into particle of mass m B < m A and a 
single photon) 


Do energy and momentum conservation allow these processes to occur? If not, show it. 
If allowed, solve for the magnitude of photon momentum, p = E/c, in terms of m A in 
(i) and in terms of mA and m B in case (ii). (Assume particle m A is at rest in both cases.) 


20. Consider so-called Compton scattering on an initially stationary electron of mass m 
(use p = Ipl , p' = Ipl , P e = |P e | ) : 



(a) Given the energy, E, and momentum, p, of the incoming photon along the 
x-direction, write all the equations which follow from momentum and energy 
conservation. 

(b) From these equations show that 


/ 

-me 

V 



1 - cos 0 . 


This determines the momentum, p’, of the scattered photon in terms of the initial 
momentum, p, and the scattering angle, 0. [Hint: Get two expressions for the square of 
the electron momentum, P e , and set them equal.] 
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ENDNOTES 


1. Notice that we could replace condition 3 with the alternate condition that the system of equations, 
(12.17), for small oscillations yield only real nonzero eigenfrequencies, to r . This was the same 
condition for stability seen in Ch.2. 

2. Note: We are using “x” here to denote coordinates, but of course in the previous treatment 
jLt = 1,2,3,4 and x 4 = ict, whereas here \i = 0,1,2,3 and x Q = ct (x = -ct). 
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